Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “couple cluster approximation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

Application of Coupled-Cluster Methods to the Prediction and Interpretation of the Spectra of Molecules of Interest in Atmospheric Chemistry

The quality of fundamental vibrational frequencies determined using the CCSD(T) method (singles and doubles coupled-cluster theory plus a perturbational estimate of the effects of connected triple excitations) is shown to be very good, usually predicting band centers to within +/-8/cm. This approach is applied to several molecules of interest in atmospheric chemistry, including HNO, NO2(+), H2CO, and HOCl. The HNO molecule displays a large and unusual anharmonicity in the H-N stretch. For the calculation of ultraviolet (UV) spectra, the linear response CCSD (LRCCSD) approach (which is equivalent to EOM-CCSD) has been shown to yield vertical excitation energies that are accurate to approximately 0.1 eV for singly excited electronic states. This method together with more approximate methods is used to examine the UV spectra of several molecules important in stratospheric chemistry, including HOCl, Cl2O, ClOOCl, ClOOH, and HOOH.

Lee, Timothy J.↗

The Singlet–Triplet Gap of Cyclobutadiene: The CIPSI-Driven CC( P ; Q ) Study

An accurate determination of singlet−triplet gaps in biradicals, including cyclobutadiene in the automerization barrier region where one has to balance the substantial nondynamical many-electron correlation effects characterizing the singlet ground state with the predominantly dynamical correlations of the lowest-energy triplet, remains a challenge for many quantum chemistry methods. High-level coupled-cluster (CC) approaches, such as the CC method with a full treatment of singly, doubly, and triply excited clusters (CCSDT), are often capable of providing reliable results, but routine application of such methods is hindered by their high computational costs. We have recently proposed a practical alternative to converging the CCSDT energetics at small fractions of the computational effort, even when electron correlations become stronger and connected triply excited clusters are larger and nonperturbative, by merging the CC(P;Q) moment expansions with the selected configuration interaction methodology abbreviated as CIPSI. We demonstrate that one can accurately approximate the highly accurate CCSDT potential surfaces characterizing the lowest singlet and triplet states of cyclobutadiene along the automerization coordinate and the gap between them using tiny fractions of triply excited cluster amplitudes identified with the help of relatively inexpensive CIPSI Hamiltonian diagonalizations.

Basis sets↗

Analysis of spin frustration in an Fe III 7 cluster using a combination of computational, experimental, and magnetostructural correlation methods

The synthesis, structure, and magnetic properties are reported for [Fe 7 O 3 (O 2 C t Bu) 9 (mda) 3 (H 2 O) 3 ] ( 1 ), where mdaH 2 is N -methyldiethanolamine. 1 was prepared from the reaction of [Fe 3 O(O 2 C t Bu) 6 (H 2 O) 3 ](NO 3 ) with mdaH 2 in a 1:~3 ratio in MeCN. The core of 1 consists of a central octahedral Fe III ion held within a non-planar Fe 6 loop by three μ 3 -O 2- and three μ 2 -RO - arms from the three mda 2- chelates. Variable-temperature dc and ac magnetic susceptibility studies revealed dominant antiferromagnetic coupling, leading to a ground state spin of S = 5 / 2 . The ground state was confirmed by a fit of magnetization data collected in the 0.1–7.0 T and 1.8–10.0 K ranges. The four Fe 2 pairwise exchange parameters ( J 1 - J 4 ) were estimated by independent methods: theoretical calculations using either broken symmetry energy differences (-46.3, -16.2, -3.9, and - 28.1 cm -1 , respectively) or Green’s function approximation methods (-41.4, -14.8, -13.2, and - 24.7 cm -1 ), and a magnetostructural correlation (MSC) previously developed for high nuclearity Fe III /O complexes (-39.5, -13.8, -6.7, and - 23.5 cm -1 ). Additionally, the J 1 - J 4 obtained from the MSC and theoretical methods were used with the program PHI to both simulate χ M T vs T as well as to serve as reasonable input values to fit the experimental data (-41.0, -11.4, -5.0, and - 27.3 cm -1 ). Analysis of the J ij led to identification of the spin frustration effects operative and the resultant spin vector alignments at each Fe III ion, thus allowing for the rationalization of the experimental ground state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Systematic improvement of trial states in phaseless auxiliary-field quantum Monte Carlo

