Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Virtualization”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Scalable Generation of High-fidelity Synthetic Population Ensembles

Used within social simulations, synthetic population ensembles enable uncertainty quantification (UQ) methods for obtaining more robust model inference and prediction. A synthetic population ensemble is a series of plausible virtual reconstructions of an area’s population at the granularity of people and residences, generated stochastically to preserve privacy of the source population survey’s respondents. In this paper, we demonstrate the production of large synthetic population ensembles for the U.S. via Oak Ridge National Laboratory’s UrbanPop framework to support modeling of high spatial resolution energy affordability metrics from nationwide social surveys in collaboration with the fusionACS project. The study involves two scenarios: creating ensembles for (1) 17 U.S. metropolitan areas in 2019 and (2) full U.S. Census Divisions in 2023, with each scenario consisting of 41 population instances (a base realization and 40 replicates). To accomplish this task at scale, we configured an integrated system within a research cloud, comprised of virtual containerizations, GPU-enhanced functionality, and orchestrated deployments of UrbanPop’s maturing Likeness Python ecosystem. Results demonstrate we maintained high-fidelity approximations of residential totals by areas of interest and the demographic characteristics of neighborhoods while reducing manual workflow burdens. Finally, we discuss plans to fine-tune and further develop our automated workflows for truly distributed job orchestration to increase computational efficiency, as well as provide an outlook for broadening applications of the ensembles.

Cluster computing↗
DOE OSTI2022

Detection of surface water temperature variations of Mongolian lakes benefiting from the spatially and temporally gap-filled MODIS data

Lakes provide critical water resources for human activities and ecosystems, particularly in the Mongolian Plateau (MP), which is characterized by a dry climate and a harsh environment. As a region that is sensitive to anthropogenic warming, tracking lake surface water temperature (LSWT) changes in Mongolian lakes is crucial for understanding the consequences of a warming climate on lake ecosystems. However, the long-term monitoring of LSWT is restricted by the spatiotemporal gaps in the raw imagery of remote sensing-based land surface temperature (LST), e.g., the commonly used Moderate Resolution Imaging Spectroradiometer (MODIS) LST products. This study applied an improved gap-filling method by utilizing the discrete cosine transform-based penalized least squares (DCT-PLS) strategy in the spatial domain combined with the linear interpolation (LI) algorithm in the temporal domain. The method was applied to fill gaps in the LSWT imagery of 12 representative lakes across MP. The randomly sampled high-quality MODIS LSWT values in the spatial and temporal domains were excavated as false data gaps and considered “virtual true” validation datasets. The spatial validation results showed that the estimated LSWT for all the lake cases were comparable with the “virtual true” LSWT values, with the average values of the coefficient of determination, mean absolute error, mean square error, and root mean square error being 0.98, 0.38 °C, 0.45 °C, and 0.59 °C, respectively. Meanwhile, the error of nighttime LSWT results was relatively lower than that of daytime LSWT. For temporal interpolation validation, the LI algorithm exhibited relatively better performance and could more objectively indicate the variation in LSWT. Benefiting from the spatially and temporally well-constrained data, we analyzed the interannual and intra-annual change characteristics of the LSWTs of the 12 lakes. The long-term variations of annual and seasonal mean LSWTs in the 12 selected lakes exhibited no evident trends in 2000–2020, while presented apparent interannual fluctuations. The slight changes in the average LSWTs of the 12 selected lakes were in excellent synchronization with the surrounding LST derived from the reanalysis datasets, confirming the widely reported phenomenon of “global warming hiatus” that occurred in the early 21st century. This study improves the understanding of the LSWT variations in Mongolian lakes in response to global climate change. It has the potential to provide an effective approach for monitoring LSWT changes in other large-scale studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

On the implementation of dislocation reactions in continuum dislocation dynamics modeling of mesoscale plasticity

