SEARCH · Engineering Papers

Results for “Physics Regularization”

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.

SourceResearch connectionFrom yearThrough year

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

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Particulate Effluent Characterization (Final Report)

Understanding particulate emissions from nuclear facilities could help differentiate between normal nuclear operations and potential nuclear accidents or nuclear weapons tests. Oak Ridge National Laboratory (ORNL) is unique in that there are several types of nuclear facilities on site: an operating production reactor, radiochemical separation facilities, and a spallation neutron source. This project deployed a high-volume particulate air filter sampler to collect airborne particulate effluent from the nuclear facilities on site at ORNL. Collections occurred regularly from November 1, 2024, through June 30, 2025, and were analyzed via gamma spectroscopy in the laboratory. The radioisotope iodine-123 ( 123 I) was detected in several samples throughout the collection period. Detailed atmospheric transport modeling was performed on all detections for source attribution, and the most likely source of the 123 I was determined to be the Spallation Neutron Source. The project demonstrated the viability of ORNL as a test bed for effluent monitoring studies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Prediction of DIII-D Pedestal Structure from Externally Controllable Parameters

The sharp increase of pressure at the edge of a high confinement mode (H-mode) plasma, the pedestal, strongly impacts overall plasma performance. Predicting the pedestal is a necessity to control and optimize tokamak operations. An experimental data-driven machine learning (ML) approach is presented that predicts the pedestal heights and widths of electron density (ne) and electron temperature (Te) profiles as well as the separatrix ne from externally controllable parameters such as the plasma shape, heating method and power, and gas puff rate and integrated gas puff. The OMFIT framework was used with DIII-D data to efficiently, robustly, and automatically build a database of pedestal parameters to train machine learning models. Database creation was enabled by the search engine tool for DIII-D data, TokSearch, which parallelizes data fetching, enabling fast searches through basic signals of thousands of DIII-D shots and selection of relevant time intervals. Principal Component Analysis (PCA) separated the database into three clusters that represent classes of plasma shapes that are regularly used in DIII-D. The most important parameters for setting the pedestal structure were plasma current (Ip), toroidal magnetic field (Bφ), neutral beam heating power (PNBI) and shaping quantities. The Deep Jointly Informed Neural Networks (DJINN) algorithm was applied to identify suitable neural network (NN) architectures that appropriately capture the features of the pedestal database. Separate NNs were implemented for each pedestal parameter, and ensembling methods were used to improve the prediction accuracy and allowed estimation of the prediction uncertainty. The pedestal predictions of the test dataset lie within the measurement uncertainties of the pedestal parameters. The NN outperformed simple Linear Regression (LR) analysis, indicating non-linear dependencies in the pedestal structure. The presented achievements illustrate a promising path for future research, using feature extraction to infer experimental trends and thereby improve pedestal models as well as deploying NN for a fast pedestal prediction in DIII-D scenario development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Controlled Self-Assembly of Gold Nanotetrahedra into Quasicrystals and Complex Periodic Supracrystals

The self-assembly of shape-anisotropic nanocrystals into large-scale structures is a versatile and scalable approach to creating multifunctional materials. The tetrahedral geometry is ubiquitous in natural and manmade materials, yet regular tetrahedra present a formidable challenge in understanding their self-assembly behavior as they do not tile space. Here, we report diverse supracrystals from gold nanotetrahedra including the quasicrystal (QC) and the dimer packing predicted more than a decade ago and hitherto unknown phases. We solve the complex three-dimensional (3D) structure of the QC by a combination of electron microscopy, tomography, and synchrotron X-ray scattering. Nanotetrahedron vertex sharpness, surface ligands, and assembly conditions work in concert to regulate supracrystal structure. We also discover that the surface curvature of supracrystals can induce structural changes of the QC tiling and eventually, for small supracrystals with high curvature, stabilize a hexagonal approximant. Finally, our findings bridge the gap between computational design and experimental realization of soft matter assemblies and demonstrate the importance of accurate control over nanocrystal attributes and the assembly conditions to realize increasingly complex nanopolyhedron supracrystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Prediction of DIII-D Pedestal Structure From Externally Controllable Parameters

