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.

168 records · Page 10

Convergent close-coupling calculations of electrons scattering on electronically excited molecular hydrogen

We use the adiabatic-nuclei molecular convergent close-coupling method to perform calculations of 0.01–1000 eV electrons scattering on the c 3 Π u , a 3 $Σ$$^{+}_{g}$, B 1 $Σ$$^{+}_{u}$, C 1 Π u , and EF 1 $Σ$$^{+}_{g}$ states of H 2 in the v=0 vibrational level. Elastic, superelastic, ionization, and grand-total cross sections are presented, as well as cross sections for excitation of the n=2–3 singlet and triplet states of H 2 (where n is the atomic-limit principle quantum number). Comparison with available theoretical results is made. Good agreement is found with the recent R-matrix results [J. Phys. B 53, 245203 (2020)] for most of the exchange and dipole-forbidden transitions, but not for the dipole-allowed transitions. Here, the sources of disagreement were found to be an unconverged partial-wave expansion and the utilization of the fixed-nuclei approximation (as opposed to adiabatic-nuclei) in the R-matrix calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗

NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts

In this work, we present NENCI-2021, a benchmark database of ~8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation–π and anion–π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion–π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (E int ) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The E int values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) E int value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent—these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)- based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation–π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the E int values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive E int values, and maximum errors of ±0.2–0.3 kcal/mol (i.e., ~±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Charge radii in covariant density functional theory: A global view

A systematic global investigation of differential charge radii has been performed within the CDFT framework for the first time. Theoretical results obtained with conventional covariant energy density functionals and the separable pairing interaction of Tian et al. [Phys. Lett. B 676, 44 (2009)] are compared with experimental differential charge radii in the regions of the nuclear chart in which available experimental data crosses the neutron shell closures at N=28,50,82, and 126. The analysis of absolute differential radii of different isotopic chains and their relative properties indicate clearly that such properties are reasonably well described in model calculations in the cases when the mean-field approximation is justified. However, while the observed clusterization of differential charge radii of different isotopic chains is well described above the N=50 and N=126 shell closures, it is more difficult to reproduce it above the N=28 and N=82 shell closures because of possible deficiencies in the underlying single-particle structure. The impact of the latter has been evaluated for spherical shapes and it was shown that the relative energies of the single-particle states and the patterns of their occupation with increasing neutron number have an appreciable impact on the evolution of the δ$\langle$r2$\rangle$N,N' values. These factors also limit the predictive power of model calculations in the regions of high densities of the single-particle states of different origin. It is shown that the kinks in the charge radii at neutron shell closures are due to the underlying single-particle structure and due to weakening or collapse of pairing at these closures. The regions of the nuclear chart in which the correlations beyond mean field are expected to have an impact on charge radii are indicated; the analysis shows that the assignment of a calculated excited prolate minimum to the experimental ground state allows us to understand the trends of the evolution of differential charge radii with neutron number in many cases of shape coexistence even at the mean-field level. It is usually assumed that pairing is a dominant contributor to odd-even staggering (OES) in charge radii. Our analysis paints a more complicated picture. It suggests a new mechanism in which the fragmentation of the single-particle content of the ground state in odd-mass nuclei due to particle-vibration coupling provides a significant contribution to OES in charge radii.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Self-Consistent Simulation of the Brownian Stage of Dust Growth

It is a widely accepted view that in proto-planetary accretion disks the collision and following sticking of dust particles embedded in the gas eventually leads to the formation of planetesimals (coagulation). For the smallest dust grains, Brownian motion is assumed to be the dominant source of their relative velocities leading to collisions between these dust grains. As the dust grains grow they eventually couple to the turbulent motion of the gas which then drives the coagulation much more efficiently. Many numerical coagulation simulations have been carried out to calculate the fractal dimension of the aggregates, which determines the duration of the ineffective Brownian stage of growth. Predominantly on-lattice and off-lattice methods were used. However, both methods require simplification of the astrophysical conditions. The aggregates found by those methods had a fractal dimension of approximately 2 which is equivalent to a constant, mass-independent friction time. If this value were valid for the conditions in an accretion disk, this would mean that the coagulation process would finally 'freeze out' and the growth of a planetesimal would be impossible within the lifetime of an accretion disk. In order to investigate whether this fractal dimension is model independent, we simulate self-consistently the Brownian stage of the coagulation by an N-particle code. This method has the advantage that no further assumptions about homogeneity of the dust have to be made. In our model, the dust grains are considered as aggregates built up of spheres. The equation of motion of the dust grains is based on the probability density for the diffusive transport within the gas atmosphere. Because of the very low number density of the dust grains, only 2-body-collisions have to be considered. As the Brownian stage of growth is very inefficient, the system is to be simulated over long periods of time. In order to find close particle pairs of the system which are most likely to undergo a collision, we use a particle-in-cell (PIC) method for the early stages of the simulation where the system is still very homogeneous and a tree method later when the particles are more clustered.

Kempf, S.↗

The RSA survey of dwarf galaxies, 1: Optical photometry

We present detailed surface photometry, based on broad B-band charge coupled device (CCD) images, of about 80 dwarf galaxies. Our sample represents approximately 10% of all dwarf galaxies identified in the vicinity of Revised Shapley-Ames (RSA) galaxies on high resolution blue photographic plates, referred to as the RSA survey of dwarf galaxies. We derive global properties and radial surface brightness profiles, and examine the morphologies. The radial surface brightness profiles of dwarf galaxies, whether early or late type, display the same varieties in shape and complexity as those of classical giant galaxies. Only a few are well described by a pure r(exp 1/4) law. Exponential profiles prevail. Features typical of giant disk galaxies, such as exponential profiles with a central depression, lenses, and even, in one case (IC 2041), a relatively prominent bulge are also found in dwarf galaxies. Our data suggest that the central region evolves from being bulge-like, with an r(exp 1/4) law profile, in bright galaxies to a lens-like structure in dwarf galaxies. We prove detailed surface photometry to be a helpful if not always sufficient tool in investigating the structure of dwarf galaxies. In many cases kinematic information is needed to complete the picture. We find the shapes of the surface brightness profiles to be loosely associated with morphological type. Our sample contains several new galaxies with properties intermediate between those of giant and dwarf ellipticals (but no M32-like objects). This shows that such intermediate galaxies exist so that at least a fraction of early-type dwarf ellipticals is structurally related to early-type giants instead of belonging to a totally unrelated, disjunct family. This supports an origin of early-type dwarf galaxies as originally more massive systems that acquired their current morphology as a result of substantial, presumable supernova-driven, mass loss. On the other hand, several early-type dwarfs in our sample are merger candidates. Merger events may lead to anisotropic velocity distributions in systems of any luminosity, including dwarfs. The RSA sample of dwarf galaxies is more likely to contain mergers because, in contrast to earlier dwarf galaxy surveys that have focused on clusters and rich groups of galaxies, the RSA dwarfs are typically located in low density environments. The occurrence of mergers among dwarf galaxies is of interest in connection with the rapid evolution of faint blue galaxy counts at redshift z less than 1 which suggests that dwarf galaxies were about five times more numerous in the recent past. Finally, our sample contains several examples of late-type dwarfs and 'transition' types that are potential precursors of nucleated early-type dwarfs. All the above processes--mass loss, mergers, astration--are likely to have contributed to the formation of the current population of diffuse early-type dwarfs. A few new redshifts of dwarf galaxies are reported in this paper.

Vader, J. Patricia↗

Perspective on Kramers symmetry breaking and restoration in relativistic electronic structure methods for open-shell systems

Without rigorous symmetry constraints, solutions to approximate electronic structure methods may artificially break symmetry. In the case of the relativistic electronic structure, if time-reversal symmetry is not enforced in calculations of molecules not subject to a magnetic field, it is possible to artificially break Kramers degeneracy in open shell systems. This leads to a description of excited states that may be qualitatively incorrect. Despite this, different electronic structure methods to incorporate correlation and excited states can partially restore Kramers degeneracy from a broken symmetry solution. For single-reference techniques, the inclusion of double and possibly triple excitations in the ground state provides much of the needed correction. Formally, however, this imbalanced treatment of the Kramers-paired spaces is a multi-reference problem, and so methods such as complete-active-space methods perform much better at recovering much of the correct symmetry by state averaging. Using multi-reference configuration interaction, any additional corrections can be obtained as the solution approaches the full configuration interaction limit. A recently proposed “Kramers contamination” value is also used to assess the magnitude of symmetry breaking.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