Search NASA⌕ Search

SEARCH · Search NASA

Results for “OCCUPATION NUMBER”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

High-Temperature Water Adsorption Isotherms and Ambient Temperature Water Diffusion Rates on Water Harvesting Metal–Organic Frameworks

Water adsorption isotherms from 25 to 125 °C were measured for three metal−organic frameworks (MOFs), MOF-303, MOF-LA2-1, and MIL-100(Fe), which are frequently studied for water harvesting applications. The results show how the step in the water adsorption isotherm varies as a function of temperature and detail the combination of pressure and temperature necessary to remove adsorbed water. Furthermore, isobaric−isothermal Gibbs ensemble Monte Carlo simulations performed for MOF-303 shed light on the change in occupation numbers of the different known water adsorption sites with increasing temperature. Additionally, the diffusion rates of water through these materials were measured using concentration swing frequency response, and micropore diffusion was identified as the controlling mechanism. The Darken relation was used to show the dependence of the diffusion rate on the concentration and the impact of the adsorption isotherm slope. The adsorption of water on MOF-LA2-1 is faster than that on MIL-100(Fe). These data show that MOF-LA2-1 with its high-water adsorption capacity, quick adsorption rate, and favorable desorption energetics is a leading candidate for atmospheric water harvesting.

Adsorption↗
DOE OSTI2019

Distinguishing artificial and essential symmetry breaking in a single determinant: approach and application to the C 60 , C 36 , and C 20 fullerenes

We present a thorough analysis of symmetry breaking observed in Hartree–Fock (HF) solutions of the C 60 , C 36 , and C 20 fullerenes in order to characterize the nature of electron correlation in them. Our analysis is based on (1) the critical regularization strength to restore symmetry breaking in the recently developed regularized orbital optimized second-order Møller–Plesset perturbation theory (κ-OOMP2), (2) singlet–triplet gaps from various MP2 methods, and (3) natural orbital occupation numbers from restricted coupled-cluster with singles and doubles (RCCSD) and coupled-cluster valence bond with singles and doubles (CCVB-SD). Based on these three independent probes, we conclude that C 36 (D6h) exhibits genuine strong correlation and symmetry breaking whereas C 60 exhibits artificial HF symmetry breaking and is not strongly correlated. Investigating the critical regularization strength, we discuss strong correlation in C 20 at the Jahn–Teller distorted geometries (C 2h , D 2h , C i , and D 3h ) and the Ih geometry. Only C 20 (I h ) was found to be strongly correlated while others exhibit artificial HF symmetry breaking. This analysis highlights a useful feature of the recommended κ-OOMP2 method. It is an electronic structure method that describes dynamic correlation, and attenuates strong correlation in MP2 towards zero by regularization. Therefore, κ-OOMP2 will exhibit symmetry breaking in its reference determinant only when correlation is strong (i.e., essential symmetry breaking). Furthermore, artificial symmetry breaking (arising in HF due to neglect of dynamic correlation) thus appears to be removed in κ-OOMP2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Developing new and understanding old approximations in TDDFT

When a system has evolved far from a ground-state, the adiabatic approximations commonly used in time-dependent density functional theory calculations completely fail in some applications, while giving qualitatively good predictions in others, and sometimes even quantitative predictions. It is not clearly understood why this is so, and developing practical approximations going beyond the adiabatic approximation remains a challenge. This paper explores three different lines of investigation. First, an expression for the exact time-dependent exchange–correlation potential suggests that the accuracy of an adiabatic approximation is intimately related to the deviation between the natural orbital occupation numbers of the physical system and those of the Kohn–Sham system, and we explore this on some exactly-solvable model systems. The exact expression further suggests a path to go beyond the adiabatic approximations, and in the second part we discuss a newly proposed class of memory-dependent approximations developed in this way. Finally, we derive a new expression for the exact exchange–correlation potential from a coupling-constant path integration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

The intrapair electron correlation in natural orbital functional theory