The sharp increase of pressure at the edge of a high confinement mode (H-mode) plasma, the pedestal, strongly impacts overall plasma performance. Predicting the pedestal is a necessity to control and optimize tokamak operations. Here, an experimental data-driven machine learning (ML) approach is presented that predicts the pedestal heights and widths of electron density (n e ) and electron temperature (T e ) profiles as well as the separatrix ne from externally controllable parameters such as the plasma shape, heating method and power, and gas puff rate and integrated gas puff. The OMFIT framework was used with DIII-D data to efficiently, robustly, and automatically build a database of pedestal parameters to train machine learning models. Database creation was enabled by the search engine tool for DIII-D data, TokSearch, which parallelizes data fetching, enabling fast searches through basic signals of thousands of DIII-D shots and selection of relevant time intervals. Principal Component Analysis (PCA) separated the database into three clusters that represent classes of plasma shapes that are regularly used in DIII-D. The most important parameters for setting the pedestal structure were plasma current (I p ), toroidal magnetic field (B Φ ), neutral beam heating power (P NBI ) and shaping quantities. The Deep Jointly Informed Neural Networks (DJINN) algorithm was applied to identify suitable neural network (NN) architectures that appropriately capture the features of the pedestal database. Separate NNs were implemented for each pedestal parameter, and ensembling methods were used to improve the prediction accuracy and allowed estimation of the prediction uncertainty. The pedestal predictions of the test dataset lie within the measurement uncertainties of the pedestal parameters. The NN outperformed simple Linear Regression (LR) analysis, indicating non-linear dependencies in the pedestal structure. The presented achievements illustrate a promising path for future research, using feature extraction to infer experimental trends and thereby improve pedestal models as well as deploying NN for a fast pedestal prediction in DIII-D scenario development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Dark Matter Reconstruction in LBAI Experiments with Imperfect Data

Long-baseline atom interferometer (LBAI) experiments offer unprecedented sensitivity to ultralight scalar dark matter (DM) [1], however reconstruction of a putative DM signal with traditional frequency-domain analysis requires ``perfect data (i.e., regularly-sampled with no missing samples). In a real LBAI experiment, there will undoubtedly be imperfections in the data leading to downtime. This downtime can arise from operational considerations (e.g., maintenance), the operational environment (motion of people and animals [2] or elevators), and robustness of the experimental apparatus (e.g., bad atom launches). In this work, we investigate the impact of various downtime models on the overall DM sensitivity of such an experiment. We compare the sensitivity for each downtime model as determined by a ``compound FFT analysis to a baseline no-downtime case. We also show how much sensitivity can be regained by moving to a Lomb-Scargle frequency analysis, as in [2]. Furthermore, we demonstrate reconstruction of the DM wave s phase as well as its frequency. [1] D. Antypas, et al, ``New Horizons: Scalar and Vector Ultralight Dark Matter (2022). arXiv:2203.14915 [2] J. Carlton and C. McCabe, ``From RATs to riches: mitigating anthropogenic and synanthropic noise in atom interferometer searches for ultra-light dark matter (2023). arXiv:2308.101731

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Relativistic Douglas–Kroll–Hess calculations of hyperfine interactions within first-principles multireference methods

A relativistic magnetic hyperfine interaction Hamiltonian based on the Douglas–Kroll–Hess (DKH) theory up to the second order is implemented within the ab initio multireference methods, including spin–orbit coupling in the Molcas/OpenMolcas package. This implementation is applied to calculate relativistic hyperfine coupling (HFC) parameters for atomic systems and diatomic radicals with valence s or d orbitals by systematically varying active space size in the restricted active space self-consistent field formalism with restricted active space state interaction for spin–orbit coupling. The DKH relativistic treatment of the hyperfine interaction reduces the Fermi contact contribution to the HFC due to the presence of kinetic factors that regularize the singularity of the Dirac delta function in the nonrelativistic Fermi contact operator. This effect is more prominent for heavier nuclei. As the active space size increases, the relativistic correction of the Fermi contact contribution converges well to the experimental data for light and moderately heavy nuclei. The relativistic correction, however, does not significantly affect the spin-dipole contribution to the hyperfine interaction. In addition to the atomic and molecular systems, the implementation is applied to calculate the relativistic HFC parameters for large trivalent and divalent Tb-based single-molecule magnets (SMMs), such as Tb(III)Pc2 and Tb(II)(CpiPr5)2 without ligand truncation using well-converged basis sets. In particular, for the divalent SMM, which has an unpaired valence 6s/5d hybrid orbital, the relativistic treatment of HFC is crucial for a proper description of the Fermi contact contribution. Even with the relativistic hyperfine Hamiltonian, the divalent SMM is shown to exhibit strong tunability of HFC via an external electric field (i.e., strong hyperfine Stark effect).