We extend the use of coupled cluster (CC) trial states in the phaseless auxiliary-field quantum Monte Carlo (AFQMC) method beyond single and double excitations to include both triple and quadruple excitations. With this AFQMC/CC hierarchy, we are able to systematically benchmark the method's performance on molecular systems as the quality of the trial is improved. Our results show that the phaseless AFQMC energy improves systematically and is typically significantly more accurate than the energy of the underlying trial state. However, the relative improvement compared to the trial CC energy decreases as we ascend the CC hierarchy. As the CC wavefunction is usually further approximated when used as an AFQMC trial, we also explore the relationship between the components of the CC wavefunction and the resulting AFQMC/CC error. Our results suggest that improving the representation of the CC wave function in the AFQMC trial does not always lower the bias even when it increases the fidelity of the trial with the exact ground state.

Chemical Physics (physics.chem-ph)↗

Ab Initio Studies of Fluorine and Chlorine Oxide and Nitrogen Oxide Species of Interest in Stratospheric Chemistry

The quality of fundamental vibrational frequencies determined using the CCSD(T) method (singles and doubles coupled-cluster theory plus a perturbational estimate of the effects of connected triple excitations) is shown to be very good, usually predicting band centers to within plus or minus 8 per centimeter. This approach is applied to several molecules of interest in atmospheric chemistry, such as HNO, cis-FONO, cis-ClONO, and ClOOH. The HNO molecule displays a large and unusual anharmonicity in the H-N stretch. For the calculation of ultraviolet (UV) spectra, the linear response CCSD (LRCCSD) approach (which is equivalent to EOM-CCSD) has been shown to yield vertical excitation energies that are accurate to approximately equal to 0.1 eV for singly excited electronic states. This method together with more approximate methods is used to examine the UV spectra of several molecules important in stratospheric chemistry, including HOCl, Cl2O, ClONO2, HONO2, ClOOCl, ClOOH, and HOOH.

Lee, Timothy J.↗

The Physical Origin of Galaxy Morphologies and Scaling Laws

We propose a numerical study designed to interpret the origin and evolution of galaxy properties revealed by space- and ground-based imaging and spectroscopical surveys. Our aim is to unravel the physical processes responsible for the development of different galaxy morphologies and for the establishment of scaling laws such as the Tully-Fisher relation for spirals and the Fundamental Plane of ellipticals. In particular, we plan to address the following major topics: (1) The morphology and observability of protogalaxies, and in particular the relationship between primordial galaxies and the z approximately 3 'Ly-break' systems identified in the Hubble Deep Field and in ground-based searches; (2) The origin of the disk and spheroidal components in galaxies, the timing and mode of their assembly, the corresponding evolution in galaxy morphologies and its sensitivity to cosmological parameters; (3) The origin and redshift evolution of the scaling laws that link the mass, luminosity size, stellar content, and metal abundances of galaxies of different morphological types. This investigation will use state-of-the-art N-body/gasdynamical codes to provide a spatially resolved description of the galaxy formation process in hierarchically clustering universes. Coupled with population synthesis techniques. our models can be used to provide synthetic 'observations' that can be compared directly with observations of galaxies both nearby and at cosmologically significant distances. This study will thus provide insight into the nature of protogalaxies and into the formation process of galaxies like our own Milky Way. It will also help us to assess the cosmological significance of these observations within the context of hierarchical theories of galaxy formation and will supply a theoretical context within which current and future observations can be interpreted.

Steinmetz, Matthias↗

Demonstration of Imaging Flow Diagnostics Using Rayleigh Scattering in Langley 0.3-Meter Transonic Cryogenic Tunnel

The feasibility of using the Rayleigh scattering technique for molecular density imaging of the free-stream flow field in the Langley 0.3-Meter Transonic Cryogenic Tunnel has been experimentally demonstrated. The Rayleigh scattering was viewed with a near-backward geometry with a frequency-doubled output from a diode-pumped CW Nd:YAG laser and an intensified charge-coupled device camera. Measurements performed in the range of free-stream densities from 3 x 10(exp 25) to 24 x 10(exp 25) molecules/cu m indicate that the observed relative Rayleigh signal levels are approximately linear with flow field density. The absolute signal levels agree (within approx. 30 percent) with the expected signal levels computed based on the well-known quantities of flow field density, Rayleigh scattering cross section for N2, solid angle of collection, transmission of the optics, and the independently calibrated camera sensitivity. These results show that the flow field in this facility is primarily molecular (i.e., not contaminated by clusters) and that Rayleigh scattering is a viable technique for quantitative nonintrusive diagnostics in this facility.