A previously proposed [M. Piris, X. Lopez, F. Ruipérez, J. M. Matxain, and J. M. Ugalde, J. Chem. Phys. 134, 164102 (2011)] formulation of the two-particle cumulant, based on an orbital-pairing scheme, is extended here for including more than two natural orbitals. This new approximation is used to reconstruct the two-particle reduced density matrix (2-RDM) constrained to the D, Q, and G positivity necessary conditions of the N-representable 2-RDM. In this way, we have derived an extended version of the Piris natural orbital functional 5 (PNOF5e). An antisymmetrized product of strongly orthogonal geminals with the expansion coefficients explicitly expressed by the occupation numbers is also used to generate the PNOF5e. The theory is applied to the homolytic dissociation of selected diatomic molecules: H{sub 2}, LiH, and Li{sub 2}. The Bader's theory of atoms in molecules is used to analyze the electron density and the presence of non-nuclear maxima in the case of a set of light atomic clusters: Li{sub 2}, Li {sub 3}{sup +}, Li {sub 4}{sup 2+}, and H{sub 3}{sup +}. The improvement of PNOF5e over PNOF5 was observed by visualizing the electron densities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Excitation energies with linear response density matrix functional theory along the dissociation coordinate of an electron-pair bond in N-electron systems

Time dependent density matrix functional theory in its adiabatic linear response formulation delivers exact excitation energies ω{sub α} and oscillator strengths f{sub α} for two-electron systems if extended to the so-called phase including natural orbital (PINO) theory. The Löwdin-Shull expression for the energy of two-electron systems in terms of the natural orbitals and their phases affords in this case an exact phase-including natural orbital functional (PILS), which is non-primitive (contains other than just J and K integrals). In this paper, the extension of the PILS functional to N-electron systems is investigated. With the example of an elementary primitive NO functional (BBC1) it is shown that current density matrix functional theory ground state functionals, which were designed to produce decent approximations to the total energy, fail to deliver a qualitatively correct structure of the (inverse) response function, due to essential deficiencies in the reconstruction of the two-body reduced density matrix (2RDM). We now deduce essential features of an N-electron functional from a wavefunction Ansatz: The extension of the two-electron Löwdin-Shull wavefunction to the N-electron case informs about the phase information. In this paper, applications of this extended Löwdin-Shull (ELS) functional are considered for the simplest case, ELS(1): one (dissociating) two-electron bond in the field of occupied (including core) orbitals. ELS(1) produces high quality ω{sub α}(R) curves along the bond dissociation coordinate R for the molecules LiH, Li{sub 2}, and BH with the two outer valence electrons correlated. All of these results indicate that response properties are much more sensitive to deficiencies in the reconstruction of the 2RDM than the ground state energy, since derivatives of the functional with respect to both the NOs and the occupation numbers need to be accurate.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2013

A non-adiabatic dynamics study of octatetraene: The radiationless conversion from S{sub 2} to S{sub 1}

Simulation of the excited state dynamics of all-trans-1,3,5,7-octatetraene has been performed to investigate the ultrafast radiationless S{sub 2} → S{sub 1} internal conversion process. Multireference configuration interaction with single excitation method has been employed to optimize the equilibrium structure of the excited states, as well as the S{sub 2}/S{sub 1} conical intersection, and to investigate the non-adiabatic molecular dynamics of the S{sub 2}/S{sub 1} state transition. At the conical intersection, the molecule is found to be distorted from the original planar trans structure to a nearly perpendicular conformation around C{sub 3}−C{sub 4} bond, with the torsion angle being about 107°. Such structural change can result in mutual approaching of states S{sub 2} and S{sub 1} in energy, and drastically increase the nonadiabatic coupling between the two states by destroying the inter-state symmetry prohibition in the electronic wavefunctions. Surface-hopping molecular dynamics simulations are performed to describe the non-adiabatic process. Upon the Franck-Condon excitation to the S{sub 2} state, the molecule quickly twists its C{sub 3}−C{sub 4} bond and approaches the conical intersection region, where it can undergo efficient internal conversion to S{sub 1}. The decay time constant (τ) of S{sub 2} state is estimated to be around 251 fs by fitting the occupation number of average fraction of trajectories using an exponential damping function. This value is reasonably consistent with previous experimental measurements of around 300–400 fs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

The effect of classical and quantum dynamics on vibrational frequency shifts of H{sub 2} in clathrate hydrates

