Search NASA⌕ Search

SEARCH · Search NASA

Results for “Periodic orbits”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Dynamical Avenues for Mercury’s Origin. I. The Lone Survivor of a Primordial Generation of Short-period Protoplanets

The absence of planets interior to Mercury continues to puzzle terrestrial-planet formation models, particularly when contrasted with the relatively high derived occurrence rates of short-period planets around Sun-like stars. Recent work proposed that the majority of systems hosting hot super-Earths attain their orbital architectures through an epoch of dynamical instability after forming in quasi-stable, tightly packed configurations. Isotopic evidence seems to suggest that the formation of objects in the super-Earth-mass regime is unlikely to have occurred in the solar system as the terrestrial-forming disk is thought to have been significantly mass deprived starting around 2 Myr after the formation of calcium-aluminum-rich inclusions—a consequence of either Jupiter’s growth or an intrinsic disk feature. Nevertheless, terrestrial-planet formation models and high-resolution investigations of planetesimal dynamics in the gas-disk phase occasionally find that quasi-stable protoplanets with mass comparable to that of Mars emerge in the vicinity of Mercury’s modern orbit. In this paper, we investigate whether it is possible for a primordial configuration of such objects to be cataclysmically destroyed in a manner that leaves Mercury behind as the sole survivor without disturbing the other terrestrial worlds. We use numerical simulations to show that this scenario is plausible. In many cases, the surviving Mercury analog experiences a series of erosive impacts, thereby boosting its Fe/Si ratio. A caveat of our proposed genesis scenario for Mercury is that Venus typically experiences at least one late giant impact.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Orbits and Masses of Binaries from Speckle Interferometry at SOAR

We present results from Speckle inteferometric observations of 15 visual binaries and one double-line spectroscopic binary, carried out with the HRCam Speckle camera of the SOAR 4.1 m telescope. These systems were observed as a part of an on-going survey to characterize the binary population in the solar vicinity, out to a distance of 250 pc. We obtained orbital elements and mass sums for our sample of visual binaries. The orbits were computed using a Markov Chain Monte Carlo algorithm that delivers maximum likelihood estimates of the parameters, as well as posterior probability density functions that allow us to evaluate their uncertainty. Their periods cover a range from 5 yr to more than 500 yr; and their spectral types go from early A to mid M, implying total system masses from slightly more than 4M{sub ⊙} down to 0.2M {sub ⊙}. They are located at distances between approximately 12 and 200 pc, mostly at low Galactic latitude. For the double-line spectroscopic binary YSC8, we present the first combined astrometric/radial-velocity orbit resulting from a self-consistent fit, leading to individual component masses of 0.897 ± 0.027 M {sub ⊙} and 0.857 ± 0.026 M {sub ⊙}; and an orbital parallax of 26.61 ± 0.29 mas, which compares very well with the Gaia DR2 trigonometric parallax (26.55 ± 0.27 mas). In combination with published photometry and trigonometric parallaxes, we place our objects on an H-R diagram and discuss their evolutionary status. We also present a thorough analysis of the precision and consistency of the photometry available for them.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Extensive Secondary Cratering From the InSight Sol 1034a Impact Event

Abstract Impact cratering is one of the fundamental processes throughout the history of the Solar System. The formation of new impact craters on planetary bodies has been observed with repeat images from orbiting satellites. However, the time gap between images is often large enough to preclude detailed analysis of smaller‐scale features such as secondary impact craters, which are often removed or buried over a short time period. Here we use a seismic event detected on Mars by the NASA InSight mission to investigate secondary cratering at a new impact crater. We strengthen the case that the seismic event that occurred on Sol 1034 (S1034a) is the result of a new impact cratering event. Using the exact timing of this event from InSight, we investigated the resulting new impact crater in orbital image data. The S1034a impact crater is approximately 9 m in diameter but is responsible for over 900 secondary impact events in the form of low albedo spots that are located at distances of up to almost 7 km from the primary crater. We suggest that the low albedo spots formed from relatively low energy ejecta, with individual ejecta block velocities less than 200 m s −1 . We estimate that the low albedo spots, the main evidence of secondary impact processes at this new impact event, fade within 200–300 days after formation.

