SEARCH · Search NASA

Results for “ESPResSo++”

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.

12 records

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Fundamental physics with ESPRESSO: a new determination of the D/H ratio towards PKS1937-101

ABSTRACT Primordial abundances of light elements are sensitive to the physics of the early Universe and can directly constrain cosmological quantities, such as the baryon-to-photon ratio $\eta _{10}$, the baryon density, and the number of neutrino families. Deuterium is especially suited for these studies: its primordial abundance is sensitive and monotonically dependent on $\eta _{10}$, allowing an independent measurement of the cosmic baryon density that can be compared, for instance, against the Planck satellite data. The primordial deuterium abundance can be measured in high H i column density absorption systems towards distant quasars. We report here a new measurement, based on high-resolution ESPRESSO data, of the primordial D i abundance of a system at redshift $z \sim 3.572$, towards PKS1937-101. Using only ESPRESSO data, we find a D /H ratio of $2.638\pm 10^{-5}$, while including the available UVES data improves the precision, leading to a ratio of $2.608 \pm 10^{-5}$. The results of this analysis agree with those of the most precise existing measurements. We find that the relatively low column density of this system ($\log {N_{\rm H_I}/ {\rm cm}^{-2}}\sim 18$) introduces modelling uncertainties, which become the main contributor to the error budget.

Guarneri, Francesco (ORCID:0000000347409762)↗

DOE OSTI2025

Validation of the ESPRESSO wavelength calibration using iodine absorption cell spectra

ABSTRACT High-quality wavelength calibration is crucial for science cases like radial-velocity studies of exoplanets, the search for a possible variation of fundamental constants, and the redshift drift experiment. However, for state-of-the-art spectrographs, it has become difficult to verify the wavelength calibration on sky because no astrophysical source provides spectra with sufficiently stable or accurate wavelength information. We therefore propose to use iodine absorption cells to validate the wavelength calibration. Observing a bright and featureless star through the iodine cell emulates an astrophysical target with exactly known spectral features that can be analysed like any other science target, allowing to verify the wavelength calibration derived from the internal calibration sources and to identify systematics in the data processing. As demonstration, we temporarily installed an $\mathrm{I_2}$ absorption cell at ESPRESSO. Employing a full forward modelling approach of the $\mathrm{I_2}$ spectrum, including the instrumental line-spread function, we demonstrate wavelength calibration accuracy at the level of a few $\mathrm{m\,s^{-1}}$. We also show that wavelength measurements do depend on the geometry of the light-injection into the spectrograph fibers. This highlights the importance of probing exactly the same light path as science targets, something not possible with internal calibration sources alone. We also demonstrate excellent radial-velocity stability at the ${{< 20\, \mathrm{cm\,s^{-1}}}}$ level in a full end-to-end fashion, from sky to data product. Our study therefore showcases the great potential of absorption cells for the verification and long-term monitoring of the wavelength calibration as well as the unique insights they can provide.

Schmidt, Tobias M. (ORCID:0000000248337273)↗

DOE OSTI2026

RIKEN TRIP Magnets Database

This dataset contains ab-initio calculation results for the temperature-dependent anomalous Hall conductivity, the anomalous Nernst effect, and the Seebeck coefficient. All calculations are based on ab-inito Quantum Espresso (PWSCF v.6.3) + Wannier90 (v.3.0.0). The dependence on carrier doping is also calculated. For all calculations a ferromagnetic order has been assumed, which might not correspond to the true ground state of the system. Tabulated values for the magnetic moments and essential input files for Quantum Espresso are available for download as attachments. This project has been supported by the RIKEN Transformative Research Innovation Platform (TRIP), Use Case: Many-body Electron Systems.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Code modernization strategies for short-range non-bonded molecular dynamics simulations