Vibrational frequency shifts of H{sub 2} in clathrate hydrates are important to understand the properties and elucidate details of the clathrate structure. Experimental spectra of H{sub 2} in clathrate hydrates have been measured for different clathrate compositions, temperatures, and pressures. In order to establish reliable relationships between the clathrate structure, dynamics, and observed frequencies, calculations of vibrational frequency shifts in different clathrate environments are required. In this study, a combination of classical molecular dynamics simulations, electronic structure calculations, and quantum dynamical simulation is used to calculate relative vibrational frequencies of H{sub 2} in clathrate hydrates. This approach allows us to assess dynamical effects and simulate the change of vibrational frequencies with temperature and pressure. The frequency distributions of the H{sub 2} vibrations in the different clathrate cage types agree favorably with experiment. Also, the simulations demonstrate that H{sub 2} in the 5{sup 12} cage is more sensitive to the details of the environment and to quantum dynamical effects, in particular when the cage is doubly occupied. We show that for the 5{sup 12} cage quantum effects lead to frequency increases and double occupation is unlikely. This is different for the 5{sup 12}6{sup 4} cages for which higher occupation numbers than one H{sub 2} per cage are likely.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Investigation of site preference of Zn doped Ba{sub 3}Co{sub 2−x}Zn{sub x}Fe{sub 24}O{sub 41} by Mössbauer spectroscopy

The polycrystalline Ba{sub 3}Co{sub 2−x}Zn{sub x}Fe{sub 24}O{sub 41} (x = 0.0, 0.5, 1.0) samples were prepared by using solid-state-reaction method. The crystal structures and magnetic properties of samples were investigated with x-ray diffractometer, vibrating sample magnetometer, and Mössbauer spectroscopy. The crystal structure of Ba{sub 3}Co{sub 2−x}Zn{sub x}Fe{sub 24}O{sub 41} (x = 0.0, 0.5, 1.0) samples was determined to be a hexagonal structure with P6{sub 3}/mmc space group at 295 K, and the saturation magnetization (M{sub s}) of Ba{sub 3}Co{sub 2−x}Zn{sub x}Fe{sub 24}O{sub 41} (x = 0.0, 0.5, 1.0) samples were found to be M{sub s} = 50.9, 53.1, 55.0 emu/g, respectively. From the temperature dependence of magnetization curves under 100 Oe between 4.2 and 740 K, we were able to observe the spin transition, and both spin transition temperature (T{sub s}) and Curie temperature (T{sub C}) decrease with increasing Zn concentration. Mössbauer spectra of all samples were obtained and analyzed at various temperatures ranging from 4.2 to 295 K. With ten-sextets for Fe sites corresponding to the Z-type hexagonal crystallographic sites, all spectra below T{sub C} were fitted by least-square method. In addition, from the site occupation numbers of Fe, calculated from the relative areas fitted to the Mössbauer spectra, we find that Zn ions preferentially occupy the tetrahedral sublattices of down sites.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Kinetic and interaction components of the exact time-dependent correlation potential

The exact exchange-correlation (xc) potential of time-dependent density functional theory has been shown to have striking features. For example, step and peak features are generically found when the system is far from its ground-state, and these depend nonlocally on the density in space and time. We analyze the xc potential by decomposing it into kinetic and interaction components and comparing each with their exact-adiabatic counterparts, for a range of dynamical situations in model one-dimensional two-electron systems. We find that often, but not always, the kinetic contribution is largely responsible for these features that are missed by the adiabatic approximation. The adiabatic approximation often makes a smaller error for the interaction component, which we write in two parts, one being the Coulomb potential due to the time-dependent xc hole. Non-adiabatic features of the kinetic component were also larger than those of the interaction component in cases that we studied when there is negligible step structure. In ground-state situations, step and peak structures arise in cases of static correlation, when more than one determinant is essential to describe the interacting state. We investigate the time-dependent natural orbital occupation numbers and find the corresponding relation between these and the dynamical step is more complex than for the ground-state case.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

A quasi-classical mapping approach to vibrationally coupled electron transport in molecular junctions

We develop a classical mapping approach suitable to describe vibrationally coupled charge transport in molecular junctions based on the Cartesian mapping for many-electron systems [B. Li and W. H. Miller, J. Chem. Phys. 137, 154107 (2012)]. To properly describe vibrational quantum effects in the transport characteristics, we introduce a simple transformation rewriting the Hamiltonian in terms of occupation numbers and use a binning function to facilitate quantization. The approach provides accurate results for the nonequilibrium Holstein model for a range of bias voltages, vibrational frequencies, and temperatures. It also captures the hallmarks of vibrational quantum effects apparent in step-like structure in the current-voltage characteristics at low temperatures as well as the phenomenon of Franck-Condon blockade.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Accurate prediction of diradical chemistry from a single-reference density-matrix method: Model application to the bicyclobutane to gauche-1,3-butadiene isomerization