Grindrod, P. M. [Natural History Museum London UK]↗
DOE OSTI2023

2D Exotic Optical Lattice via a Digital‐Coding Circular Airy Beam

Optical lattices have been widely used from classical to quantum physics. The tunable and scalable fabrication of lattices would be of great significance in lattice‐based multipartite applications. This work demonstrates first that a circular Airy beam (CAB), which has the peculiar properties of self‐healing and abrupt autofocusing, can be used to generate two‐dimensional (2D) optical lattices in propagation when encoded by a programmable spatial mask, resulting in the formation of large‐scale and tunable optical lattices with both axis and axial symmetry, and even high‐orbital kaleidoscope shapes. The efficient diffraction of CAB during the spatial crosstalk with the mask enables the realization of tunable lattices with rich periodicity and complexity. The study shows a flexible method to manipulate lattices with large‐scale and versatile structures for potential applications in integrated and scalable optical and photonic devices.

Sun, Peisheng↗
DOE OSTI2023

Density streams in the disc winds of Classical T Tauri stars

ABSTRACT Spectral and photometric variability of the Classical T Tauri stars RY Tau and SU Aur from 2013 to 2022 is analysed. We find that in SU Aur the H α line’s flux at radial velocity RV = −50 ± 7 km s−1 varies with a period P = 255 ± 5 d. A similar effect previously discovered in RY Tau is confirmed with these new data: P = 21.6 d at RV = −95 ± 5 km s. In both stars, the radial velocity of these variations, the period, and the mass of the star turn out to be related by Kepler’s law, suggesting structural features on the disc plane orbiting at radii of 0.2 au in RY Tau and 0.9 au in SU Aur, respectively. Both stars have a large inclination of the accretion disc to the line of sight – so that the line of sight passes through the region of the disc wind. We propose there is an azimuthal asymmetry in the disc wind, presumably in the form of ‘density streams,’ caused by substructures of the accretion disc surface. These streams cannot dissipate until they go beyond the Alfven surface in the disc’s magnetic field. These findings open up the possibility to learn about the structure of the inner accretion disc of CTTS on scales less than 1 au and to reveal the orbital distances related to the planet’s formation.

Petrov, P. P. (ORCID:0000000155514963)↗
DOE OSTI2020

Block orbital-selective Mott insulators: A spin excitation analysis

We present a comprehensive study of the spin excitations—as measured by the dynamical spin structure factor S(q,ω)—of the so-called block-magnetic state of low-dimensional orbital-selective Mott insulators. We realize this state via both a multi-orbital Hubbard model and a generalized Kondo-Heisenberg Hamiltonian. Due to various competing energy scales present in the models, the system develops periodic ferromagnetic islands of various shapes and sizes, which are antiferromagnetically coupled. The 2×2 particular case was already found experimentally in the ladder material BaFe2Se3 that becomes superconducting under pressure. Here we discuss the electronic density as well as Hubbard and Hund coupling dependence of S(q,ω) using density matrix renormalization group method. Several interesting features were identified: (1) An acoustic (dispersive spin-wave) mode develops. (2) The spin-wave bandwidth establishes a new energy scale that is strongly dependent on the size of the magnetic island and becomes abnormally small for large clusters. (3) Optical (dispersionless spin excitation) modes are present for all block states studied here. In addition, a variety of phenomenological spin Hamiltonians have been investigated but none matches entirely our results that were obtained primarily at intermediate Hubbard U strengths. Our comprehensive analysis provides theoretical guidance and motivation to crystal growers to search for appropriate candidate materials to realize the block states, and to neutron scattering experimentalists to confirm the exotic dynamical magnetic properties unveiled here, with a rich mixture of acoustic and optical features.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Accurate frozen core approximation for all-electron density-functional theory