The continuum dislocation dynamics framework for mesoscale plasticity is intended to capture the dislocation density evolution and the deformation of crystals when subjected to mechanical loading. It does so by solving a set of transport equations for dislocations concurrently with crystal mechanics equations, with the latter being cast in the form of an eigenstrain problem. Incorporating dislocation reactions in the dislocation transport equations is essential for making such continuum dislocation dynamics predictive. A formulation is proposed to incorporate dislocation reactions in the transport equations of the vector density-based continuum dislocation dynamics. This formulation aims to rigorously enforce dislocation line continuity using the concept of virtual dislocations that close all dislocation loops involved in cross slip, annihilation, and glissile and sessile junction reactions. The addition of virtual dislocations enables us to accurately enforce the divergence free condition upon the numerical solution of the dislocation transport equations for all slip systems individually. A set of tests were performed to illustrate the accuracy of the formulation and the solution of the transport equations within the vector density-based continuum dislocation dynamics. Comparing the results from these tests with an earlier approach in which the divergence free constraint was enforced on the total dislocation density tensor or the sum of two densities when only cross slip is considered shows that the new approach yields highly accurate results. Bulk simulations were performed for a face centered cubic crystal based on the new formulation and the results were compared with discrete dislocation dynamics predictions of the same. The microstructural features obtained from continuum dislocation dynamics were also analyzed with reference to relevant experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A multidimensional study of the structure function ratio $\sigma_{LT'}$/$\sigma_0$ from hard exclusive $\pi^+$ electro-production off protons in the GPD regime

A multidimensional extraction of the structure function ratio σ LT' /σ 0 from the hard exclusive $\vec{e}$p → e'nπ+ reaction above the resonance region has been performed. The study was done based on beam-spin asymmetry measurements using a 10.6 GeV incident electron beam on a liquid-hydrogen target and the CLAS12 spectrometer at Jefferson Lab. The measurements focus on the very forward regime (t/Q 2 $\ll$ 1) with a wide kinematic range of x B in the valence regime (0.17 < x B < 0.55), and virtualities Q 2 ranging from 1.5 GeV 2 up to 6 GeV 2 . The results and their comparison to theoretical models based on Generalized Parton Distributions demonstrate the sensitivity to chiral-odd GPDs and the directly related tensor charge of the nucleon. In addition, the data is compared to an extension of a Regge formalism at high photon virtualities. It was found that the Regge model provides a better description at low Q 2 , while the GPD model is more appropriate at high Q 2 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Evidence of off-shell Higgs boson production from ZZ leptonic decay channels and constraints on its total width with the ATLAS detector

This Letter reports on a search for off-shell production of the Higgs boson using 139 fb —1 of $pp$ collision data at $\sqrt{s}$ = 13 TeV collected by the ATLAS detector at the Large Hadron Collider. The signature is a pair of $\overline{Z}$ bosons, with contributions from both the production and subsequent decay of a virtual Higgs boson and the interference of that process with other processes. The two observable final states are $\overline{Z}$ $\overline{Z}$ → 4$\ell$ and $\overline{Z}$ $\overline{Z}$ → 2$\ell$2$ν$ with $\ell$ = $e$ or $μ$. In the $\overline{Z}$ $\overline{Z}$ → 4$\ell$ final state, a dense Neural Network is used to enhance analysis sensitivity with respect to matrix element-based discrimination. The backgroundonly hypothesis is rejected with an observed (expected) significance of 3.3 (2.2) standard deviations, representing experimental evidence for off-shell Higgs boson production. Assuming that no new particles enter the production of the virtual Higgs boson, its total width can be deduced from the measurement of its off-shell production cross-section. The measured total width of the Higgs boson is ${4.5}_{—2.5}^{+3.3}$ MeV, and the observed (expected) upper limit on the total width is found to be 10.5 (10.9) MeV at 95% confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Direct-photon production in inelastic and high-multiplicity proton–proton collisions at $\sqrt{s}$ = 13 TeV

In this letter, we present the first measurement of direct photons at the transverse momentum of 1 < 𝑝 T < 6 GeV/𝑐 at midrapidity 𝜂 < 0.8 in inelastic and high-multiplicity proton–proton collisions at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV. The fraction of virtual direct photons in the inclusive virtual photon spectrum is obtained from a fit to the dielectron invariant mass spectrum. In the limit of zero invariant mass, this fraction is equal to the relative contribution of real direct photons in the inclusive real photon spectrum. Contributions from decays of light- flavour neutral mesons are estimated using independent measurements in proton–proton collisions at the same energy and the same event class. For the first time at the LHC energies, a direct-photon signal is observed at low 𝑝T in both inelastic and high-multiplicity event classes, with a significance of 3.2𝜎 and 1.9𝜎 in terms of standard deviations, correspondingly. The yield of direct photons in inelastic pp collisions is compared to perturbative QCD calculations. The integrated photon yield is studied as a function of charged-particle multiplicity and is compared to the results from other experiments and theoretical calculations. The results show a significant increase of direct- photon yield with charged-particle multiplicity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Multi-decade high-resolution regional hindcasts for wave energy resource characterization in U.S. coastal waters

Long-term, high-resolution, regional wave hindcast datasets were generated using unstructured-grid Simulating WAves Nearshore (SWAN) models for the U.S. coastal waters to support nearshore wave energy development in the U.S. including those bordering U.S. territorial islands. The model domains resolved the entire U.S. exclusive economic zones, with a spatial resolution of approximately 200 m nearshore. The regional SWAN models were driven by the global WAVEWATCH III® model outputs and run for a 42-year period from 1979 to 2020. Extensive model validations were performed using buoy observations and altimeter data. Regional resource characterization was performed based on hindcast data points at 2 km from shore and along the 100 m isobath. Aggregations of wave resource parameters were produced, and spatial and seasonal variations were analyzed for all the regions. Wave resource metrics recommended by international standards, including a 3-hour time series of six resource parameters, hourly frequency- and directionally resolved wave spectra at selected “virtual buoy” locations, and average-annual values of omni-directional wave power, significant wave height, and energy period are publicly disseminated through an Amazon Web Service and a Marine Energy Atlas web application tool to facilitate wave energy research and a wide range of coastal ocean applications.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Pose Classification Using Three-Dimensional Atomic Structure-Based Neural Networks Applied to Ion Channel–Ligand Docking

The identification of promising lead compounds showing pharmacological activities toward a biological target is essential in early stage drug discovery. With the recent increase in available small-molecule databases, virtual high-throughput screening using physics-based molecular docking has emerged as an essential tool in assisting fast and cost-efficient lead discovery and optimization. However, the best scored docking poses are often suboptimal, resulting in incorrect screening and chemical property calculation. We address the pose classification problem by leveraging data-driven machine learning approaches to identify correct docking poses from AutoDock Vina and Glide screens. To enable effective classification of docking poses, we present two convolutional neural network approaches: a three-dimensional convolutional neural network (3D-CNN) and an attention-based point cloud network (PCN) trained on the PDBbind refined set. We demonstrate the effectiveness of our proposed classifiers on multiple evaluation data sets including the standard PDBbind CASF-2016 benchmark data set and various compound libraries with structurally different protein targets including an ion channel data set extracted from Protein Data Bank (PDB) and an in-house KCa3.1 inhibitor data set. Our experiments show that excluding false positive docking poses using the proposed classifiers improves virtual high-throughput screening to identify novel molecules against each target protein compared to the initial screen based on the docking scores.

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

Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. More nefariously, intensity borrowing by the low-energy charge-transfer states can affect intensities of the valence excitations even if those excitation energies are accurate. At higher excitation energies, it is difficult to distinguish spurious CT states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization scheme that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces the concept of “natural charge-transfer orbitals”, which provide a means to isolate characteristic pairs of orbitals that are most likely to participate in a CTTS excitation. The projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation. We apply this method to analyze an ab initio molecular dynamics (MD) trajectory of I-(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with experimental measurements for bulk I-(aq), and only one-third of the water molecules in the I-(H2O)96 simulation cell need to be described with LR-TDDFT in order to obtain excitation energies that are converged to < 0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments

Carter-Fenk, Kevin D.↗
DOE OSTI2020

Phase Behavior of Diblock Copolymer–Homopolymer Ternary Blends with a Compositionally Asymmetric Diblock Copolymer

The phase behavior of ternary polymer blends comprising poly(cyclohexylethylene) (C) and polyethylene (E) homopolymers and a compositionally asymmetric CE diblock copolymer with f C = 0.67 was investigated, where f C is the volume fraction of C. Morphology was established in the phase prism (volume fractions of C, E, and CE vs. temperature) by optical transmission, small-angle X-ray scattering (SAXS), and small-angle neutron scattering (SANS) measurements. The locations of lamellar (LAM), hexagonally packed cylinders (HEX), and gyroid (GYR) ordered phases are shifted significantly towards lower fractions of C homopolymer compared to previous results obtained from ternary polymer blends with a symmetric diblock copolymer (f C = 0.5). Conversely, the Scott line of critical points, which delineates the boundary between single-phase disorder and macroscopic phase separation, remains virtually unchanged, coincident with the fraction of C homopolymer associated with the binary homopolymer blend critical composition. A central finding of this study is that the line of nearly congruent order-disorder transitions, where the LAM phase melts virtually directly into the disordered state, is decoupled in composition from that of the Scott line of critical points. A wide range of phase space between the compositions associated with the congruent transition and Scott line was identified as containing a microemulsion morphology. This study demonstrates that diblock copolymer compositional asymmetry significantly impacts the ordered phase regime but has a marginal effect on the region displaying macroscopic phase separation. It also provides useful guidance for tuning interfacial curvature, a crucial factor in the formation of bicontinuous microemulsions (BµEs).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Advancing Cross-Disciplinary Understanding of Land-Atmosphere Interactions

The evolution of disciplinary silos and increasingly narrow disciplinary boundaries have together resulted in one-sided approaches to the study of land-atmosphere interactions—a field that requires a bi-directional approach to understand the complex feedbacks and interactions that occur. The integration of surface flux and atmospheric boundary layer measurements is therefore essential to advancing our understanding. The Land-Atmosphere 2021 workshop (held virtually, June 10-11, 2021) involved almost 300 participants from around the world and promoted cross-discipline collaboration by way of talks from invited speakers, moderated discussions, breakout sessions, and a virtual poster session. The workshop focused on five main theme areas: “big picture” overview, instrumentation and remote sensing, modeling, water, and aerosols and clouds. In talks and breakout groups, there were frequent calls for more AmeriFlux sites to be instrumented for boundary layer height measurements, and for the development of some “super sites” where profiling instruments would be deployed. There was further agreement on the need for the standardization of various datasets. There was also a consensus that funding agencies need to be willing to support the sorts of large projects (including associated instrumentation) which can drive interdisciplinary work. Early-career scientists, in particular, expressed enthusiasm for working across disciplinary boundaries but noted that there need to be more financial support and training opportunities so they would be better prepared for interdisciplinary work. Investment in these career development opportunities would enable today's cohort of early-career scientists to advance the frontiers of interdisciplinary work over the next couple of decades.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Snapshot multispectral imaging using a diffractive optical network

Abstract Multispectral imaging has been used for numerous applications in e.g., environmental monitoring, aerospace, defense, and biomedicine. Here, we present a diffractive optical network-based multispectral imaging system trained using deep learning to create a virtual spectral filter array at the output image field-of-view. This diffractive multispectral imager performs spatially-coherent imaging over a large spectrum, and at the same time, routes a pre-determined set of spectral channels onto an array of pixels at the output plane, converting a monochrome focal-plane array or image sensor into a multispectral imaging device without any spectral filters or image recovery algorithms. Furthermore, the spectral responsivity of this diffractive multispectral imager is not sensitive to input polarization states. Through numerical simulations, we present different diffractive network designs that achieve snapshot multispectral imaging with 4, 9 and 16 unique spectral bands within the visible spectrum, based on passive spatially-structured diffractive surfaces, with a compact design that axially spans ~72 λ m , where λ m is the mean wavelength of the spectral band of interest. Moreover, we experimentally demonstrate a diffractive multispectral imager based on a 3D-printed diffractive network that creates at its output image plane a spatially repeating virtual spectral filter array with 2 × 2 = 4 unique bands at terahertz spectrum. Due to their compact form factor and computation-free, power-efficient and polarization-insensitive forward operation, diffractive multispectral imagers can be transformative for various imaging and sensing applications and be used at different parts of the electromagnetic spectrum where high-density and wide-area multispectral pixel arrays are not widely available.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Prediction on X-ray output of free electron laser based on artificial neural networks

Abstract Knowledge of x-ray free electron lasers’ (XFELs) pulse characteristics delivered to a sample is crucial for ensuring high-quality x-rays for scientific experiments. XFELs’ self-amplified spontaneous emission process causes spatial and spectral variations in x-ray pulses entering a sample, which leads to measurement uncertainties for experiments relying on multiple XFEL pulses. Accurate in-situ measurements of x-ray wavefront and energy spectrum incident upon a sample poses challenges. Here we address this by developing a virtual diagnostics framework using an artificial neural network (ANN) to predict x-ray photon beam properties from electron beam properties. We recorded XFEL electron parameters while adjusting the accelerator’s configurations and measured the resulting x-ray wavefront and energy spectrum shot-to-shot. Training the ANN with this data enables effective prediction of single-shot or average x-ray beam output based on XFEL undulator and electron parameters. This demonstrates the potential of utilizing ANNs for virtual diagnostics linking XFEL electron and photon beam properties.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Stochastic 3D reconstruction of cracked polycrystalline NMC particles using 2D SEM data

Li-ion battery performance is strongly influenced by the 3D microstructure of its cathode particles. Cracks within these particles develop during calendaring and cycling, reducing connectivity but increasing reactive surface, making their impact on battery performance complex. Understanding these contradictory effects requires a quantitative link between particle morphology and battery performance. However, informative 3D imaging techniques are time-consuming, costly and rarely available, such that analyses often have to rely on 2D image data. This paper presents a novel stereological approach for generating virtual 3D cathode particles exhibiting crack networks that are statistically equivalent to those observed in 2D sections of experimentally measured particles. Consequently, 2D image data suffices for deriving a full 3D characterization of cracked cathodes particles. Such virtually generated 3D particles could serve as geometry input for spatially resolved electro-chemo-mechanical simulations to enhance our understanding of structure-property relationships of cathodes in Li-ion batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Vertical full-colour micro-LEDs via 2D materials-based layer transfer

Micro-LEDs (µLEDs) have been explored for augmented and virtual reality display applications that require extremely high pixels per inch and luminance. However, conventional manufacturing processes based on the lateral assembly of red, green and blue (RGB) µLEDs have limitations in enhancing pixel density. Recent demonstrations of vertical µLED displays have attempted to address this issue by stacking freestanding RGB LED membranes and fabricating top-down, but minimization of the lateral dimensions of stacked µLEDs has been difficult. Here we report full-colour, vertically stacked µLEDs that achieve, to our knowledge, the highest array density (5,100 pixels per inch) and the smallest size (4 µm) reported to date. This is enabled by a two-dimensional materials-based layer transfer technique that allows the growth of RGB LEDs of near-submicron thickness on two-dimensional material-coated substrates via remote or van der Waals epitaxy, mechanical release and stacking of LEDs, followed by top-down fabrication. The smallest-ever stack height of around 9 µm is the key enabler for record high µLED array density. We also demonstrate vertical integration of blue µLEDs with silicon membrane transistors for active matrix operation. Furthermore, these results establish routes to creating full-colour µLED displays for augmented and virtual reality, while also offering a generalizable platform for broader classes of three-dimensional integrated devices.

42 ENGINEERING↗
DOE OSTI2026

An agentic artificially intelligent X-ray scientist

Executing experimental tasks in both normal research laboratories and large-scale scientific facilities often requires extensive human supervision and remains a key challenge on the path to fully autonomous, artificial intelligence (AI)-driven science. Here we demonstrate a large language model-driven agent that autonomously performs X-ray sample alignment on a synchrotron beamline by planning actions, executing instrumental commands, interpreting observations and iterating towards experimental goals. Based on existing large language models with structured tool-use via the model context protocol, our AI X-ray scientist was guided and tested using an in-house-built virtual experimental setup that mirrors a six-circle diffractometer at an operational synchrotron beamline. The agentic workflow developed in the virtual environment was directly deployed on a real beamline, where it correctly identified reference reflections and determined the orientation matrix, an essential first step in any type of single-crystal scattering experiment. Our AI X-ray scientist responded effectively to unexpected experimental conditions, demonstrating adaptive problem-solving and readiness for addressing practical experimental situations. Our study provides a step towards autonomous operation across diverse experimental environments at large-scale scattering facilities.

Chen, Zhantao (ORCID:0000000319543868)↗
DOE OSTI2023

Advancing energy storage through solubility prediction: leveraging the potential of deep learning

Solubility prediction plays a crucial role in energy storage applications, such as redox flow batteries, because it directly affects the efficiency and reliability. Researchers have developed various methods that utilize quantum calculations and descriptors to predict the aqueous solubilities of organic molecules. Notably, machine learning models based on descriptors have shown promise for solubility prediction. As deep learning tools, graph neural networks (GNNs) have emerged to capture complex structure–property relationships for material property prediction. Specifically, MolGAT, a type of GNN model, was designed to incorporate n-dimensional edge attributes, enabling the modeling of intricacies in molecular graphs and enhancing the prediction capabilities. In a previous study, MolGAT successfully screened 23 467 promising redox-active molecules from a database of over 500 000 compounds, based on redox potential predictions. This study focused on applying the MolGAT model to predict the aqueous solubility (log S) of a broad range of organic compounds, including those previously screened for redox activity. The model was trained on a diverse sample of 8494 organic molecules from AqSolDB and benchmarked against literature data, demonstrating superior accuracy compared with other state of the art graph-based and descriptor-based models. Subsequently, the trained MolGAT model was employed to screen redox-active organic compounds identified in the first phase of high-throughput virtual screening, targeting favorable solubility in energy storage applications. The second round of screening, which considered solubility, yielded 12 332 promising redox-active and soluble organic molecules suitable for use in aqueous redox flow batteries. Thus, the two-phase high-throughput virtual screening approach utilizing MolGAT, specifically trained for redox potential and solubility, is an effective strategy for selecting suitable intrinsically soluble redox-active molecules from extensive databases, potentially advancing energy storage through reliable material development. This indicates that the model is reliable for predicting the solubility of various molecules and provides valuable insights for energy storage, pharmaceutical, environmental, and chemical applications.

25 ENERGY STORAGE↗
DOE OSTI2022

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