Multireference correlation in diradical molecules can be captured by a single-reference 2-electron reduced-density-matrix (2-RDM) calculation with only single and double excitations in the 2-RDM parametrization. The 2-RDM parametrization is determined by N-representability conditions that are non-perturbative in their treatment of the electron correlation. Conventional single-reference wave function methods cannot describe the entanglement within diradical molecules without employing triple- and potentially even higher-order excitations of the mean-field determinant. In the isomerization of bicyclobutane to gauche-1,3-butadiene the parametric 2-RDM (p2-RDM) method predicts that the diradical disrotatory transition state is 58.9 kcal/mol above bicyclobutane. This barrier is in agreement with previous multireference calculations as well as recent Monte Carlo and higher-order coupled cluster calculations. The p2-RDM method predicts the Nth natural-orbital occupation number of the transition state to be 0.635, revealing its diradical character. The optimized geometry from the p2-RDM method differs in important details from the complete-active-space self-consistent-field geometry used in many previous studies including the Monte Carlo calculation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Benchmark calculations on the lowest-energy singlet, triplet, and quintet states of the four-electron harmonium atom

For a wide range of confinement strengths ω, explicitly-correlated calculations afford approximate energies E(ω) of the ground and low-lying excited states of the four-electron harmonium atom that are within few μhartree of the exact values, the errors in the respective energy components being only slightly higher. This level of accuracy constitutes an improvement of several orders of magnitude over the previously published data, establishing a set of benchmarks for stringent calibration and testing of approximate electronic structure methods. Its usefulness is further enhanced by the construction of differentiable approximants that allow for accurate computation of E(ω) and its components for arbitrary values of ω. The diversity of the electronic states in question, which involve both single- and multideterminantal first-order wavefunctions, and the availability of the relevant natural spinorbitals and their occupation numbers make the present results particularly useful in research on approximate density-matrix functionals. The four-electron harmonium atom is found to possess the {sup 3}P{sub +} triplet ground state at strong confinements and the {sup 5}S{sub −} quintet ground state at the weak ones, the energy crossing occurring at ω ≈ 0.0240919.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Approximate spin projected spin-unrestricted density functional theory method: Application to diradical character dependences of second hyperpolarizabilities

We develop a novel method for the calculation and the analysis of the one-electron reduced densities in open-shell molecular systems using the natural orbitals and approximate spin projected occupation numbers obtained from broken symmetry (BS), i.e., spin-unrestricted (U), density functional theory (DFT) calculations. The performance of this approximate spin projection (ASP) scheme is examined for the diradical character dependence of the second hyperpolarizability (γ) using several exchange-correlation functionals, i.e., hybrid and long-range corrected UDFT schemes. It is found that the ASP-LC-UBLYP method with a range separating parameter μ = 0.47 reproduces semi-quantitatively the strongly-correlated [UCCSD(T)] result for p-quinodimethane, i.e., the γ variation as a function of the diradical character.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Reduced density-matrix functional theory: Correlation and spectroscopy

In this work, we explore the performance of approximations to electron correlation in reduced density-matrix functional theory (RDMFT) and of approximations to the observables calculated within this theory. Our analysis focuses on the calculation of total energies, occupation numbers, removal/addition energies, and spectral functions. We use the exactly solvable Hubbard dimer at 1/4 and 1/2 fillings as test systems. This allows us to analyze the underlying physics and to elucidate the origin of the observed trends. For comparison, we also report the results of the GW approximation, where the self-energy functional is approximated, but no further hypothesis is made concerning the approximations of the observables. In particular, we focus on the atomic limit, where the two sites of the dimer are pulled apart and electrons localize on either site with equal probability, unless a small perturbation is present: this is the regime of strong electron correlation. In this limit, using the Hubbard dimer at 1/2 filling with or without a spin-symmetry-broken ground state allows us to explore how degeneracies and spin-symmetry breaking are treated in RDMFT. We find that, within the used approximations, neither in RDMFT nor in GW, the signature of strong correlation is present, when looking at the removal/addition energies and spectral function from the spin-singlet ground state, whereas both give the exact result for the spin-symmetry broken case. Moreover, we show how the spectroscopic properties change from one spin structure to the other.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Orbitals from local RDMFT: Are they Kohn-Sham or natural orbitals?

