-
An "ultimate" coupled cluster method based entirely on $T_2$
Authors:
Zachary W. Windom,
Ajith Perera,
Rodney J. Bartlett
Abstract:
Electronic structure methods built around double-electron excitations have a rich history in quantum chemistry. However, it seems to be the case that such methods are only suitable in particular situations and are not naturally equipped to simultaneously handle the variety of electron correlations that might be present in chemical systems. To this end, the current work seeks a computationally effi…
▽ More
Electronic structure methods built around double-electron excitations have a rich history in quantum chemistry. However, it seems to be the case that such methods are only suitable in particular situations and are not naturally equipped to simultaneously handle the variety of electron correlations that might be present in chemical systems. To this end, the current work seeks a computationally efficient, low-rank, "ultimate" coupled cluster method based exclusively on $T_2$ and its products which can effectively emulate more "complete" methods that explicitly consider higher-rank, $T_{2m}$ operators. We introduce a hierarchy of methods designed to systematically account for higher, even order cluster operators - like $T_4, T_6, \cdots, T_{2m}$ - by invoking tenets of the factorization theorem of perturbation theory and expectation-value coupled cluster theory. It is shown that each member within this methodological hierarchy is defined such that both the wavefunction and energy are correct through some order in many-body perturbation theory (MBPT), and can be extended up to arbitrarily high orders in $T_2$. The efficacy of such approximations are determined by studying the potential energy surface of several prototypical systems that are chosen to represent both non-dynamic, static, and dynamic correlation regimes. We find that the proposed hierarchy of augmented $T_2$ methods essentially reduce to standard CCD for problems where dynamic electron correlations dominate, but offer improvements in situations where non-dynamic and static correlations become relevant. A notable highlight of this work is that the cheapest methods in this hierarchy - which are correct through fifth-order in MBPT - consistently emulate the behavior of the $\mathcal{O}(N^{10})$ CCDQ method, yet only require a $\mathcal{O}(N^{6})$ algorithm by virtue of factorized intermediates.
△ Less
Submitted 11 July, 2024;
originally announced July 2024.
-
Factorized Quadruples and a Predictor of Higher-Level Correlation in Thermochemistry
Authors:
James H. Thorpe,
Zachary W. Windom,
Rodney J. Bartlett,
Devin A. Matthews
Abstract:
Coupled cluster theory has had a momentous impact on the ab initio prediction of molecular properties, and remains a staple ingratiate in high-accuracy thermochemical model chemistries. However, these methods require inclusion of at least some connected quadruple excitations, which generally scale at best as $\mathcal{O}(N^9)$ with the number of basis functions. It very difficult to predict, a pri…
▽ More
Coupled cluster theory has had a momentous impact on the ab initio prediction of molecular properties, and remains a staple ingratiate in high-accuracy thermochemical model chemistries. However, these methods require inclusion of at least some connected quadruple excitations, which generally scale at best as $\mathcal{O}(N^9)$ with the number of basis functions. It very difficult to predict, a priori, the effect correlation past CCSD(T) has on a give reaction energies. The purpose of this work is to examine cost-effective quadruple corrections based on the factorization theorem of many-body perturbation theory that may address these challenges. We show that the $\mathcal{O}(N^7)$, factorized CCSD(TQ${}_\text{f}$) method introduces minimal error to predicted correlation and reaction energies as compared to the $\mathcal{O}(N^9)$ CCSD(TQ). Further, we examine the performance of Goodson's continued fraction method in the estimation of CCSDT(Q)${}_Λ$ contributions to reaction energies, as well as a "new" method related to %TAE[(T)] that we refer to as a scaled perturbation estimator. We find that the scaled perturbation estimator based upon CCSD(TQ${}_\text{f}$)/cc-pVDZ is capable of predicting CCSDT(Q)${}_Λ$/cc-pVDZ contributions to reaction energies with an average error of 0.07 kcal mol${}^{-1}$ and a RMST of 0.52 kcal mol${}^{-1}$ when applied to a test-suite of nearly 3000 reactions. This offers a means by which to reliably ballpark how important post-CCSD(T) contributions are to reaction energies while incurring no more than CCSD(T) formal cost and a little mental math.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
An attractive way to correct for missing singles excitations in unitary coupled cluster doubles theory
Authors:
Zachary W. Windom,
Daniel Claudino,
Rodney J. Bartlett
Abstract:
Coupled cluster methods based exclusively on double excitations are comparatively "cheap" and interesting model chemistries, as they are typically able to capture the bulk of the dynamical electron correlation effects. The trade-off in such approximations is that the effect of neglected excitations, particularly single excitations, can be considerable. Using standard and electron pair-restricted…
▽ More
Coupled cluster methods based exclusively on double excitations are comparatively "cheap" and interesting model chemistries, as they are typically able to capture the bulk of the dynamical electron correlation effects. The trade-off in such approximations is that the effect of neglected excitations, particularly single excitations, can be considerable. Using standard and electron pair-restricted $T_2$ operators to define two flavors of unitary coupled cluster doubles (UCCD) methods, we investigate the extent in which missing single excitations can be recovered from low-order corrections in many-body perturbation theory (MBPT) within the unitary coupled cluster (UCC) formalism. Our analysis includes the derivations of finite-order, UCC energy functionals which are used as a basis to define perturbative estimates of missed single excitations. This leads to the novel UCCD[4S] and UCCD[6S] methods, which consider energy corrections for missing singles excitations through fourth- and sixth-order in MBPT, respectively. We also apply the same methodology to the electron pair-restricted ansatz, but the improvements are only marginal. Our findings show that augmenting UCCD with these post hoc perturbative corrections can lead to UCCSD-quality results.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
-
An assessment of frozen natural orbitals and band gaps using equation of motion coupled cluster theory: a case study on polyacene and trans-polyacetylene
Authors:
Zachary W. Windom,
AV Lam,
Ajith Perera,
Rodney J. Bartlett
Abstract:
Frozen natural orbitals (FNOs) are used to augment IP/EA-EOM-CCSD calculations targeting the band gap of trans-polyacetylene and polyacene. We show the resulting electron affinities (EAs), ionization potentials (IPs), and extrapolated band gaps incur errors that are largely tunable to a desired accuracy, yet require many orders of magnitude fewer core-hours as compared to the corresponding full ca…
▽ More
Frozen natural orbitals (FNOs) are used to augment IP/EA-EOM-CCSD calculations targeting the band gap of trans-polyacetylene and polyacene. We show the resulting electron affinities (EAs), ionization potentials (IPs), and extrapolated band gaps incur errors that are largely tunable to a desired accuracy, yet require many orders of magnitude fewer core-hours as compared to the corresponding full calculation. The relationship between various FNO truncation schemes and (cc-pV$n$Z) basis set is also examined.
△ Less
Submitted 13 February, 2024;
originally announced February 2024.
-
A new "gold standard": perturbative triples corrections in unitary coupled cluster theory and prospects for quantum computing
Authors:
Zachary W. Windom,
Daniel Claudino,
Rodney J. Bartlett
Abstract:
A major difficulty in quantum simulation is the adequate treatment of a large collection of entangled particles, synonymous with electron correlation in electronic structure theory, with coupled cluster (CC) theory being the leading framework in dealing with this problem. Augmenting computationally affordable low-rank approximations in CC theory with a perturbative account of higher-rank excitatio…
▽ More
A major difficulty in quantum simulation is the adequate treatment of a large collection of entangled particles, synonymous with electron correlation in electronic structure theory, with coupled cluster (CC) theory being the leading framework in dealing with this problem. Augmenting computationally affordable low-rank approximations in CC theory with a perturbative account of higher-rank excitations is a tractable and effective way of accounting for the missing electron correlation in those approximations. This is perhaps best exemplified by the "gold standard" CCSD(T) method, which bolsters the baseline CCSD with effects of triple excitations using considerations from many-body perturbation theory (MBPT). Despite this established success, such a synergy between MBPT and the unitary analog of CC theory (UCC) has not been explored. In this work, we propose a similar approach wherein converged UCCSD amplitudes, which can be obtained on a quantum computer, are leveraged by a classical computer to evaluate energy corrections associated with triple excitations - leading to the UCCSD[T] and UCCSD(T*) methods. The rationale behind these choices is shown to be rigorous by studying the properties of finite-order UCC energy functionals. Although our efforts do not support the addition of the fifth-order contribution as in the (T) correction, comparisons are nevertheless made using a hybrid UCCSD(T) approach. We assess the performance of these approaches on a collection of small molecules, and demonstrate the benefits of harnessing the inherent synergy between MBPT and UCC theories.
△ Less
Submitted 11 January, 2024;
originally announced January 2024.
-
Nonconvergence of the Feynman-Dyson diagrammatic perturbation expansion of propagators
Authors:
So Hirata,
Ireneusz Grabowski,
J. V. Ortiz,
Rodney J. Bartlett
Abstract:
Using a general-order ab initio many-body Green's function method, we numerically illustrate several pathological behaviors of the Feynman-Dyson diagrammatic perturbation expansion of one-particle many-body Green's functions as electron Feynman propagators. (i) The perturbation expansion of the frequency-dependent self-energy is not convergent at the exact self-energy in many frequency domains. (i…
▽ More
Using a general-order ab initio many-body Green's function method, we numerically illustrate several pathological behaviors of the Feynman-Dyson diagrammatic perturbation expansion of one-particle many-body Green's functions as electron Feynman propagators. (i) The perturbation expansion of the frequency-dependent self-energy is not convergent at the exact self-energy in many frequency domains. (ii) An odd-perturbation-order self-energy has a qualitatively wrong shape and, as a result, many roots of the corresponding Dyson equation are nonphysical in that the poles may be complex or residues can exceed unity or be negative. (iii) A higher even-order self-energy consists of vertical lines at many frequencies, predicting numerous phantom poles with zero residues. (iv) Infinite partial resummations of diagrams by vertex or edge renormalization tend to exacerbate these pathologies. (v) The nonconvergence is caused by the nonanalyticity of the rational-function form of the exact Green's function at many frequencies, where the radius of convergence of its Taylor expansion is zero. This is consistent with the fact that (vi) Padé approximants (power-series expansions of a rational function) can largely restore the correct shape and poles of the Green's function. Nevertheless, not only does the nonconvergence render higher-order Feynman-Dyson diagrammatic perturbation theory useless for many lower-lying ionization or higher-lying electron-attachment states, but it also calls into question the validity of its combined use with the ansätze requiring the knowledge of all poles and residues. Such ansätze include the Galitskii-Migdal identity, the self-consistent Green's function methods, and some models of the algebraic diagrammatic construction.
△ Less
Submitted 25 April, 2024; v1 submitted 5 December, 2023;
originally announced December 2023.
-
Point Absorber Limits to Future Gravitational-Wave Detectors
Authors:
W. Jia,
H. Yamamoto,
K. Kuns,
A. Effler,
M. Evans,
P. Fritschel,
R. Abbott,
C. Adams,
R. X. Adhikari,
A. Ananyeva,
S. Appert,
K. Arai,
J. S. Areeda,
Y. Asali,
S. M. Aston,
C. Austin,
A. M. Baer,
M. Ball,
S. W. Ballmer,
S. Banagiri,
D. Barker,
L. Barsotti,
J. Bartlett,
B. K. Berger,
J. Betzwieser
, et al. (176 additional authors not shown)
Abstract:
High-quality optical resonant cavities require low optical loss, typically on the scale of parts per million. However, unintended micron-scale contaminants on the resonator mirrors that absorb the light circulating in the cavity can deform the surface thermoelastically, and thus increase losses by scattering light out of the resonant mode. The point absorber effect is a limiting factor in some hig…
▽ More
High-quality optical resonant cavities require low optical loss, typically on the scale of parts per million. However, unintended micron-scale contaminants on the resonator mirrors that absorb the light circulating in the cavity can deform the surface thermoelastically, and thus increase losses by scattering light out of the resonant mode. The point absorber effect is a limiting factor in some high-power cavity experiments, for example, the Advanced LIGO gravitational wave detector. In this Letter, we present a general approach to the point absorber effect from first principles and simulate its contribution to the increased scattering. The achievable circulating power in current and future gravitational-wave detectors is calculated statistically given different point absorber configurations. Our formulation is further confirmed experimentally in comparison with the scattered power in the arm cavity of Advanced LIGO measured by in-situ photodiodes. The understanding presented here provides an important tool in the global effort to design future gravitational wave detectors that support high optical power, and thus reduce quantum noise.
△ Less
Submitted 17 September, 2021;
originally announced September 2021.
-
LIGOs Quantum Response to Squeezed States
Authors:
L. McCuller,
S. E. Dwyer,
A. C. Green,
Haocun Yu,
L. Barsotti,
C. D. Blair,
D. D. Brown,
A. Effler,
M. Evans,
A. Fernandez-Galiana,
P. Fritschel,
V. V. Frolov,
N. Kijbunchoo,
G. L. Mansell,
F. Matichard,
N. Mavalvala,
D. E. McClelland,
T. McRae,
A. Mullavey,
D. Sigg,
B. J. J. Slagmolen,
M. Tse,
T. Vo,
R. L. Ward,
C. Whittle
, et al. (172 additional authors not shown)
Abstract:
Gravitational Wave interferometers achieve their profound sensitivity by combining a Michelson interferometer with optical cavities, suspended masses, and now, squeezed quantum states of light. These states modify the measurement process of the LIGO, VIRGO and GEO600 interferometers to reduce the quantum noise that masks astrophysical signals; thus, improvements to squeezing are essential to furth…
▽ More
Gravitational Wave interferometers achieve their profound sensitivity by combining a Michelson interferometer with optical cavities, suspended masses, and now, squeezed quantum states of light. These states modify the measurement process of the LIGO, VIRGO and GEO600 interferometers to reduce the quantum noise that masks astrophysical signals; thus, improvements to squeezing are essential to further expand our gravitational view of the universe. Further reducing quantum noise will require both lowering decoherence from losses as well more sophisticated manipulations to counter the quantum back-action from radiation pressure. Both tasks require fully understanding the physical interactions between squeezed light and the many components of km-scale interferometers. To this end, data from both LIGO observatories in observing run three are expressed using frequency-dependent metrics to analyze each detector's quantum response to squeezed states. The response metrics are derived and used to concisely describe physical mechanisms behind squeezing's simultaneous interaction with transverse-mode selective optical cavities and the quantum radiation pressure noise of suspended mirrors. These metrics and related analysis are broadly applicable for cavity-enhanced optomechanics experiments that incorporate external squeezing, and -- for the first time -- give physical descriptions of every feature so far observed in the quantum noise of the LIGO detectors.
△ Less
Submitted 25 May, 2021;
originally announced May 2021.
-
Point absorbers in Advanced LIGO
Authors:
Aidan F. Brooks,
Gabriele Vajente,
Hiro Yamamoto,
Rich Abbott,
Carl Adams,
Rana X. Adhikari,
Alena Ananyeva,
Stephen Appert,
Koji Arai,
Joseph S. Areeda,
Yasmeen Asali,
Stuart M. Aston,
Corey Austin,
Anne M. Baer,
Matthew Ball,
Stefan W. Ballmer,
Sharan Banagiri,
David Barker,
Lisa Barsotti,
Jeffrey Bartlett,
Beverly K. Berger,
Joseph Betzwieser,
Dripta Bhattacharjee,
Garilynn Billingsley,
Sebastien Biscans
, et al. (176 additional authors not shown)
Abstract:
Small, highly absorbing points are randomly present on the surfaces of the main interferometer optics in Advanced LIGO. The resulting nano-meter scale thermo-elastic deformations and substrate lenses from these micron-scale absorbers significantly reduces the sensitivity of the interferometer directly though a reduction in the power-recycling gain and indirect interactions with the feedback contro…
▽ More
Small, highly absorbing points are randomly present on the surfaces of the main interferometer optics in Advanced LIGO. The resulting nano-meter scale thermo-elastic deformations and substrate lenses from these micron-scale absorbers significantly reduces the sensitivity of the interferometer directly though a reduction in the power-recycling gain and indirect interactions with the feedback control system. We review the expected surface deformation from point absorbers and provide a pedagogical description of the impact on power build-up in second generation gravitational wave detectors (dual-recycled Fabry-Perot Michelson interferometers). This analysis predicts that the power-dependent reduction in interferometer performance will significantly degrade maximum stored power by up to 50% and hence, limit GW sensitivity, but suggests system wide corrections that can be implemented in current and future GW detectors. This is particularly pressing given that future GW detectors call for an order of magnitude more stored power than currently used in Advanced LIGO in Observing Run 3. We briefly review strategies to mitigate the effects of point absorbers in current and future GW wave detectors to maximize the success of these enterprises.
△ Less
Submitted 25 March, 2021; v1 submitted 14 January, 2021;
originally announced January 2021.
-
A route to improving RPA excitation energies through its connection to equation-of-motion coupled cluster theory
Authors:
Varun Rishi,
Ajith Perera,
Rodney J. Bartlett
Abstract:
We revisit the connection between equation-of-motion coupled cluster (EOM-CC) and random phase approximation (RPA) explored recently by Berkelbach [J. Chem. Phys. 149, 041103 (2018)] and unify various methodological aspects of these diverse treatment of ground and excited states. The identity of RPA and EOM-CC based on the ring coupled cluster doubles is established with numerical results which wa…
▽ More
We revisit the connection between equation-of-motion coupled cluster (EOM-CC) and random phase approximation (RPA) explored recently by Berkelbach [J. Chem. Phys. 149, 041103 (2018)] and unify various methodological aspects of these diverse treatment of ground and excited states. The identity of RPA and EOM-CC based on the ring coupled cluster doubles is established with numerical results which was proved previously on theoretical grounds. We then introduce new approximations in EOM-CC and RPA family of methods, assess their numerical performance and explore a way to reap the benefits of such a connection to improve on excitation energies. Our results suggest that addition of perturbative corrections to account for double excitations and missing exchange effects could result in significantly improved estimates.
△ Less
Submitted 6 November, 2020; v1 submitted 1 August, 2020;
originally announced August 2020.
-
Improving the Robustness of the Advanced LIGO Detectors to Earthquakes
Authors:
Eyal Schwartz,
A Pele,
J Warner,
B Lantz,
J Betzwieser,
K L Dooley,
S Biscans,
M Coughlin,
N Mukund,
R Abbott,
C Adams,
R X Adhikari,
A Ananyeva,
S Appert,
K Arai,
J S Areeda,
Y Asali,
S M Aston,
C Austin,
A M Baer,
M Ball,
S W Ballmer,
S Banagiri,
D Barker,
L Barsotti
, et al. (174 additional authors not shown)
Abstract:
Teleseismic, or distant, earthquakes regularly disrupt the operation of ground--based gravitational wave detectors such as Advanced LIGO. Here, we present \emph{EQ mode}, a new global control scheme, consisting of an automated sequence of optimized control filters that reduces and coordinates the motion of the seismic isolation platforms during earthquakes. This, in turn, suppresses the differenti…
▽ More
Teleseismic, or distant, earthquakes regularly disrupt the operation of ground--based gravitational wave detectors such as Advanced LIGO. Here, we present \emph{EQ mode}, a new global control scheme, consisting of an automated sequence of optimized control filters that reduces and coordinates the motion of the seismic isolation platforms during earthquakes. This, in turn, suppresses the differential motion of the interferometer arms with respect to one another, resulting in a reduction of DARM signal at frequencies below 100\,mHz. Our method greatly improved the interferometers' capability to remain operational during earthquakes, with ground velocities up to 3.9\,$μ\mbox{m/s}$ rms in the beam direction, setting a new record for both detectors. This sets a milestone in seismic controls of the Advanced LIGO detectors' ability to manage high ground motion induced by earthquakes, opening a path for further robust operation in other extreme environmental conditions.
△ Less
Submitted 24 July, 2020;
originally announced July 2020.
-
Rigid platform for applying large tunable strains to mechanically delicate samples
Authors:
Joonbum Park,
Jack M. Bartlett,
Hilary M. L. Noad,
Alexander Stern,
Mark E. Barber,
Markus König,
Suguru Hosoi,
Takasada Shibauchi,
Andrew P. Mackenzie,
Alexander Steppke,
Clifford W. Hicks
Abstract:
Response to uniaxial stress has become a major probe of electronic materials. Tuneable uniaxial stress may be applied using piezoelectric actuators, and so far two methods have been developed to couple samples to actuators. In one, actuators apply force along the length of a free, beam-like sample, allowing very large strains to be achieved. In the other, samples are affixed directly to piezoelect…
▽ More
Response to uniaxial stress has become a major probe of electronic materials. Tuneable uniaxial stress may be applied using piezoelectric actuators, and so far two methods have been developed to couple samples to actuators. In one, actuators apply force along the length of a free, beam-like sample, allowing very large strains to be achieved. In the other, samples are affixed directly to piezoelectric actuators, allowing study of mechanically delicate materials. Here, we describe an approach that merges the two: thin samples are affixed to a substrate, that is then pressurized uniaxially using piezoelectric actuators. Using this approach, we demonstrate application of large elastic strains to mechanically delicate samples: the van der Waals-bonded material FeSe, and a sample of CeAuSb$_2$ that was shaped with a focused ion beam.
△ Less
Submitted 25 March, 2020;
originally announced March 2020.
-
Robust Modeling of Acoustic Phonon Transmission in Nanomechanical Structures
Authors:
John Bartlett,
Karwan Rostem,
Edward Wollack
Abstract:
The transmission of acoustic phonons is an important element in the design and performance of nano-mechanical devices operating in the mesoscopic limit. Analytic expressions for the power transmission coefficient, T, exist only in the low-frequency (quasi-static) limit described by thin-plate elastic theory, and for well-defined elastic wave-guiding geometries. We compare two numerical techniques…
▽ More
The transmission of acoustic phonons is an important element in the design and performance of nano-mechanical devices operating in the mesoscopic limit. Analytic expressions for the power transmission coefficient, T, exist only in the low-frequency (quasi-static) limit described by thin-plate elastic theory, and for well-defined elastic wave-guiding geometries. We compare two numerical techniques based on finite-element computations to determine the frequency dependence of T for arbitrary phonon scattering structures. Both methods take into account acoustic mode conversion to acoustic and optical modes. In one case, phase and amplitude of complex-valued reflected waves are determined and related to transmission through a Fresnel equation, while in the other the magnitude of the transmitted mechanical power is directly calculated. The numerical robustness of these methods is demonstrated by considering the transmission across an abrupt junction in a rectangular elastic beam, a well-known problem of considerable importance in mesoscopic device physics. The simulations presented extend the standard results for acoustic phonon transmission at an abrupt junction, and are in good agreement with analytic predictions from thin-plate elastic theory in the long-wavelength limit. More generally, the numerical methods developed provide an effective tool for calculating acoustic mode energy loss in nano-mechanical resonators through mode conversion and heat transfer in arbitrary mesoscopic structures.
△ Less
Submitted 6 December, 2018;
originally announced December 2018.
-
Review: Far-Infrared Instrumentation and Technology Development for the Next Decade
Authors:
Duncan Farrah,
Kimberly Ennico Smith,
David Ardila,
Charles M. Bradford,
Michael Dipirro,
Carl Ferkinhoff,
Jason Glenn,
Paul Goldsmith,
David Leisawitz,
Thomas Nikola,
Naseem Rangwala,
Stephen A. Rinehart,
Johannes Staguhn,
Michael Zemcov,
Jonas Zmuidzinas,
James Bartlett,
Sean Carey,
William J. Fischer,
Julia Kamenetzky,
Jeyhan Kartaltepe,
Mark Lacy,
Dariusz C. Lis,
Lisa Locke,
Enrique Lopez-Rodriguez,
Meredith MacGregor
, et al. (11 additional authors not shown)
Abstract:
Far-infrared astronomy has advanced rapidly since its inception in the late 1950's, driven by a maturing technology base and an expanding community of researchers. This advancement has shown that observations at far-infrared wavelengths are important in nearly all areas of astrophysics, from the search for habitable planets and the origin of life, to the earliest stages of galaxy assembly in the f…
▽ More
Far-infrared astronomy has advanced rapidly since its inception in the late 1950's, driven by a maturing technology base and an expanding community of researchers. This advancement has shown that observations at far-infrared wavelengths are important in nearly all areas of astrophysics, from the search for habitable planets and the origin of life, to the earliest stages of galaxy assembly in the first few hundred million years of cosmic history. The combination of a still developing portfolio of technologies, particularly in the field of detectors, and a widening ensemble of platforms within which these technologies can be deployed, means that far-infrared astronomy holds the potential for paradigm-shifting advances over the next decade. In this review, we examine current and future far-infrared observing platforms, including ground-based, sub-orbital, and space-based facilities, and discuss the technology development pathways that will enable and enhance these platforms to best address the challenges facing far-infrared astronomy in the 21st century.
△ Less
Submitted 8 January, 2019; v1 submitted 7 September, 2017;
originally announced September 2017.
-
Computing solvated excited states using fragment-effective-field coupled-cluster perturbation theory with application to the electronic spectra of nucleobases in water
Authors:
Jason N. Byrd,
Robert W. Molt, Jr.,
Rodney J. Bartlett,
Beverly A. Sanders
Abstract:
Molecular fragment or embedding methods are powerful techniques for overcoming scalability limitations in electronic structure theory by dividing large molecular systems into individual units that are small enough to be treated using standard techniques. Our group has developed a linear-scaling coupled-cluster based fragment effective field method, molecular cluster perturbation theory (MCPT), tha…
▽ More
Molecular fragment or embedding methods are powerful techniques for overcoming scalability limitations in electronic structure theory by dividing large molecular systems into individual units that are small enough to be treated using standard techniques. Our group has developed a linear-scaling coupled-cluster based fragment effective field method, molecular cluster perturbation theory (MCPT), that provides a systematically improvable description of the electronic wavefunction and is scalable to thousands of atoms using the Aces4 massively parallel software package. We will present a systematic computational benchmarking study of small organic molecule solvation free energies obtained from MCPT using various perturbation corrections using experimental and large-scale ab initio reference values. We will also present recent work on the relative crystal lattice polymorph energies of several molecular crystals where the roles of many-body dispersion and intermol ecular electron exchange will be investigated. The molecular solvent and crystal systems presented here are carefully chosen so as to provide a robust test set for the fragment electronic structure community.
△ Less
Submitted 10 August, 2017;
originally announced August 2017.
-
Exploring Cosmic Origins with CORE: The Instrument
Authors:
P. de Bernardis,
P. A. R. Ade,
J. J. A. Baselmans,
E. S. Battistelli,
A. Benoit,
M. Bersanelli,
A. Bideaud,
M. Calvo,
F. J. Casas,
G. Castellano,
A. Catalano,
I. Charles,
I. Colantoni,
F. Columbro,
A. Coppolecchia,
M. Crook,
G. D'Alessandro,
M. De Petris,
J. Delabrouille,
S. Doyle,
C. Franceschet,
A. Gomez,
J. Goupy,
S. Hanany,
M. Hills
, et al. (104 additional authors not shown)
Abstract:
We describe a space-borne, multi-band, multi-beam polarimeter aiming at a precise and accurate measurement of the polarization of the Cosmic Microwave Background. The instrument is optimized to be compatible with the strict budget requirements of a medium-size space mission within the Cosmic Vision Programme of the European Space Agency. The instrument has no moving parts, and uses arrays of diffr…
▽ More
We describe a space-borne, multi-band, multi-beam polarimeter aiming at a precise and accurate measurement of the polarization of the Cosmic Microwave Background. The instrument is optimized to be compatible with the strict budget requirements of a medium-size space mission within the Cosmic Vision Programme of the European Space Agency. The instrument has no moving parts, and uses arrays of diffraction-limited Kinetic Inductance Detectors to cover the frequency range from 60 GHz to 600 GHz in 19 wide bands, in the focal plane of a 1.2 m aperture telescope cooled at 40 K, allowing for an accurate extraction of the CMB signal from polarized foreground emission. The projected CMB polarization survey sensitivity of this instrument, after foregrounds removal, is 1.7 μK$\cdot$arcmin. The design is robust enough to allow, if needed, a downscoped version of the instrument covering the 100 GHz to 600 GHz range with a 0.8 m aperture telescope cooled at 85 K, with a projected CMB polarization survey sensitivity of 3.2 μK$\cdot$arcmin.
△ Less
Submitted 22 May, 2017; v1 submitted 5 May, 2017;
originally announced May 2017.
-
Quantum correlation measurements in interferometric gravitational wave detectors
Authors:
D. V. Martynov,
V. V. Frolov,
S. Kandhasamy,
K. Izumi,
H. Miao,
N. Mavalvala,
E. D. Hall,
R. Lanza,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
C. Adams,
R. X. Adhikari,
S. B. Anderson,
A. Ananyeva,
S. Appert,
K. Arai,
S. M. Aston,
S. W. Ballmer,
D. Barker,
B. Barr,
L. Barsotti,
J. Bartlett,
I. Bartos,
J. C. Batch
, et al. (177 additional authors not shown)
Abstract:
Quantum fluctuations in the phase and amplitude quadratures of light set limitations on the sensitivity of modern optical instruments. The sensitivity of the interferometric gravitational wave detectors, such as the Advanced Laser Interferometer Gravitational wave Observatory (LIGO), is limited by quantum shot noise, quantum radiation pressure noise, and a set of classical noises. We show how the…
▽ More
Quantum fluctuations in the phase and amplitude quadratures of light set limitations on the sensitivity of modern optical instruments. The sensitivity of the interferometric gravitational wave detectors, such as the Advanced Laser Interferometer Gravitational wave Observatory (LIGO), is limited by quantum shot noise, quantum radiation pressure noise, and a set of classical noises. We show how the quantum properties of light can be used to distinguish these noises using correlation techniques. Particularly, in the first part of the paper we show estimations of the coating thermal noise and gas phase noise, hidden below the quantum shot noise in the Advanced LIGO sensitivity curve. We also make projections on the observatory sensitivity during the next science runs. In the second part of the paper we discuss the correlation technique that reveals the quantum radiation pressure noise from the background of classical noises and shot noise. We apply this technique to the Advanced LIGO data, collected during the first science run, and experimentally estimate the quantum correlations and quantum radiation pressure noise in the interferometer for the first time.
△ Less
Submitted 10 February, 2017;
originally announced February 2017.
-
Predictive coupled-cluster isomer orderings for some Si${}_n$C${}_m$ ($m, n\le 12$) clusters; A pragmatic comparison between DFT and complete basis limit coupled-cluster benchmarks
Authors:
Jason N. Byrd,
Jesse J. Lutz,
Yifan Jin,
Duminda S. Ranasinghe,
John A. Montgomery Jr.,
Ajith Perera,
Xiaofeng F. Duan,
Larry W. Burggraf,
Beverly A. Sanders,
Rodney J. Bartlett
Abstract:
The accurate determination of the preferred ${\rm Si}_{12}{\rm C}_{12}$ isomer is important to guide experimental efforts directed towards synthesizing SiC nano-wires and related polymer structures which are anticipated to be highly efficient exciton materials for opto-electronic devices. In order to definitively identify preferred isomeric structures for silicon carbon nano-clusters, highly accur…
▽ More
The accurate determination of the preferred ${\rm Si}_{12}{\rm C}_{12}$ isomer is important to guide experimental efforts directed towards synthesizing SiC nano-wires and related polymer structures which are anticipated to be highly efficient exciton materials for opto-electronic devices. In order to definitively identify preferred isomeric structures for silicon carbon nano-clusters, highly accurate geometries, energies and harmonic zero point energies have been computed using coupled-cluster theory with systematic extrapolation to the complete basis limit for set of silicon carbon clusters ranging in size from SiC$_3$ to ${\rm Si}_{12}{\rm C}_{12}$. It is found that post-MBPT(2) correlation energy plays a significant role in obtaining converged relative isomer energies, suggesting that predictions using low rung density functional methods will not have adequate accuracy. Utilizing the best composite coupled-cluster energy that is still computationally feasible, entailing a 3-4 SCF and CCSD extrapolation with triple-$ζ$ (T) correlation, the {\it closo} ${\rm Si}_{12}{\rm C}_{12}$ isomer is identified to be the preferred isomer in support of previous calculations [J. Chem. Phys. 2015, 142, 034303]. Additionally we have investigated more pragmatic approaches to obtaining accurate silicon carbide isomer energies, including the use of frozen natural orbital coupled-cluster theory and several rungs of standard and double-hybrid density functional theory. Frozen natural orbitals as a way to compute post MBPT(2) correlation energy is found to be an excellent balance between efficiency and accuracy.
△ Less
Submitted 21 June, 2016;
originally announced June 2016.
-
The Sensitivity of the Advanced LIGO Detectors at the Beginning of Gravitational Wave Astronomy
Authors:
D. V. Martynov,
E. D. Hall,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
C. Adams,
R. X. Adhikari,
R. A. Anderson,
S. B. Anderson,
K. Arai,
M. A. Arain,
S. M. Aston,
L. Austin,
S. W. Ballmer,
M. Barbet,
D. Barker,
B. Barr,
L. Barsotti,
J. Bartlett,
M. A. Barton,
I. Bartos,
J. C. Batch,
A. S. Bell,
I. Belopolski,
J. Bergman
, et al. (239 additional authors not shown)
Abstract:
The Laser Interferometer Gravitational Wave Observatory (LIGO) consists of two widely separated 4 km laser interferometers designed to detect gravitational waves from distant astrophysical sources in the frequency range from 10 Hz to 10 kHz. The first observation run of the Advanced LIGO detectors started in September 2015 and ended in January 2016. A strain sensitivity of better than…
▽ More
The Laser Interferometer Gravitational Wave Observatory (LIGO) consists of two widely separated 4 km laser interferometers designed to detect gravitational waves from distant astrophysical sources in the frequency range from 10 Hz to 10 kHz. The first observation run of the Advanced LIGO detectors started in September 2015 and ended in January 2016. A strain sensitivity of better than $10^{-23}/\sqrt{\text{Hz}}$ was achieved around 100 Hz. Understanding both the fundamental and the technical noise sources was critical for increasing the observable volume in the universe. The average distance at which coalescing binary black hole systems with individual masses of 30 $M_\odot$ could be detected was 1.3 Gpc. Similarly, the range for binary neutron star inspirals was about 75 Mpc. With respect to the initial detectors, the observable volume of Universe increased respectively by a factor 69 and 43. These improvements allowed Advanced LIGO to detect the gravitational wave signal from the binary black hole coalescence, known as GW150914.
△ Less
Submitted 10 February, 2018; v1 submitted 1 April, 2016;
originally announced April 2016.
-
Calibration of the Advanced LIGO detectors for the discovery of the binary black-hole merger GW150914
Authors:
The LIGO Scientific Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
K. Ackley,
C. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
N. Aggarwal,
O. D. Aguiar,
A. Ain,
P. Ajith,
B. Allen,
P. A. Altin,
D. V. Amariutei,
S. B. Anderson,
W. G. Anderson,
K. Arai,
M. C. Araya,
C. C. Arceneaux,
J. S. Areeda,
K. G. Arun
, et al. (702 additional authors not shown)
Abstract:
In Advanced LIGO, detection and astrophysical source parameter estimation of the binary black hole merger GW150914 requires a calibrated estimate of the gravitational-wave strain sensed by the detectors. Producing an estimate from each detector's differential arm length control loop readout signals requires applying time domain filters, which are designed from a frequency domain model of the detec…
▽ More
In Advanced LIGO, detection and astrophysical source parameter estimation of the binary black hole merger GW150914 requires a calibrated estimate of the gravitational-wave strain sensed by the detectors. Producing an estimate from each detector's differential arm length control loop readout signals requires applying time domain filters, which are designed from a frequency domain model of the detector's gravitational-wave response. The gravitational-wave response model is determined by the detector's opto-mechanical response and the properties of its feedback control system. The measurements used to validate the model and characterize its uncertainty are derived primarily from a dedicated photon radiation pressure actuator, with cross-checks provided by optical and radio frequency references. We describe how the gravitational-wave readout signal is calibrated into equivalent gravitational-wave-induced strain and how the statistical uncertainties and systematic errors are assessed. Detector data collected over 38 calendar days, from September 12 to October 20, 2015, contain the event GW150914 and approximately 16 of coincident data used to estimate the event false alarm probability. The calibration uncertainty is less than 10% in magnitude and 10 degrees in phase across the relevant frequency band 20 Hz to 1 kHz.
△ Less
Submitted 28 February, 2017; v1 submitted 11 February, 2016;
originally announced February 2016.
-
Approximating electronically excited states with equation-of-motion linear coupled-cluster theory
Authors:
Jason N. Byrd,
Varun Rishi,
Ajith Perera,
Rodney J. Bartlett
Abstract:
A new perturbative approach to canonical equation-of-motion coupled-cluster theory is presented using coupled-cluster perturbation theory. A second-order Møller-Plesset partitioning of the Hamiltonian is used to obtain the well known equation-of-motion many-body perturbation theory (EOM-MBPT(2)) equations and two new equation-of-motion methods based on the linear coupled-cluster doubles (EOM-LCCD)…
▽ More
A new perturbative approach to canonical equation-of-motion coupled-cluster theory is presented using coupled-cluster perturbation theory. A second-order Møller-Plesset partitioning of the Hamiltonian is used to obtain the well known equation-of-motion many-body perturbation theory (EOM-MBPT(2)) equations and two new equation-of-motion methods based on the linear coupled-cluster doubles (EOM-LCCD) and linear coupled-cluster singles and doubles (EOM-LCCSD) wavefunctions. This is achieved by performing a short-circuiting procedure on the MBPT(2) similarity transformed Hamiltonian. These new methods are benchmarked against very accurate theoretical and experimental spectra from 25 small organic molecules. It is found that the proposed methods have excellent agreement with canonical EOM-CCSD state for state orderings and relative excited state energies as well as acceptable quantitative agreement for absolute excitation energies compared with the best estimate theory and experimental spectra.
△ Less
Submitted 7 July, 2015;
originally announced July 2015.
-
Molecular Cluster Perturbation Theory. I. Formalism
Authors:
Jason N. Byrd,
Nakul Jindal,
Robert W. Molt, Jr.,
Rodney J. Bartlett,
Beverly A. Sanders,
Victor F. Lotrich
Abstract:
We present second-order molecular cluster perturbation theory (MCPT(2)), a linear scaling methodology to calculate arbitrarily large systems with explicit calculation of individual wavefunctions in a coupled-cluster framework. This new MCPT(2) framework uses coupled-cluster perturbation theory and an expansion in terms of molecular dimer interactions to obtain molecular wavefunctions that are infi…
▽ More
We present second-order molecular cluster perturbation theory (MCPT(2)), a linear scaling methodology to calculate arbitrarily large systems with explicit calculation of individual wavefunctions in a coupled-cluster framework. This new MCPT(2) framework uses coupled-cluster perturbation theory and an expansion in terms of molecular dimer interactions to obtain molecular wavefunctions that are infinite-order in both the electronic fluctuation operator and all possible dimer (and products of dimers) interactions. The MCPT(2) framework has been implemented in the new SIA/Aces4 parallel architecture, making use of the advanced dynamic memory control and fine grained parallelism to perform very large explicit molecular cluster calculations. To illustrate the power of this method, we have computed energy shifts, lattice site dipole moments, and harmonic vibrational frequencies via explicit calculation of the bulk system for the polar and non-polar polymorphs of solid hydrogen fluoride. The explicit lattice size (without using any periodic boundary conditions) was expanded up to 1,000 HF molecules, with 32,000 basis functions and 10,000 electrons. Our obtained HF lattice site dipole moments and harmonic vibrational frequencies agree well with the existing literature.
△ Less
Submitted 23 March, 2015; v1 submitted 6 September, 2014;
originally announced September 2014.
-
Correlation correction to configuration interaction singles from coupled cluster perturbation theory
Authors:
Jason N. Byrd,
Victor F. Lotrich,
Rodney J. Bartlett
Abstract:
A new state specific correlation correction to configuration interaction singles (CIS) excitation energies is preseted using coupled cluster perturbation theory (CCPT). General expressions for CIS-CCPT are derived and expanded explicitly to first order in the wavefunction and second order in the energy. By virtue of the nature of CCPT this method is {\it a priori} size extensive and incorporates i…
▽ More
A new state specific correlation correction to configuration interaction singles (CIS) excitation energies is preseted using coupled cluster perturbation theory (CCPT). General expressions for CIS-CCPT are derived and expanded explicitly to first order in the wavefunction and second order in the energy. By virtue of the nature of CCPT this method is {\it a priori} size extensive and incorporates infinite order effects into the wavefunction. This results in a balanced singles space excited state theory that at second order is an improvement over the ubiquitous CIS(D) method and comparable in quality to equation of motion coupled cluster (EOM-CC). A modest test set composed of the first four excited states from nine small organic molecules was used to quantify the accuracy and consistency of the CIS-CCPT2 excitation energies and density of states. We find that CIS-CCPT2 has a standard deviation error of 0.18 eV for excitation energies and 0.14 eV for density of states compared to EOM-CC, a factor of two better than CIS(D) with a significant reduction in the maximum deviation as well
△ Less
Submitted 23 April, 2014;
originally announced April 2014.