Modern HPC systems are increasingly relying on greater core counts and wider vector registers. Thus, applications need to be adapted to fully utilize these hardware capabilities. One class of applications that can benefit from this increase in parallelism are molecular dynamics simulations. In this paper, we describe our efforts at modernizing the ESPResSo++ simulation package for molecular dynamics by restructuring its particle data layout for efficient memory accesses and applying vectorization techniques to benefit the calculation of short-range non-bonded forces, which results in an overall three times speedup and serves as a baseline for further optimizations. We also implement fine-grained parallelism for multi-core CPUs through HPX, a C++ runtime system which uses lightweight threads and an asynchronous many-task approach to maximize concurrency. Our goal is to evaluate the performance of an HPX-based approach compared to the bulk-synchronous MPI-based implementation. This requires the introduction of an additional layer to the domain decomposition scheme that defines the task granularity. On spatially inhomogeneous systems, which impose a corresponding load-imbalance in traditional MPI-based approaches, we demonstrate that by choosing an optimal task size, the efficient work-stealing mechanisms of HPX can overcome the overhead of communication resulting in an overall 1.4 times speedup compared to the baseline MPI version.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

SMATool: Strength of materials analysis toolkit

The study of the strength of materials is a cornerstone in material science and engineering, playing a critical role in shaping the progress and application of materials in diverse industrial sectors. The strength of a material is meticulously examined to understand the behavior of the material under different stress conditions and environments, thereby guiding material selection and structural design. Herein, we introduce the SMATool, a computational toolkit for the efficient calculation and analysis of material strength at both zero and finite temperatures for 3D, 2D, 1D, and tubular 2D-based nanostructures and nanotubes, as well as 1D nanoribbons. The toolkit is capable of calculating tensile, shear, ultimate, yield, and indentation (Vickers' hardness) strengths in various dimensions, as well as the energy storage capacity. We conducted several calculations both at zero and finite temperatures to validate the accuracy and reliability of the developed software. Here, the results show that the SMATool package provides accurate predictions that align with existing data on material strength. SMATool integrates seamlessly with widely used electronic structure codes like VASP and Quantum Espresso, providing a user-friendly interface catering to academic researchers and industry professionals.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Low-energy pathways lead to self-healing defects in CsPbBr 3

Self-regulation of free charge carriers in perovskites via Schottky defect formation has been posited as the origin of the well-known defect-tolerance of metal halide perovskite materials. Understanding the mechanisms of self-regulation promises to lead to the fabrication of better performing solar cell materials with higher efficiencies. We investigate many such mechanisms here for CsPbBr 3 , a popular representative of a more commercially viable all-inorganic metal halide perovskite. We investigate different atomic-level mechanisms and pathways of the diffusion and recombination of neutral and charged interstitials and vacancies (Schottky pairs) in CsPbBr 3 . We use nudged elastic band calculations and ab initio-derived pseudopotentials within quantum ESPRESSO to determine energies of formation and migration and hence the activation energies for these defects. While halide vacancies are known to exhibit low formation energies, the migration of interstitials is less studied. Our calculations uncover interstitial defect pathways capable of producing an activation energy at, or below, the single experimental value of 0.53 eV observed for the slow, temperature-dependent recovery of light-induced conductivity in bulk CsPbBr 3 . Our work reveals the existence of a low-energy diffusion pathway involving a concerted “domino effect” of interstitials, with the net result that interstitials can diffuse more readily over longer distances than expected. This observation suggests that defect self-healing can be promoted if the “domino effect” strategy can be engaged.

14 SOLAR ENERGY↗

DOE OSTI2024

TOI–757 b: an eccentric transiting mini–Neptune on a 17.5–d orbit

ABSTRACT We report the spectroscopic confirmation and fundamental properties of TOI$-$757 b, a mini$-$Neptune on a 17.5$-$d orbit transiting a bright star ($V\, =\, 9.7$ mag) discovered by the TESS mission. We acquired high$-$precision radial velocity measurements with the HARPS, ESPRESSO, and PFS spectrographs to confirm the planet detection and determine its mass. We also acquired space$-$borne transit photometry with the CHEOPS space telescope to place stronger constraints on the planet radius, supported with ground$-$based LCOGT photometry. WASP and KELT photometry were used to help constrain the stellar rotation period. We also determined the fundamental parameters of the host star. We find that TOI$-$757 b has a radius of $R_{\mathrm{p}} = 2.5 \pm 0.1 R_{\oplus }$ and a mass of $M_{\mathrm{p}} = 10.5^{+2.2}_{-2.1} M_{\oplus }$, implying a bulk density of $\rho _{\text{p}} = 3.6 \pm 0.8$ g cm$^{-3}$. Our internal composition modelling was unable to constrain the composition of TOI$-$757 b, highlighting the importance of atmospheric observations for the system. We also find the planet to be highly eccentric with e = 0.39$^{+0.08}_{-0.07}$, making it one of the very few highly eccentric planets among precisely characterized mini$-$Neptunes. Based on comparisons to other similar eccentric systems, we find a likely scenario for TOI$-$757 b’s formation to be high eccentricity migration due to a distant outer companion. We additionally propose the possibility of a more intrinsic explanation for the high eccentricity due to star$-$star interactions during the earlier epoch of the Galactic disc formation, given the low metallicity and older age of TOI$-$757.

Alqasim, A. (ORCID:0000000151025505)↗

DOE OSTI2024

A possible misaligned orbit for the young planet AU Mic c

ABSTRACT The AU Microscopii planetary system is only 24 Myr old, and its geometry may provide clues about the early dynamical history of planetary systems. Here, we present the first measurement of the Rossiter–McLaughlin effect for the warm sub-Neptune AU Mic c, using two transits observed simultaneously with the European Southern Observatory's (ESO's) Very Large Telescope (VLT)/Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations (ESPRESSO), CHaracterising ExOPlanet Satellite (CHEOPS), and Next-Generation Transit Survey (NGTS). After correcting for flares and for the magnetic activity of the host star, and accounting for transit-timing variations, we find the sky-projected spin–orbit angle of planet c to be in the range $\lambda _{\mathrm{c}}=67.8_{-49.0}^{+31.7}$ degrees (1$\sigma$). We examine the possibility that planet c is misaligned with respect to the orbit of the inner planet b ($\lambda _{\mathrm{b}}=-2.96_{-10.30}^{+10.44}$), and the equatorial plane of the host star, and discuss scenarios that could explain both this and the planet’s high density, including secular interactions with other bodies in the system or a giant impact. We note that a significantly misaligned orbit for planet c is in some degree of tension with the dynamical stability of the system, and with the fact that we see both planets in transit, though these arguments alone do not preclude such an orbit. Further observations would be highly desirable to constrain the spin–orbit angle of planet c more precisely.