Shirinzadeh, B.↗

A review on recent machine learning applications for imaging mass spectrometry studies

Imaging mass spectrometry (IMS) is a powerful analytical technique widely used in biology, chemistry, and materials science fields that continue to expand. IMS provides a qualitative compositional analysis and spatial mapping with high chemical specificity. The spatial mapping information can be 2D or 3D depending on the analysis technique employed. Due to the combination of complex mass spectra coupled with spatial information, large high-dimensional datasets (hyperspectral) are often produced. Therefore, the use of automated computational methods for an exploratory analysis is highly beneficial. The fast-paced development of artificial intelligence (AI) and machine learning (ML) tools has received significant attention in recent years. These tools, in principle, can enable the unification of data collection and analysis into a single pipeline to make sampling and analysis decisions on the go. There are various ML approaches that have been applied to IMS data over the last decade. Here, in this review, we discuss recent examples of the common unsupervised (principal component analysis, non-negative matrix factorization, k-means clustering, uniform manifold approximation and projection), supervised (random forest, logistic regression, XGboost, support vector machine), and other methods applied to various IMS datasets in the past five years. The information from this review will be useful for specialists from both IMS and ML fields since it summarizes current and representative studies of computational ML-based exploratory methods for IMS.

47 OTHER INSTRUMENTATION↗

Disruption of the Globular Cluster Pal 5

Orbit calculations suggest that the sparse globular cluster, Pal 5, will pass within 7 kpc of the Galactic center the next time it crosses the plane, where it might be destroyed by tidal stresses. We study this problem, treating Pal 5 as a self-consistent dynamical system orbiting through an external potential that represents the Galaxy. The first part of the problem is to find suitable analytic approximations to the Galactic potential. They must be valid in all regions the cluster is likely to explore. Observed velocity and positional data for Pal 5 are used as initial conditions to determine the orbit. Methods we used for a different problem some 12 years ago have been adapted to this problem. Three experiments have been run, with M/L= 1, 3, and 10, for the cluster model. The cluster blew up shortly after passing through the Galactic plane (about 130 Myrs after the beginning of the run) with M/L=1. At M/L = 3 and 10 the cluster survived, although it got quite a kick in the fundamental mode on passing through the plane. But the fundamental mode oscillation died out in a couple of oscillation cycles at M/L=10. Pal 5 will probably be destroyed on its next crossing of the Galactic plane if M/L=1, but it can survive (albeit with fairly heavy damage) if NI/L=3. We haven't tried to trap the mass limits more closely than that. Pal 5 comes through pretty well unscathed at M/L=10. An interesting follow-up experiment would be to back the cluster up along its orbit to look at its previous passage through the Galactic plane, to see what kind of object it might have been at earlier times.

Miller, R. H.↗

Efficient Implementation for Unitary Coupled Cluster State Preparation for Near-Term Quantum Computers

Unitary coupled cluster theory (UCC) is a common wave function ansatz for quantum simulation of molecular electronic structure using the variational quantum eigenvalue solver (VQE). Even for small molecules using a double-ζ basis, the number of variational parameters required to minimize the electronic energy (i.e., optimize the circuit) is large and beyond the reach of current quantum computers. For example, a circuit simulating C2 using the UCCSD ansatz and the cc-pVDZ basis set with frozen-core will require over 10,000 variational parameters and a Hilbert space of over 10^8 determinants. To make progress on simulating such molecular systems on near-term quantum computers, we explore how much of the optimization can be approximately prepared with classical simulation while reducing the number of optimization steps performed on a quantum device. Recently, Chen, Cheng, and Freericks [J. Chem. Theory Comput. 2021, 17, 841-847] presented an algorithm for the factorized form of the UCC ansatz that allows for efficient UCC optimizations on classical hardware. We flip the algorithm around and use it to prepare approximate quantum circuits for systems that require a large number of qubits to represent. We will present results from our implementation and discuss strategies for incorporating this implementation for algorithms involving near-term quantum computers.

J Wayne Mullinax↗

Triple Excitations in Green’s Function Coupled Cluster Solver for Studies of Strongly Correlated Systems in the Framework of Self-Energy Embedding Theory