Chemistry↗

DOE OSTI2024

Investigation of Fe-Ni Battery/Module for Grid Service Duty Cycles

Iron–nickel (Fe-Ni) batteries are renowned for their durability and resilience against overcharging and operating temperatures. However, they encounter challenges in achieving widespread adoption for energy storage applications due to their low efficiency and the need for regular maintenance and electrolyte replacement, which adds to maintenance costs. This study evaluates and demonstrates the capabilities of Fe-Ni batteries for participating in grid energy storage applications. Stable performance was observed frequency regulation (FR) testing at 100% and 50% state of charge (SOC)s, while at 50% SOC, there was a 14% increase in efficiency compared to 100% SOC. Although 25% SOC achieved higher efficiency, limited cyclability was observed due to reaching the discharge cutoff voltage. Optimal SOC selection, battery monitoring, maintenance, and appropriate charging strategies of Fe-Ni batteries seem to be crucial for their FR applications. Fe-Ni batteries exhibit stable peak shaving (PS) results, indicating their suitability and reliability under various load conditions for PS testing. Extended cycling tests confirm their potential for long-term grid-scale energy storage, enhancing their appeal for PS and FR applications.

25 ENERGY STORAGE↗

DOE OSTI2022

Hybrid quantum-classical algorithms for approximate graph coloring

We show how to apply the recursive quantum approximate optimization algorithm (RQAOA) to MAX- k -CUT, the problem of finding an approximate k -vertex coloring of a graph. We compare this proposal to the best known classical and hybrid classical-quantum algorithms. First, we show that the standard (non-recursive) QAOA fails to solve this optimization problem for most regular bipartite graphs at any constant level p : the approximation ratio achieved by QAOA is hardly better than assigning colors to vertices at random. Second, we construct an efficient classical simulation algorithm which simulates level- 1 QAOA and level- 1 RQAOA for arbitrary graphs. In particular, these hybrid algorithms give rise to efficient classical algorithms, and no benefit arising from the use of quantum mechanics is to be expected. Nevertheless, they provide a suitable testbed for assessing the potential benefit of hybrid algorithm: We use the simulation algorithm to perform large-scale simulation of level- 1 QAOA and RQAOA with up to 300 qutrits applied to ensembles of randomly generated 3 -colorable constant-degree graphs. We find that level- 1 RQAOA is surprisingly competitive: for the ensembles considered, its approximation ratios are often higher than those achieved by the best known generic classical algorithm based on rounding an SDP relaxation. This suggests the intriguing possibility that higher-level RQAOA may be a potentially useful algorithm for NISQ devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Lattice realizations of topological defects in the critical (1+1)-d three-state Potts model

Topological/perfectly-transmissive defects play a fundamental role in the analysis of the symmetries of two dimensional conformal field theories (CFTs). In the present work, spin chain regularizations for these defects are proposed and analyzed in the case of the three-state Potts CFT. In particular, lattice versions for all the primitive defects are presented, with the remaining defects obtained from the fusion of the primitive ones. The defects are obtained by introducing modified interactions around two given sites of an otherwise homogeneous spin chain with periodic boundary condition. The various primitive defects are topological on the lattice except for one, which is topological only in the scaling limit. The lattice models are analyzed using a combination of exact diagonalization and density matrix renormalization group techniques. Low-lying energy spectra for different defect Hamiltonians as well as entanglement entropy of blocks located symmetrically around the defects are computed. The latter provides a convenient way to compute the g-function which characterizes various defects. Finally, the eigenvalues of the line operators in the “crossed channel” and fusion of different defect lines are also analyzed. The results are all in agreement with expectations from conformal field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Ligand Controls Excited Charge Carrier Dynamics in Metal-Rich CdSe Quantum Dots: Computational Insights