Yu, H. (ORCID:0000000209716078)↗

DOE OSTI2026

matsim-agents v1.0

matsim-agents is a multi-agent AI framework for atomistic materials simulation and discovery. It orchestrates large language models (LLMs), machine-learned interatomic potentials (MLIPs), and DFT codes into a single agentic loop running on laptops and DOE leadership-class supercomputers. MULTI-AGENT ORCHESTRATION A LangGraph state machine with three nodes: a Planner that converts a natural-language research objective into structured tasks; an Executor that dispatches atomistic tools and loops until the queue is empty; and an Analyst that summarizes results into a human-readable report. State is checkpointed after every step and human-in-the-loop gates can be inserted at any edge. HYPOTHESIS-DRIVEN DISCOVERY CHAT An interactive REPL (matsim-agents chat) that couples LLM dialogue with atomistic simulation. Chemical formulas are automatically detected in conversation turns and trigger a full crystal-phase exploration: structure generation → relaxation → stability scoring → result injection back into the conversation, creating a closed hypothesis-refinement loop. CRYSTAL PHASE ENUMERATION Given a composition, the phase explorer enumerates prototypes by stoichiometry: elemental (fcc/bcc/hcp/sc/diamond), binary 1:1 (rocksalt/CsCl/zincblende/ wurtzite/fluorite/rutile), ternary 1:1:3 (cubic perovskite), ternary 1:2:4 (perovskite + spinel), quaternary 1:1:2:6 (Fm-3m double perovskite). 2-D prototypes (graphene, h-BN, MoS2 2H/1T) and multilayer stacking are also supported via --include-2d and --num-layers. SUPERCELL GENERATION AND SITE DECORATION Auto-tiling to a minimum atom count (--min-atoms), explicit NxNxN tiling (--supercell), symmetry-distinct site decorations (--n-orderings), and isotropic lattice-scale sweeps (--lattice-scales) for volume bracketing. MLFF RELAXATION AND STABILITY SCORING HydraGNN (multi-headed GNN) drives structure relaxation via ASE with FIRE, BFGS, or BFGSLineSearch. Stability output: delta-E/atom ranking across phases and a max-residual-force dynamical-stability proxy. Other MLIPs (MACE, NequIP, Orb) can be plugged in through the same interface. DFT BACKENDS Quantum ESPRESSO pw.x and VASP 6.6 are first-class labellers. Both have validated GPU builds and SLURM/PBS launchers for three DOE platforms: Frontier (AMD MI250X, ROCm), Aurora (Intel PVC, oneAPI), Perlmutter (NVIDIA A100, CUDA). QE produces ~100 binaries (pw.x, ph.x, epw.x, ...). VASP supports scf, relax, vc-relax, and vc-relax-shape run types. ACTIVE-LEARNING LOOP matsim-agents al run CONFIG.yaml drives an iterative HydraGNN-DFT loop: MD generates candidates → ensemble/MC-dropout uncertainty selects the most informative → DFT labels them in parallel inside one allocation → dataset grows → HydraGNN retrains → repeat. DFT backend is a single YAML toggle (dft.backend: vasp | qe). LLM-generated seed structures are supported (no curated POSCAR library needed). Config uses ${VAR}, ${VAR:-default}, ${VAR:?msg} shell-style substitution for cross-user/cross-site portability. LLM BACKENDS Ollama (local, default), vLLM (HPC multi-GPU serving), OpenAI, Anthropic, HuggingFace Transformers+Accelerate. Selected at runtime via flag or env var with no code changes. HPC PORTABILITY Same Python entry points run on Frontier (ROCm 7.2), Aurora (oneAPI), and Perlmutter (CUDA 12). DFT and ML stacks are never co-loaded in the same shell; they couple through the scheduler and filesystem. Advanced multi-node launchers (serve, discovery-chat, single-relaxation, active-learning, QE warm-start) are provided for all three platforms. CODABENCH COMPETITION BUNDLE A self-contained benchmark: 159 atomistic test structures across 11 material classes, 5 tasks (formation energy, forces, ML relaxation, AI-DFT relaxation, phase stability ranking), public/private leaderboard split (30/70), and four ready-to-run baselines: MACE-MP-0, HydraGNN, UMA, AllScAIP.

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗

DOE OSTI2024

HTESP (High-throughput electronic structure package): A package for high-throughput ab initio calculations

High-throughput ab initio calculations are the indispensable parts of data-driven discovery of new materials with desirable properties, as reflected in the establishment of several online material databases. The accumulation of extensive theoretical data through computations enables data-driven discovery by constructing machine learning and artificial intelligence models to predict novel compounds and forecast their properties. Efficient usage and extraction of data from these existing online material databases can accelerate the next stage materials discovery that targets different and more advanced properties, such as electron–phonon coupling for phonon-mediated superconductivity. However, extracting data from these databases, generating tailored input files for different ab initio calculations, performing such calculations, and analyzing new results can be demanding tasks. Here, in this work, we introduce a software package named “HTESP” (High-Throughput Electronic Structure Package) written in Python and Bash languages, which automates the entire workflow including data extraction, input file generation, calculation submission, result collection and plotting. Our HTESP will help speed up future computational materials discovery processes.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Violation of the Cauchy–Born rule in multi-principal element alloys

Multi-principal element alloys are a novel class of materials that are formed by combining multiple elements in high concentrations and show exceptional properties compared to conventional alloys. These alloys have high configurational entropy due to inherent atomic disorder. The Cauchy–Born rule is a popular homogenization method for linking atomistic to continuum length scales. In this Letter, we use ab initio density functional theory calculations to report that the Cauchy–Born rule, which holds in conventional alloys in the absence of defects, is not valid in multi-principal element alloys. Further, the violation of the Cauchy–Born rule in these alloys is attributed to the presence of atomic disorder due to which the deformations are inhomogeneous. Our results also show the deviations in deformation are related to the magnitude of the stretch and shear by power laws.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Li1−xNiO2 Many-body DMC Benchmark Dataset

The dataset contains all numerical data generated in support of the manuscript “Many‑body Benchmark of Electronic Charge and Spin Densities for Li1–xNiO2” (Journal of Chemical Theory and Computation, DOI: 10.1021/acs.jctc.5c02097, URL: https://pubs.acs.org/doi/10.1021/acs.jctc.5c02097). The materials included in this repository are: 1. Data files used to produce all figures and tables in the main manuscript and supporting information. 2. Benchmark density‑functional theory (DFT) datasets used for the charge‑ and spin‑density analyses. 3. Reference many‑body diffusion Monte Carlo (DMC) calculations and associated input/output files.

36 MATERIALS SCIENCE↗