Embedding theories became important approaches used for accurate calculations of both molecules and solids. In these theories, a small chosen subset of orbitals is treated with an accurate method, called an impurity solver, capable of describing higher correlation effects. Ideally, such a chosen fragment should contain multiple orbitals responsible for the chemical and physical behavior of the compound. Handling a large number of chosen orbitals presents a very significant challenge for the current generation of solvers used in the physics and chemistry community. Here, we develop a Green’s function coupled cluster singles doubles and triples (GFCCSDT) solver that can be used for a quantitative description in both molecules and solids. This solver allows us to treat orbital spaces that are inaccessible to other accurate solvers. At the same time, GFCCSDT maintains high accuracy of the resulting self-energy. Moreover, in conjunction with the GFCCSD solver, it allows us to test the systematic convergence of computational studies. Developing the CC family of solvers paves the road to fully systematic Green’s function embedding calculations in solids. In this paper, we focus on the investigation of GFCCSDT self-energies for a strongly correlated problem of SrMnO 3 solid. Subsequently, we apply this solver to solid MnO showing that an approximate variant of GFCCSDT is capable of yielding a high accuracy orbital resolved spectral function.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Global Disease Outbreaks Associated with the 2015-2016 El Niño Event

Interannual climate variability patterns associated with the El Niño-Southern Oscillation phenomenon result in climate and environmental anomaly conditions in specific regions worldwide that directly favor outbreaks and/or amplification of variety of diseases of public health concern including chikungunya, hantavirus, Rift Valley fever, cholera, plague, and Zika. We analyzed patterns of some disease outbreaks during the strong 2015–2016 El Niño event in relation to climate anomalies derived from satellite measurements. Disease outbreaks in multiple El Niño-connected regions worldwide (including Southeast Asia, Tanzania, western US, and Brazil) followed shifts in rainfall, temperature, and vegetation in which both drought and flooding occurred in excess (14–81% precipitation departures from normal). These shifts favored ecological conditions appropriate for pathogens and their vectors to emerge and propagate clusters of diseases activity in these regions. Our analysis indicates that intensity of disease activity in some ENSO-teleconnected regions were approximately 2.5–28% higher during years with El Niño events than those without. Plague in Colorado and New Mexico as well as cholera in Tanzania were significantly associated with above normal rainfall (p < 0.05); while dengue in Brazil and southeast Asia were significantly associated with above normal land surface temperature (p < 0.05). Routine and ongoing global satellite monitoring of key climate variable anomalies calibrated to specific regions could identify regions at risk for emergence and propagation of disease vectors. Such information can provide sufficient lead-time for outbreak prevention and potentially reduce the burden and spread of ecologically coupled diseases.

El Nino↗

A CCD-based search for very low mass members of the Pleiades cluster

We have obtained deep charge coupled device (CCD)V and I images of a number of fields near the center of the Pleiades open cluster. We have also obtained imaging data for Praesepe, a very similar cluster in terms of distance and richness but nearly 10 times older than the Pleiades. Because brown dwarfs are predicted to become much fainter and cooler between Pleiades and Praesepe ages, this provides a powerful differential technique for placing constraints on the brown dwarf population in open clusters. Combined with our previously reported observations, we now have data for about 0.4 sq deg in the Pleiades, corresponding roughly to 5% of the area of that cluster. We have searched the new CCD frames for additional Pleiades brown dwarf candidates. Two possible candidates are present, the faintest of which has V approximately equal to 22.5, (V-I)(sub K) approximately equal to 4.6. Because we do not have proper motion data and the colors of these objects are not redder than the reddest known field stars, it is possible that some or all of our candidates are somewhat higher mass field stars rather than Pleiades-age brown dwarfs. Even if all six of the proposed brown dwarf candidates in our 0.4 sq deg field are Pleiades members, the relatively small number found suggests that low mass stars or brown dwarfs do not contribute significantly to the total mass of the cluster.

Stauffer, John R.↗

Rapidly convergent quantum Monte Carlo using a Chebyshev projector

The multireference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multireference problem. These include a logarithmically scaling metric-tree-based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants with essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.

Zhao, Zijun↗

Optical studies of gamma-ray bursts. 2: The superbowl burst GRB 930131