We implement and benchmark the frozen core approximation, a technique commonly adopted in electronic structure theory to reduce the computational cost by means of mathematically fixing the chemically inactive core electron states. The accuracy and efficiency of this approach are well controlled by a single parameter, the number of frozen orbitals. Explicit corrections for the frozen core orbitals and the unfrozen valence orbitals are introduced, safeguarding against seemingly minor numerical deviations from the assumed orthonormality conditions of the basis functions. A speedup of over twofold can be achieved for the diagonalization step in all-electron density-functional theory simulations containing heavy elements, without any accuracy degradation in terms of the electron density, total energy, and atomic forces. This is demonstrated in a benchmark study covering 103 materials across the Periodic Table and a large-scale simulation of CsPbBr3 with 2560 atoms. Our study provides a rigorous benchmark of the precision of the frozen core approximation (sub-meV per atom for frozen core orbitals below −200 eV) for a wide range of test cases and for chemical elements ranging from Li to Po. The algorithms discussed here are implemented in the open-source Electronic Structure Infrastructure software package.

Yu, Victor Wen-zhe↗
DOE OSTI2022

A detailed study of the barium central star of the planetary nebula Abell 70

ABSTRACT We present a detailed study of the barium star at the heart of the planetary nebula Abell 70. Time-series photometry obtained over a period of more than 10 yr demonstrates that the barium-contaminated companion is a rapid rotator with temporal variability due to spots. The amplitude and phasing of the photometric variability change abruptly; however, there is no evidence for a change in the rotation period (P = 2.06 d) over the course of the observations. The co-addition of 17 high-resolution spectra obtained with Ultraviolet and Visual Échelle Spectrograph mounted on the Very Large Telescope allows us to measure the physical and chemical properties of the companion, confirming it to be a chromospherically active, late G-type sub-giant with more than +1 dex of barium enhancement. We find no evidence of radial velocity variability in the spectra, obtained over the course of approximately 130 d with a single additional point some 8 yr later, with the radial velocities of all epochs approximately −10 km s −1 from the previously measured systemic velocity of the nebula. This is perhaps indicative that the binary has a relatively long period (P ≳ 2 yr) and high eccentricity (e ≳ 0.3), and that all the observations were taken around radial velocity minimum. However, unless the binary orbital plane is not aligned with the waist of the nebula or the systemic velocity of the binary is not equal to the literature value for the nebula, this would imply an unfeasibly large mass for the nebular progenitor.

Jones, David (ORCID:0000000339475946)↗
DOE OSTI2022

Three Dimensional Electromagnetic High Frequency Convex Cavity Scars

This report examines the localization of high frequency electromagnetic fields in general three-dimensional convex walled cavities along periodic paths between opposing sides of the cavity. The report examines the three-dimensional case where the mirrors at the end of the orbit have two different radii of curvature. The cases where these orbits lead to unstable localized modes are known as scars.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Periodic Gamma-Ray Modulation of the Blazar PG 1553+113 Confirmed by Fermi-LAT and Multiwavelength Observations

A 2.1 yr periodic oscillation of the gamma-ray flux from the blazar PG 1553+113 has previously been tentatively identified in ∼7 yr of data from the Fermi Large Area Telescope. After 15 yr of Fermi sky-survey observations, doubling the total time range, we report >7 cycle gamma-ray modulation with an estimated significance of 4σ against stochastic red noise. Independent determinations of oscillation period and phase in the earlier and the new data are in close agreement (chance probability <0.01). Pulse timing over the full light curve is also consistent with a coherent periodicity. Multiwavelength new data from Swift X-Ray Telescope, Burst Alert Telescope, and UVOT, and from KAIT, Catalina Sky Survey, All-Sky Automated Survey for Supernovae, and Owens Valley Radio Observatory ground-based observatories as well as archival Rossi X-Ray Timing Explorer satellite-All Sky Monitor data, published optical data of Tuorla, and optical historical Harvard plates data are included in our work. Optical and radio light curves show clear correlations with the gamma-ray modulation, possibly with a nonconstant time lag for the radio flux. We interpret the gamma-ray periodicity as possibly arising from a pulsational accretion flow in a sub-parsec binary supermassive black hole system of elevated mass ratio, with orbital modulation of the supplied material and energy in the jet. Other astrophysical scenarios introduced include instabilities, disk and jet precession, rotation or nutation, and perturbations by massive stars or intermediate-mass black holes in polar orbit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