Small metal-rich semiconducting quantum dots (QDs) are promising for solid-state lighting and single-photon emission due to their highly tunable yet narrow emission line widths. Nonetheless, the anionic ligands commonly employed to passivate these QDs exert a substantial influence on the optoelectronic characteristics, primarily owing to strong electron–phonon interactions. In this work, we combine time-domain density functional theory and nonadiabatic molecular dynamics to investigate the excited charge carrier dynamics of Cd 28 Se 17 X 22 QDs (X = HCOO – , OH – , Cl – , and SH – ) at ambient conditions. These chemically distinct but regularly used molecular groups influence the dynamic surface-ligand interfacial interactions in Cd-rich QDs, drastically modifying their vibrational characteristics. The strong electron–phonon coupling leads to substantial transient variations at the band edge states. The strength of these interactions closely depends on the physicochemical characteristics of passivating ligands. Consequently, the ligands largely control the nonradiative recombination rates and emission characteristics in these QDs. Our simulations indicate that Cd 28 Se 17 (OH) 22 has the fastest nonradiative recombination rate due to the strongest electron–phonon interactions. Conversely, QDs passivated with thiolate or chloride exhibit considerably longer carrier lifetimes and suppressed nonradiative processes. The ligand-controlled electron–phonon interactions further give rise to the broadest and narrowest intrinsic optical line widths for OH and Cl-passivated single QDs, respectively. Finally, obtained computational insights lay the groundwork for designing appropriate passivating ligands on metal-rich QDs, making them suitable for a wide range of applications, from blue LEDs to quantum emitters.

36 MATERIALS SCIENCE↗

DOE OSTI2023

5ω Optical Thomson Scattering Report