Recently, an approximate theoretical framework was introduced, called local reduced density matrix functional theory (local-RDMFT), where functionals of the one-body reduced density matrix (1-RDM) are minimized under the additional condition that the optimal orbitals satisfy a single electron Schrödinger equation with a local potential. In the present work, we focus on the character of these optimal orbitals. In particular, we compare orbitals obtained by local-RDMFT with those obtained with the full minimization (without the extra condition) by contrasting them against the exact NOs and orbitals from a density functional calculation using the local density approximation (LDA). We find that the orbitals from local-RMDFT are very close to LDA orbitals, contrary to those of the full minimization that resemble the exact NOs. Since local RDMFT preserves the good quality of the description of strong static correlation, this finding opens the way to a mixed density/density matrix scheme, where Kohn-Sham orbitals obtain fractional occupations from a minimization of the occupation numbers using 1-RDM functionals. This will allow for a description of strong correlation at a cost only minimally higher than a density functional calculation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Thermodynamic stability and structural properties of cluster crystals formed by amphiphilic dendrimers

We pursue the goal of finding real-world examples of macromolecular aggregates that form cluster crystals, which have been predicted on the basis of coarse-grained, ultrasoft pair potentials belonging to a particular mathematical class [B. M. Mladek et al., Phys. Rev. Lett. 46, 045701 (2006)]. For this purpose, we examine in detail the phase behavior and structural properties of model amphiphilic dendrimers of the second generation by means of monomer-resolved computer simulations. On augmenting the density of these systems, a fluid comprised of clusters that contain several overlapping and penetrating macromolecules is spontaneously formed. Upon further compression of the system, a transition to multi-occupancy crystals takes place, the thermodynamic stability of which is demonstrated by means of free-energy calculations, and where the FCC is preferred over the BCC-phase. Contrary to predictions for coarse-grained theoretical models in which the particles interact exclusively by effective pair potentials, the internal degrees of freedom of these molecules cause the lattice constant to be density-dependent. Furthermore, the mechanical stability of monodisperse BCC and FCC cluster crystals is restricted to a bounded region in the plane of cluster occupation number versus density. The structural properties of the dendrimers in the dense crystals, including their overall sizes and the distribution of monomers are also thoroughly analyzed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Self-consistent implementation of ensemble density functional theory method for multiple strongly correlated electron pairs

Here, the spin-restricted ensemble-referenced Kohn-Sham (REKS) method is based on an ensemble representation of the density and is capable of correctly describing the non-dynamic electron correlation stemming from (near-)degeneracy of several electronic configurations. The existing REKS methodology describes systems with two electrons in two fractionally occupied orbitals. In this work, the REKS methodology is extended to treat systems with four fractionally occupied orbitals accommodating four electrons and self-consistent implementation of the REKS(4,4) method with simultaneous optimization of the orbitals and their fractional occupation numbers is reported. The new method is applied to a number of molecular systems where simultaneous dissociation of several chemical bonds takes place, as well as to the singlet ground states of organic tetraradicals 2,4-didehydrometaxylylene and 1,4,6,9-spiro[4.4]nonatetrayl.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Next generation extended Lagrangian first principles molecular dynamics

Extended Lagrangian Born-Oppenheimer molecular dynamics [Phys. Rev. Lett., 100, 123004 (2008)] is formulated for general Hohenberg-Kohn density functional theory and compared to the extended Lagrangian framework of rst principles molecular dynamics by Car and Parrinello [Phys. Rev. Lett. 55, 2471 (1985)]. It is shown how extended Lagrangian Born-Oppenheimer molecular dynamics overcomes several shortcomings of regular, direct Born-Oppenheimer molecular dynamics, while improving or maintaining important features of Car-Parrinello simulations. The accuracy of the electronic degrees of freedom in extended Lagrangian Born-Oppenheimer molecular dynamics, with respect to the exact Born-Oppenheimer solution, is of second order in the size of the integration time step and of fourth order in the potential energy surface. Improved stability over recent formulations of extended Lagrangian Born-Oppenheimer molecular dynamics is achieved by generalizing the theory to finite temperature ensembles, using fractional occupation numbers in the calculation of the inner-product kernel of the extended harmonic oscillator that appears as a preconditioner in the electronic equations of motion. Materials systems that normally exhibit slow self-consistent field convergence can be simulated using integration time steps of the same order as in direct Born- Oppenheimer molecular dynamics, but without the requirement of an iterative, non-linear electronic ground state optimization prior to the force evaluations and without a systematic drift in the total energy. In combination with proposed low-rank and on-the- y updates of the kernel, this formulation provides an efficient and general framework for quantum based Born-Oppenheimer molecular dynamics simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