cyquandyn

CyQuanDyn is a collection of codes developed for modelling the phonon-induced carrier pumping in ZrTe5 system. Codes in the collection include 1) TLS: time-dependent two-orbital spinless fermion Bogoliubov-de Gennes Hamiltonian simulation, 2) TDSE: ab initio time-dependent Schrödinger equation simulation of the ZrTe5 system under the periodic modulation of A1g Raman mode. Associated publication is T. Jiang, P. P. Orth, L. Luo, L.-L. Wang, F. Zhang, C.-Z. Wang, J. Zhao, K.-M. Ho, J. Wang, and Y.-X. Yao, Ab-Initio Simulations of Coherent Phonon-Induced Pumping of Carriers in ZrTe5. [arXiv:2304.08449 (2023).](https://arxiv.org/abs/2304.08449)

Yao, Yongxin↗
DOE OSTI2025

Datasets for Custom-trained Machine-learning Interatomic Potentials: Nitric Acid Aqueous Solution

This dataset was generated using an iterative active learning strategy with the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials (MLIPs) for aqueous nitric acid. Each active-learning cycle consisted of three stages: (1) training, (2) exploration, and (3) labeling. The initial training set comprised approximately 800 randomly selected configurations from a previous study by Lewis et al. (https://doi.org/10.1021/jp205510q), which investigated nitric acid solutions at 2, 3, 4, and 5 mol/L. For all configurations, single-point calculations of atomic forces and total energies were performed at the quantum density functional theory BLYP-D2 and PBE-D3 levels of theory using the CP2K Quickstep module. Valence electrons were treated explicitly, while core electrons on all atoms were represented by norm-conserving Goedecker–Teter–Hutter (GTH) pseudopotentials. Long-range dispersion interactions were accounted for using Grimme dispersion corrections. Wave functions were expanded in a mixed Gaussian-and-plane-wave scheme using TZV2P-MOLOPT basis sets for all elements and an 800 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent field convergence was accelerated using orbital transformation and Direct Inversion in the Iterative Subspace, with a convergence threshold of 10^{-6}. All single-point calculations were carried out in periodic orthorhombic cells whose dimensions match those of the molecular configurations sampled from earlier trajectories. The CELL_REF keyword in CP2K was used to define a fixed reference cell, ensuring consistency in the reference data used for MLIP training, particularly when cell fluctuations are present in NpT simulations. The resulting high-fidelity energies and forces constitute the ground-truth labels used to train the MLIPs contained in this dataset.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2026

Custom-trained Machine-learning Interatomic Potentials: ZnCl2 Aqueous Solution

This dataset was generated using an iterative active-learning strategy implemented in the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials for aqueous ZnCl2 solutions. Each active-learning cycle consisted of three stages: training, exploration, and labeling. The initial training set combined configurations generated in this work from enhanced-sampling ab initio molecular dynamics simulations with configurations from a previously reported neural-network-potential study of aqueous ZnCl2. The enhanced-sampling ab initio molecular dynamics simulations involved Zn–Cl separation and the chloride coordination number around Zn²? as collective variables. These configurations served as the seed dataset. Subsequent active-learning cycles expanded the training set by identifying and labeling configurations that were poorly represented by the current models, thereby improving coverage of ion-association states and changes in local coordination and charge-state environments relevant to the solution free-energy landscape. For all selected configurations, single-point calculations of the total energies and atomic forces were performed within density functional theory using the CP2K Quickstep module. Reference calculations employed the revPBE-D3 and r2SCAN exchange-correlation functionals. Motivated by recent work on aqueous Zn²?, the main revPBE calculations omitted D3 dispersion contributions involving Zn²?, while retaining the D3 correction for water and chloride. For comparison, fully dispersion-corrected revPBE-D3 reference calculations were also performed, with D3 applied to all species, including Zn²?. Valence electrons were treated explicitly, while core electrons were represented using norm-conserving Goedecker–Teter–Hutter pseudopotentials. The wave functions were expanded using the mixed Gaussian-and-plane-wave scheme with TZV2P-MOLOPT basis sets for all elements and a 600 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent-field convergence was accelerated using the orbital-transformation and Direct Inversion in the Iterative Subspace algorithms, with a convergence threshold of 10?6. All single-point calculations were performed in periodic orthorhombic cells. The CELL_REF keyword in CP2K was used to define a fixed reference cell with a box length of 25 Å. This treatment ensured a consistent reference for configurations extracted from NpT trajectories with fluctuating cell dimensions. The resulting DFT energies and atomic forces constitute the ground-truth labels used to train the MLIPs. The resulting MLIP was trained for aqueous ZnCl2 solutions spanning concentrations from 0 to 30 molal and a broad pH range, from strongly acidic to strongly basic conditions. Representative examples of configurations included in the MLIP training dataset are provided below. These include 1) Representative configurations from the dataset labeled at the revPBE-D3 level, with D3 dispersion interactions involving Zn2+ excluded (revPBE-wo-D3). 2) Representative configurations from the dataset labeled at the fully dispersion-corrected revPBE-D3 level, with D3 interactions applied to all species, including Zn2+ (revPBE-D3). 3) Representative configurations from the dataset labeled at the r2SCAN level of theory (r2SCAN).

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2023