The 5ω Optical Thomson Scattering system on the National Ignition Facility (NIF) is a diagnostic designed to measure temporally and spatially resolved plasma conditions in Inertial Confinement Fusion (ICF) Hohlraums. The system was proposed in 2014 and a phased approach to implementation was developed. In phase one the collection system was designed, built, and fielded on the NIF to be used for 3ω Thomson scattering measurements and background measurements near 211 nm (5ω. The initial commissioning experiment for the collection system was completed in Oct. 2016. Commissioning of the collection system continued through 2017 and the system is now fielded for a range of user experiments and regularly produces publication quality data. A dedicated 5ω probe laser was designed and built to allow Thomson scattering measurements in the presence of 2 MJ of 3ω drive energy. Due to scattered light (spectral reflections, laser-plasma instabilities, and unconverted drive energy) from the 3w drive lasers typical wavelengths (2ω and 4ω) fielded at other laser facilities like the Nova Laser Facility and the Omega Laser Facility were unable to meeting the signal to background requirements in design studies. A 10 Joule 5ω probe was proposed as a solution that met requirements. This 10 Joule laser was more energetic than previously fielded 5w lasers by 2-3 orders of magnitude. As part of a risk reduction plan, phase two of the project was to develop a 1 Joule, 5ω laser to demonstrate conversion efficiency >20% from 1ω to 5ω and make initial Thomson scattering measurements in the first few nanoseconds of an ICF laser pulse. Initial tests of the 5w conversion were completed at the Laboratory for Laser Energetics (LLE) and produced record 5ω energies. Based on these results from LLE, a 1 Joule, 5w laser system was designed for the NIF in 2018 and commissioning began in 2019. Commissioning of the 5ω laser system continued through 2023 and was eventually paused in Q1 of FY24 due to completing resource constraints. Commissioning of the 5w laser system proved incredibly challenging. Multiple issues were identified during commissioning and resolved, but issues remain. Currently there is not a clear understanding of why 5ω scattered has not been detected on the OTS collection system. Based on offline measurements, calculations, and preshot measurements the system appears to meet all requirements. Potential target physics issues have been investigated and do not appear to be an issue. The current hypothesis is that there is an error in one or more aspects of the offline testing not translating to expected performance when the system is fielded on the NIF. Additional full system testing in-situ utilizing the complete OTS system and NIF target chamber center time is needed to further test potential failure modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Dark matter or regular matter in neutron stars? How to tell the difference from the coalescence of compact objects

The mirror twin Higgs model is a candidate for (strongly-interacting) complex dark matter, which mirrors SM interactions with heavier quark masses. A consequence of this model are mirror neutron stars—exotic stars made entirely of mirror matter, which are significantly smaller than neutron stars and electromagnetically dark. This makes mergers of two mirror neutron stars detectable and distinguishable in gravitational wave observations, but can we observationally distinguish between regular neutron stars and those that may contain some mirror matter? This is the question we study in this paper, focusing on two possible realizations of mirror matter coupled to standard model matter within a compact object: (i) mirror matter captured by a neutron star and (ii) mirror neutron star-neutron star coalescences. Regarding (i), here we find that (nonrotating) mirror-matter-admixed neutron stars no longer have a single mass-radius sequence, but rather exist in a two-dimensional mass-radius plane. Regarding (ii), we find that binary systems with mirror neutron stars would span a much wider range of chirp masses and completely different binary Love relations, allowing merger remnants to be very light black holes. The implications of this are that gravitational wave observations with advanced LIGO and Virgo, and x-ray observations with NICER, could detect or constrain the existence of mirror matter through searches with wider model and parameter priors.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Covalent bond shortening and distortion induced by pressurization of thorium, uranium, and neptunium tetrakis aryloxides

Covalency involving the 5f orbitals is regularly invoked to explain the reactivity, structure and spectroscopic properties of the actinides, but the ionic versus covalent nature of metal-ligand bonding in actinide complexes remains controversial. The tetrakis 2,6-di- tert -butylphenoxide complexes of Th, U and Np form an isostructural series of crystal structures containing approximately tetrahedral MO 4 cores. We show that up to 3 GPa the Th and U crystal structures show negative linear compressibility as the OMO angles distort. At 3 GPa the angles snap back to their original values, reverting to a tetrahedral geometry with an abrupt shortening of the M-O distances by up to 0.1 Å. The Np complex shows similar but smaller effects, transforming above 2.4 GPa. Electronic structure calculations associate the M-O bond shortening with a change in covalency resulting from increased contributions to the M-O bonding by the metal 6d and 5f orbitals, the combination promoting MO 4 flexibility at little cost in energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Synthesis and Self-Assembly of Monodisperse Graphene Nanoribbons: Access to Submicron Architectures with Long-Range Order and Uniform Orientation

Fabricating organic semiconducting materials into large-scale well-organized architectures is critical for building high performance molecular electronics. While graphene nanoribbons (GNRs) hold enormous promise for various device applications, their assembly into a well-structured monolayer or multilayer architecture poses a substantial challenge. Here we report the preparation of length-defined monodisperse GNRs via the integrated iterative binomial synthesis (IIBS) strategy and their self-assembly into submicron-architectures with long-range order, uniform orientation as well as regular layers. Further, the use of short alkyl side chains benefits forming stable multi-layers through interlocking structures. By changing the length and backbone shapes of these monodisperse GNRs, various three-dimensional assemblies including multilayer stripes, monolayer stripes, and nanowires, can be achieved, leading to different photophysical properties and band gaps. The discovery of these intriguing self-assembly behaviors of length-defined GNRs is expected to open the door for various future applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Assessing the effect of regularization on the molecular properties predicted by SCAN and self-interaction corrected SCAN meta-GGA

Recent regularization of the SCAN meta-GGA functional (rSCAN) has simplified the numerical complexities of the SCAN functional, alleviating SCAN's stringent demand on the numerical integration grids to some extent. The regularization of rSCAN, however, results in the breaking of some constraints such as the uniform electron gas limit, the slowly varying density limit, and coordinate scaling of the iso-orbital indicator. Here, we assess the effects of regularization on the electronic, structural, vibrational, and magnetic properties of molecules by comparing the SCAN and rSCAN predictions. The properties studied include atomic energies, atomization energies, ionization potentials, electron affinities, barrier heights, infrared intensities, dissociation and reaction energies, spin moments of molecular magnets, and isomer ordering of water clusters. Our results show that rSCAN requires less dense numerical grids and gives very similar results to those of SCAN for all properties examined with the exception of atomization energies, which are worsened in rSCAN. We also examine the performance of self-interaction-corrected (SIC) rSCAN with respect to SIC-SCAN using the Perdew–Zunger (PZ) SIC method. The PZSIC method uses orbital densities to compute one-electron self-interaction errors and places an even more stringent demand on numerical grids. Our results show that SIC-rSCAN gives marginally better performance than SIC-SCAN for almost all properties studied in this work with numerical grids that are on average half or less as dense as that needed for SIC-SCAN.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

BPS Wilson loop in $\mathcal{N}$ = 2 superconformal SU(N) “orientifold” gauge theory and weak-strong coupling interpolation

We consider the expectation value $ \left\langle \mathcal{W}\right\rangle $ of the circular BPS Wilson loop in $ \mathcal{N} $ = 2 superconformal SU( N ) gauge theory containing a vector multiplet coupled to two hypermultiplets in rank-2 symmetric and antisymmetric representations. This theory admits a regular large N expansion, is planar-equivalent to $ \mathcal{N} $ = 4 SYM theory and is expected to be dual to a certain orbifold/orientifold projection of AdS 5 × S 5 superstring theory. On the string theory side $ \left\langle \mathcal{W}\right\rangle $ is represented by the path integral expanded near the same AdS 2 minimal surface as in the maximally supersymmetric case. Following the string theory argument in [5], we suggest that as in the $ \mathcal{N} $ = 4 SYM case and in the $ \mathcal{N} $ = 2 SU( N ) × SU( N ) superconformal quiver theory discussed in [19], the coefficient of the leading non-planar 1/ N 2 correction in $ \left\langle \mathcal{W}\right\rangle $ should have the universal λ 3/2 scaling at large ’t Hooft coupling. We confirm this prediction by starting with the localization matrix model representation for $ \left\langle \mathcal{W}\right\rangle $ . We complement the analytic derivation of the λ 3/2 scaling by a numerical high-precision resummation and extrapolation of the weak-coupling expansion using conformal mapping improved Padé analysis.

1/N Expansion↗

DOE OSTI2020

Utilizing Essential Symmetry Breaking in Auxiliary-Field Quantum Monte Carlo: Application to the Spin Gaps of the C 36 Fullerene and an Iron Porphyrin Model Complex

We present three distinct examples where phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) can be reliably performed with a single-determinant trial wave function with essential symmetry breaking. Essential symmetry breaking was first introduced by Lee and Head-Gordon [ Phys. Chem. Chem. Phys. 2019, 21, 4763-4778, 10.1039/C8CP07613H]. We utilized essential complex and time-reversal symmetry breaking with ph-AFQMC to compute the triplet–singlet energy gap in the TS12 set. We found statistically better performance of ph-AFQMC with complex-restricted orbitals than with spin-unrestricted orbitals. We then showed the utilization of essential spin symmetry breaking when computing the singlet–triplet gap of a known biradicaloid, C 36 . ph-AFQMC with spin-unrestricted Hartree–Fock (ph-AFQMC+UHF) fails catastrophically even with spin-projection and predicts no biradicaloid character. With approximate Brueckner orbitals obtained from regularized orbital-optimized second-order Møller–Plesset perturbation theory (κ-OOMP2), ph-AFQMC quantitatively captures strong biradicaloid character of C 36 . Lastly, we applied ph-AFQMC to the computation of the quintet–triplet gap in a model iron porphyrin complex where brute-force methods with a small active space fail to capture the triplet ground state. We show unambiguously that neither triplet nor quintet is strongly correlated using UHF, κ-OOMP2, and coupled-cluster with singles and doubles (CCSD) performed on UHF and κ-OOMP2 orbitals. There is no essential symmetry breaking in this problem. By virtue of this, we were able to perform UHF+ph-AFQMC reliably with a cc-pVTZ basis set and predicted a triplet ground state for this model geometry. The largest ph-AFQMC in this work correlated 186 electrons in 956 orbitals. Our work highlights the utility, scalability, and accuracy of ph-AFQMC with a single-determinant trial wave function with essential symmetry breaking for systems mainly dominated by dynamical correlation with little static correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

r 2 SCAN-D4: Dispersion corrected meta-generalized gradient approximation for general chemical applications

We combine a regularized variant of the strongly constrained and appropriately normed semilocal density functional with the latest generation semi-classical London dispersion correction. The resulting density functional approximation r 2 SCAN-D4 has the speed of generalized gradient approximations while approaching the accuracy of hybrid functionals for general chemical applications. We demonstrate its numerical robustness in real-life settings and benchmark molecular geometries, general main group and organo-metallic thermochemistry, and non-covalent interactions in supramolecular complexes and molecular crystals. Main group and transition metal bond lengths have errors of just 0.8%, which is competitive with hybrid functionals for main group molecules and outperforms them for transition metal complexes. The weighted mean absolute deviation (WTMAD2) on the large GMTKN55 database of chemical properties is exceptionally small at 7.5 kcal/mol. This also holds for metal organic reactions with an MAD of 3.3 kcal/mol. The versatile applicability to organic and metal–organic systems transfers to condensed systems, where lattice energies of molecular crystals are within the chemical accuracy (errors <1 kcal/mol).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