We report Charge Couple Devices (CCD) observations of the central region of the Compton Telescope (COMPTEL) error box for the 31 Jan 1993 gamma-ray burst. One optical transient was found. Its distance from the IPN track and Energy Gamma Ray Experiment Telescope (EGRET) error box make the association of this object with GRB 930131 doubtful. We present new CCD observations of two objects identified by Schaefer et al. (1993) that are located within a few arcminutes of the interplanetary network (IPN) that were reported to be possibly variable. We also report observations on the ROSAT All-sky Survey Extreme Ultraviolet (EUV) source noted by those investigators and the RS CVn star HD 106255. The ROSAT All-Sky Survey EUV source is suggested to arise from a faint, but relatively rich cluster of galaxies. HD 106255 is located approximately 7 min from the IPN. Our CCD results show that this star was bluer than previously reported, but these observations were obtained significantly after the occurrence of GRB 930131. It is not possible to estimate the color of HD 106255 when the burst occurred.

Mcnamara, Bernard↗

The evolution of the intergalactic medium and the origin of the galaxy luminosity function

The coupling of the Press and Schechter prescription with the CDM scenario and the Hoyle-Rees-Ostriker cooling criterion leads to a galaxy formation scenario in which galaxies are overproduced by a large factor. Although star formation might be suppressed in the smaller halos, a large amount of energy per galactic mass is needed to account for the present number density of galaxies. The evolution of the intergalactic medium (IGM) provides a simple criterion to prevent galaxy formation without requiring feedback, since halos with small virial temperatures are not able to retain the infalling hot gas of the IGM. If the ionizing background has decreased since z is approximately 1 - 2, then this criterion explains the slope of the luminosity function at the faint end. In addition, this scenario predicts two populations of dwarf galaxies, well differentiated in age, gas content, stellar populations, and clustering properties, which can be identified with dE and dIm galaxies.

Valls-Gabaud, David↗

On the composition of neuKREEP: QMD contamination at Apollo 11?

The Group A basalts of Apollo 11 differ in many respects from other high-Ti basalts of the region. Chemically, they are the only high-K (greater than 2000 ppm K) variety of high-Ti basalt and are enriched in incompatible trace elements relative to other basalts from both the Apollo 11 and Apollo 17 sites. In addition, Group A basalts are the youngest of all high-Ti basalts, with an age of 3.56 +/- 0.02 Ga. The cluster of compositions is consistent with the Apollo 11 Group A basalts representing a single flow. Papanastassiou et al. have also indicated the uniqueness of these basalts, based particularly on relatively young Rb-Sr model ages (3.8 - 3.9 Ga). A model for the formation of the Group A basalts was presented by Jerde et al., wherein the Apollo 17 orange volcanic glass is the parent liquid. Fractionation of this composition, coupled with the assimilation of incompatible-element-rich material, results in compositions akin to those of the Apollo 11 Group A basalt population. Orange glass of similar major-element composition is present at the Apollo 11 site as well, although complete trace element analyses are not available. New modelling results using the Apollo 11 orange glass major elements are grossly similar to those obtained using the Apollo 17 orange glass, indicating approximately 30 percent fractionation.

Jerde, Eric A.↗

Do We Know How to Model Reionization?

I compare the power spectra of the radiation fields from two recent sets of fully-coupled simulations that model cosmic reionization: “Cosmic Reionization On Computers” (CROC) and “Thesan”. While both simulations have similar power spectra of the radiation sources, the power spectra of the photoionization rate are significantly different at the same values of cosmic time or the same values of the mean neutral hydrogen fraction. However, the power spectra of the photoionization rate can be matched at large scales for the two simulations when the matching snapshots are allowed to vary independently. I.e., on large scales, the clustering of the radiation field in two simulations evolves similarly, but the exact timing of this evolution is different in different simulations and is not parameterized by an easily interpretable physical quantity like the mean neutral fraction or the mean free path. On small scales, large differences are present and remain partially unexplained. Both CROC and Thesan use the Variable Eddington Tensor approximation for modeling radiative transfer, but adopt different closure relations (optically thin OTVET versus M1). The role of this key difference is tested by using smaller simulations with a new cosmological simulation code that implements both closure relations in a controlled environment (the same hydro, cooling, and gravity solvers and the star formation recipe). In these controlled tests, both the M1 closure and the OTVET ansatz follow the expected behavior from a simple analytical approximation, demonstrating that the differences in the 2-point function of the radiation field induced by the choice of the Eddington tensor are not dominant.

79 ASTRONOMY AND ASTROPHYSICS↗