Band structures and Z 2 invariants of two-dimensional transition metal dichalcogenide monolayers from fully relativistic Dirac-Kohn-Sham theory using Gaussian-type orbitals

Two-dimensional (2D) materials exhibit a wide range of remarkable phenomena, many of which owe their existence to the relativistic spin-orbit coupling (SOC) effects. To understand and predict properties of materials containing heavy elements, such as the transition-metal dichalcogenides (TMDs), relativistic effects must be taken into account in first-principles calculations. We present an all-electron method based on the four-component Dirac Hamiltonian and Gaussian-type orbitals (GTOs) that overcomes complications associated with linear dependencies and ill-conditioned matrices that arise when diffuse functions are included in the basis. Until now, there has been no systematic study of the convergence of GTO basis sets for periodic solids either at the nonrelativistic or the relativistic level. Here, in this study, we provide such a study of relativistic band structures of the 2D TMDs in the hexagonal (2H), tetragonal (1T), and distorted tetragonal (1T') structures, along with a discussion of their SOC-driven properties (Rashba splitting and Z 2 topological invariants). We demonstrate the viability of our approach even when large basis sets with multiple basis functions involving various valence orbitals (denoted triple- and quadruple-ζ) are used in the relativistic regime. Our method does not require the use of pseudopotentials and provides access to all electronic states within the same framework. Our study paves the way for direct studies of material properties, such as the parameters in spin Hamiltonians, that depend heavily on the electron density near atomic nuclei where relativistic and SOC effects are the strongest.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Preserving Symmetry and Degeneracy in the Localized Orbital Scaling Correction Approach

Symmetry is a fundamental concept that plays a critical role in many chemical and physical phenomena and processes, which highlights the importance of theoretical methods to correctly handle symmetry. The recently developed localized orbital scaling correction (LOSC1) shows great improvement on the description of band gaps, photoemission spectra, and dissociation limits of cationic species. However, issues remain with LOSC1 in dealing with the symmetry and degeneracy of electronic states, which are also relevant to other methods using localization. In this work, we utilize a new method that deals with the physical-space and the energy-space localization on an equal footing. The resulting localized orbitals, i.e., orbitalets, are able to maintain more symmetry and the desired state degeneracy, which is important in calculating the electronic structure of both molecules and periodic bulk systems. Furthermore, the curvature matrix is redefined to improve potential energy curves for systems with stretched bonds, while retaining the correct dissociation limits. This new approach, termed LOSC2, includes only two fitting parameters. It maintains accuracy similar to that of LOSC1 over many properties, while overcoming LOSC1’s deficiencies in symmetry and degeneracy. Our tests have shown that LOSC2 orbitalets possess the full- or subgroup of molecular symmetry if allowed, which preserves the state degeneracy. Tests on differently sized planar annulenes, odd-numbered allenes, and triphenylene again verify that LOSC2 is able to maintain the state degeneracy, while LOSC1 cannot. Finally, all the tests demonstrate the advantage of LOSC2 in the calculation of molecular systems and its potential for application to periodic bulk systems.

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

Woody Feedstock 2022 State of Technology Report

The U.S. Department of Energy promotes production of advanced liquid transportation fuels from lignocellulosic biomass by funding fundamental and applied research that advances the state of technology (SOT). As part of its involvement in this mission, Idaho National Laboratory completes an annual SOT report for n th -plant and 1 st -plant woody biomass feedstock logistics. The purpose of the SOT is to provide the status of feedstock supply system technology development for woody biomass to biofuels relative to technical targets and cost goals from specific design cases, based on data and experimental results. Conventional feedstock supply systems need to be modified to meet the demands of conversion pathways, specifically to have the ability to adjust the quality of the raw biomass materials. Advanced systems incorporate innovative methods of material handling, preprocessing and supply chain configuration. In advanced designs, variability of the raw biomass can be reduced to produce feedstocks of a uniform format, moving toward biomass commoditization. Against this backdrop, the 2022 Woody SOT for low-ash woody feedstocks utilizes feedstock fractionation by incorporating technologies that can separate the biomass into its anatomical fractions (wood, bark, needle, and extrinsic ash) to reduce impurities and attempt to maximize the retention of usable fractions that satisfy downstream quality considerations. By using a series of air classification steps, this strategy can reduce the extrinsic ash in forest residues, separate out a majority of the incoming needles (which can be supplied to alternate markets), and maximize the retention of whitewood in the usable fraction. The fractionated forest residues are then mixed with clean-pine chips in a 50-50 blend to prepare the feedstock for the desired conversion pathway. The n th -plant analysis estimated the delivered cost for the feedstock at $\$$69.23/dry ton (2016$\$$) which represents a $\$$6.64/dry ton decrease compared to the cost estimate of the 2021 Woody SOT supply system for low-ash woody feedstocks. The quality requirements in the 2022 Woody SOT were identical to those of the 2021 Woody SOT at = 1.00 wt % ash and = 50.51 wt% carbon. The cost savings derive primarily from reductions in dry matter losses during air classification. The GHG emissions for the n th -plant analysis were estimated at 178.39 kg CO2e/dry ton compared to 178.71 kg CO 2 e/dry ton in the 2021 Woody SOT, a decrease of 0.32 kg CO2e/dry ton. The small change stems from an increase in emissions attributed to preprocessing and slightly larger savings in emissions from transportation. In the 1 st -plant analysis of the 2022 Woody SOT system, the average throughput was estimated to be approximately 2,128 dry tons/day or 96.51% of the name plate capacity. During the simulation the daily throughput ranged from 1,090 dry tons/day to 2,200 dry tons/day, or 49.43% to 99.75% of the daily nameplate capacity. After the year of operation 722,403 tons of processed feedstock were produced in total without regard to quality considerations (99.64% of the annual nameplate capacity). The variability in throughput was primarily caused by equipment failures in the system. Regular failures, downtime caused by routine maintenance per manufacturer guidelines, contributed to a majority 62.50% of failures and 62.60% of downtime. Failures due to wear were the other cause of disruption within the system, impacting the rotary shear and orbital screen and accounting for 37.50% of the failures and 37.40% of the total downtime. Ultimately the system was on stream for 87.84% during the simulation period, which is only 2.16 percentage points below the nth-plant assumption for on-stream time. The production cost of the system averaged $\$$71.66/dry ton. The costs ranged from a minimum of $\$$71.23/dry ton to a maximum of $\$$2,115.30/dry ton. When dry matter losses (disposed low-quality fractions as well as other losses such as in grinders) were considered the costs increased to an average of $\$$75.11/dry ton with a minimum of $\$$74.69/dry ton and a maximum of $\$$2,136.86/dry ton...

09 BIOMASS FUELS↗
DOE OSTI2020

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

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

Extended Lagrangian Born–Oppenheimer molecular dynamics for orbital-free density-functional theory and polarizable charge equilibration models

We report extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) is formulated for orbital-free Hohenberg–Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born–Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