SEARCH · Engineering Papers
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.
We discuss the current progress in applications of machine learning (ML) and artificial intelligence (AI) to meet the challenges of computational drug discovery. We identify several areas where existing methods have the potential to accelerate pharmaceutical research and disrupt more traditional approaches.
Forming a hetero-interface is a materials-design strategy that can access an astronomically large phase space. However, the immense phase space necessitates a high-throughput approach for an optimal interface design. Here we introduce a high-throughput computational framework, InterMatch, for efficiently predicting charge transfer, strain, and superlattice structure of an interface by leveraging the databases of individual bulk materials. Specifically, the algorithm reads in the lattice vectors, density of states, and the stiffness tensors for each material in their isolated form from the Materials Project. From these bulk properties, InterMatch estimates the interfacial properties. We benchmark InterMatch predictions for the charge transfer against experimental measurements and supercell density-functional theory calculations. We then use InterMatch to predict promising interface candidates for doping transition metal dichalcogenide MoSe 2 . Finally, we explain experimental observation of factor of 10 variation in the supercell periodicity within a few microns in graphene/α-RuCl 3 by exploring low energy superlattice structures as a function of twist angle using InterMatch. We anticipate our open-source InterMatch algorithm accelerating and guiding ever-growing interfacial design efforts. Moreover, the interface database resulting from the InterMatch searches presented in this paper can be readily accessed online.
In the 25 years since LLNL broke ground on the National Ignition Facility (NIF) - a cornerstone of the National Nuclear Security Administration’s Stockpile Stewardship Program - the NIF and Photon Science (NIF&PS) Directorate has steadily pushed the boundaries of laser physics, nonlinear optics, and photonics in service of inertial confinement fusion experiments, and advanced photon sources development. Beyond the 192 beamlines of the primary laser facility with its world-record energy and fusion yields, NIF&PS capabilities have expanded to include chirped-pulse amplification, kilojoule petawatt-class short-pulse systems to generate hard x-ray radiographic probes, high-average power lasers as secondary sources for generation or wake field particle acceleration along directed energy applications, and other scientific pursuits. The Computing Directorate partners with NIF&PS on numerous projects and technologies from diagnostic measurements and high performance control systems to data analysis, information technology infrastructure, and scientific simulation code development. “This strong collaboration enables cutting-edge science, particularly with regard to physics-based modeling and simulation,” says physicist Jean-Michel Di Nicola. One crucial capability for NIF&PS is the Virtual Beamline (VBL) laser simulation code. It can model all the major laser physics and technology involved in the design optimization, commissioning, and operations of advanced laser architectures from tabletop to NIF-scale. VBL provides researchers with high-fidelity models and high-resolution calculations of laser performance predictions—including for the entire NIF laser system, the Advanced Radiographic Capability, parts of the High Repetition-Rate Advanced Petawatt Laser System laser delivered to the Czech Republic, and the Optical Science Laser. After more than two decades of experimentally verified physics and computing enhancements, this unique workhorse code is wrapping up another major milestone: migration from Java to C++ with a flurry of user interface and optimization features, as well as laser physics and high-resolution enhancements thanks to parallel execution on Livermore Computing platforms.
Exciton-polaritons (EPs) are an emerging approach to achieve strong light-matter coupling and robust quantum phenomena, with application to areas such as photonic neuromorphic computing for instance reservoir computing. 2D transition metal dichalcogenides (TMDs) are promising candidates to host EPs due to their potential for strong coupling as evidenced by large Rabi splitting energies, which in the case of few-layer MoS2 has been demonstrated to reach 293 meV when coupled to the C exciton, placing it in the ultrastrong coupling regime, at room temperature. Here we examine how the Rabi splitting of CVD-grown monolayer MoS2 in a cavity can be optimized by tuning two variables: cavity geometry (thickness of silver and dielectric spacers) and the targeted exciton for coupling (A at about 1.9 eV, B at about 2.06 eV, or C at about 2.9 eV). With the cavity geometry, there is a tradeoff between optical field confinement and transmission. With the targeted exciton, there may be a tradeoff between Rabi splitting and lifetime. Transfer matrix simulations are used to inform how thick the cavity layers need to be, and at what angle measurements should be performed, to achieve Rabi splitting at a certain exciton energy. Experimentally, the Rabi splitting is measured by placing the cavity at a specified angle and measuring the optical reflectance and transmittance. By fitting the data, we find Rabi splittings of 124 meV, 194 meV, and 288 meV for the A, B, and C excitons, respectively, which places them in the strong, and on the edge of the ultrastrong, coupling regime.
Simulating the Hubbard model is of great interest to a wide range of applications within condensed matter physics, however its solution on classical computers remains challenging in dimensions larger than one. The relative simplicity of this model, embodied by the sparseness of the Hamiltonian matrix, allows for its efficient implementation on quantum computers, and for its approximate solution using variational algorithms such as the variational quantum eigensolver. While these algorithms have been shown to reproduce the qualitative features of the Hubbard model, their quantitative accuracy in terms of producing true ground state energies and other properties, and the dependence of this accuracy on the system size and interaction strength, the choice of variational ansatz, and the degree of spatial inhomogeneity in the model, remains unknown. Here we present a rigorous classical benchmarking study, demonstrating the potential impact of these factors on the accuracy of the variational solution of the Hubbard model on quantum hardware, for systems with up to 32 qubits. We find that even when using the most accurate wavefunction ansätze for the Hubbard model, the error in its ground state energy and wavefunction plateaus for larger lattices, while stronger electronic correlations magnify this issue. Concurrently, spatially inhomogeneous parameters and the presence of off-site Coulomb interactions only have a small effect on the accuracy of the computed ground state energies. Our study highlights the capabilities and limitations of current approaches for solving the Hubbard model on quantum hardware, and we discuss potential future avenues of research.
Time-accurate Navier Stokes computations were performed to study self-excited oscillatory shear layer behavior over a two-dimensional rectangular cavity and to gain insight into the vortex-acoustic coupling that sustains these oscillations. Results were obtained solving the full, unsteady compressible Navier-Stokes equations for Reynolds-averaged turbulent air flow past a cavity with L/d = 2. Conditions duplicated test conditions in experiments by Rossiter at Mach numbers varied from 0.6 to 1.2. The computational scheme was the USA 2D series code developed by Rockwell's Science Center with an algebraic turbulence model for the shear layer and recirculating region. Two physical conditions were simulated: an isolated cavity situated (1) in an infinite plane wall and (2) in an infinite channel with an opposing parallel wall. Computed results for (1) showed excellent agreement with Rossiter's measured amplitudes and frequencies. A tuned acoustic duct resonance effect was observed in (2), demonstrating a controlled change in the acoustic feedback mechanism and resultant amplification in computed oscillations.
Inverse problems play a central role in experimentation and theory/data comparisons for many areas of modern Nuclear Physics (NP) and High-Energy Physics (HEP). Bayes’s Theorem is a powerful tool for solving Inverse Problems, providing conceptually transparent and unbiased constraints on theoretical parameters and their uncertainties (“Bayesian Inference”) and enabling the quantification of agreement or tension between models and data. However, analyses based on Bayesian Inference are often challenging for NP and HEP applications, either because of the large number of parameters in the problem, the high computational cost, or both. We propose a multi-institutional collaboration to develop and deploy novel Bayesian analysis tools that advance the scientific scope of a broad range of current and future NP experiments. This project brings together NP domain scientists working on several high-profile NP projects for which new, high-performance Bayesian Uncertainty Quantification (“Bayesian UQ”) methods are essential to carry out the science, and data scientists who are developing state-of-the-art methods applicable to these problems. The NP projects in this proposal comprise measurements of the mass and fundamental nature of the neutrino; study of the Quark-Gluon Plasma that filled the early universe; and mapping of natural and anthropogenic radiation environments. While these NP projects have very different scientific goals, with datasets and analysis approaches that differ significantly, they share common requirements for improving computationally intensive Bayesian analyses using advanced Machine Learning algorithms and will benefit strongly from a coherent effort to develop general solutions. This proposal brings together these projects and forefront ML-based data science algorithms to develop such general solutions. The methods developed in this project will also be more widely applicable, thereby advancing science in the larger Nuclear Physics portfolio.
Topological superconductor (TSC) hosting Majorana fermions has been established as a milestone that may shift our scientific trajectory from research to applications in topological quantum computing. Recently, superconducting Pd-Bi binaries have attracted great attention as a possible medium for the TSC phase as a result of their large spin-orbit coupling strength. Here, we report a systematic high-resolution angle-resolved photoemission spectroscopy (ARPES) study on the normal state electronic structure of superconducting α - PdBi 2 (T c = 1.7 K). Our results show the presence of Dirac states at higher-binding energy with the location of the Dirac point at 1.26 eV below the chemical potential at the zone center. Furthermore, the ARPES data indicate multiple band crossings at the chemical potential, consistent with the metallic behavior of α - PdBi 2 . Our detailed experimental studies are complemented by first-principles calculations, which reveal the presence of surface Rashba states residing in the vicinity of the chemical potential. The obtained results provide an opportunity to investigate the relationship between superconductivity, topology, and the Majorana fermion, as well as explore pathways to possible future platforms for topological quantum computing.
We describe a highly optimized implementation of MPI domain decomposition in a GPU-enabled, general-purpose molecular dynamics code, HOOMD-blue (Anderson and Glotzer, 2013). Here, our approach is inspired by a traditional CPU-based code, LAMMPS (Plimpton, 1995), but is implemented within a code that was designed for execution on GPUs from the start (Anderson et al., 2008). The software supports short-ranged pair force and bond force fields and achieves optimal GPU performance using an autotuning algorithm. We are able to demonstrate equivalent or superior scaling on up to 3375 GPUs in Lennard-Jones and dissipative particle dynamics (DPD) simulations of up to 108 million particles. GPUDirect RDMA capabilities in recent GPU generations provide better performance in full double precision calculations. For a representative polymer physics application, HOOMD-blue 1.0 provides an effective GPU vs. CPU node speed-up of 12.5x.
Monte Carlo (MC) integration is an important calculational technique in the physical sciences. Practical considerations require that the calculations are performed as accurately as possible for a given set of computational resources. To improve the accuracy of MC integration, a number of useful variance reduction algorithms have been developed, including importance sampling and control variates. In this work, we demonstrate how these two methods can be applied simultaneously, thus combining their benefits. We provide a python wrapper, named CoVVVR, which implements our approach in the VEGAS program. The improvements are quantified with several benchmark examples from the literature.
This discussion will include a short review of the challenges that must be overcome if computational physics technology is to have a larger impact on the design cycles of U.S. aerospace companies. Some of the potential solutions to these challenges may come from the information sciences fields. A few examples of potential computational physics/information sciences synergy will be presented, as motivation and inspiration for the Improving The CFD Applications Process Workshop.
Abstract Objective We describe a detailed solution for maintaining high-capacity, data-intensive network flows (eg, 10, 40, 100 Gbps+) in a scientific, medical context while still adhering to security and privacy laws and regulations. Materials and Methods High-end networking, packet-filter firewalls, network intrusion-detection systems. Results We describe a “Medical Science DMZ” concept as an option for secure, high-volume transport of large, sensitive datasets between research institutions over national research networks, and give 3 detailed descriptions of implemented Medical Science DMZs. Discussion The exponentially increasing amounts of “omics” data, high-quality imaging, and other rapidly growing clinical datasets have resulted in the rise of biomedical research “Big Data.” The storage, analysis, and network resources required to process these data and integrate them into patient diagnoses and treatments have grown to scales that strain the capabilities of academic health centers. Some data are not generated locally and cannot be sustained locally, and shared data repositories such as those provided by the National Library of Medicine, the National Cancer Institute, and international partners such as the European Bioinformatics Institute are rapidly growing. The ability to store and compute using these data must therefore be addressed by a combination of local, national, and industry resources that exchange large datasets. Maintaining data-intensive flows that comply with the Health Insurance Portability and Accountability Act (HIPAA) and other regulations presents a new challenge for biomedical research. We describe a strategy that marries performance and security by borrowing from and redefining the concept of a Science DMZ, a framework that is used in physical sciences and engineering research to manage high-capacity data flows. Conclusion By implementing a Medical Science DMZ architecture, biomedical researchers can leverage the scale provided by high-performance computer and cloud storage facilities and national high-speed research networks while preserving privacy and meeting regulatory requirements.
Compared to oxides, the nitrides are relatively unexplored, making them a promising chemical space for novel materials discovery. Of particular interest are nitrogen-rich nitrides, which often possess useful semiconducting properties for electronic and optoelectronic applications. However, such nitrogen-rich compounds are generally metastable, and the lack of a guiding theory for their synthesis has limited their exploration. Here, we review the remarkable metastability of observed nitrides, and examine the thermodynamics of how reactive nitrogen precursors can stabilize metastable nitrogen-rich compositions during materials synthesis. We map these thermodynamic strategies onto a predictive computational search, training a data-mined ionic substitution algorithm specifically for nitride discovery, which we combine with grand-canonical DFT-SCAN phase stability calculations to compute stabilizing nitrogen chemical potentials. We identify several new nitrogen-rich binary nitrides for experimental investigation, notably the transition metal nitrides Mn3N4, Cr3N4, V3N4, and Nb3N5, the main group nitride SbN, and the pernitrides FeN2, CrN2, and Cu2N2. By formulating rational thermodynamic routes to metastable compounds, we expand the search space for functional technological materials beyond equilibrium phases and compositions.
Optoelectronic devices based on hybrid perovskites have demonstrated outstanding performance within a few years of intense study. However, commercialization of these devices requires barriers to their development to be overcome, such as their chemical instability under operating conditions. In order to investigate this instability and its consequences, the electric field applied to single crystals of methylammonium lead bromide (CH 3 NH 3 PbBr 3 ) is varied, and changes are mapped in both their elemental composition and photo-luminescence. Synchrotron-based nanoprobe X-ray fluorescence (nano-XRF) with 250 nm resolution reveals quasi-reversible field-assisted halide migration, with corresponding changes in photoluminescence. It is observed that higher local bromide concentration is correlated to superior optoelectronic performance in CH 3 NH 3 PbBr 3 . A lower limit on the electromigration rate is calculated from these experiments and the motion is interpreted as vacancy-mediated migration based on nudged elastic band density functional theory (DFT) simulations. The XRF mapping data provide direct evidence of field-assisted ionic migration in a model hybrid-perovskite thin single crystal, while the link with photoluminescence proves that the halide stoichiometry plays a key role in the optoelectronic properties of the perovskite.
The longest-wavelength π-to-π* electronic excitations of rhodamine-like dyes (RDs) with different group16 heteroatoms (O, S, Se, Te) have been investigated. Time-dependent Kohn–Sham theory (TDKST) calculations were compared with coupled-cluster (CC) and equations-of-motion (EOM) CC results for π-to-π* singlet and triplet excitations. The RDs exhibit characteristics in the TDKST calculations that are very similar to previously investigated cyanine dyes, in the sense that the singlet energies obtained with nonhybrid functionals are too high compared with the CC results at the SD(T) level. The errors became increasingly larger for functionals with increasing amounts of exact exchange. TDKST with all tested functionals led to severe underestimations of the corresponding triplet excitations and overestimations of the singlet--triplet gaps. Long-range-corrected range-separated exchange and "optimal tuning" of the range separation parameter did not significantly improve the TDKST results. A detailed analysis suggests that the problem is differential electron correlation between the ground and excited states, which is not treated sufficiently by the relatively small integrals over the exchange-correlation response kernel that enters the excitation energy expression. As a result, numerical criteria are suggested that may help identify "cyanine-like" problems in TDKST calculations of excitation spectra.
Reduced-order structure–property (S-P) linkages play a pivotal role in the tailored design of materials for advanced engineering components. There is a critical need to distill these from the simulation datasets aggregated using sophisticated, computationally expensive, physics-based numerical tools (e.g., finite element methods). The recent emergence of materials data science approaches has opened new avenues for addressing this challenge. In this paper, we critically compare the relative merits of the application of four distinct machine learning approaches for their efficacy in extracting microstructure-property linkages from the finite element simulation data aggregated on high-contrast elastic composites with different microstructures. The machine learning approaches selected for the study have included different combinations of local/global and parametric/nonparametric approaches. Furthermore, the nonparametric approaches selected for this study are based on Gaussian Process (GP) models that allow for a formal treatment of uncertainty quantification in the predicted values. The predictive performances of these different approaches have been compared against each other using rigorous cross-validation error metrics. Furthermore, their sensitivity to both the dataset size and dimensionality has been investigated.
We present an adaptive sampling method supplemented by a distributed database and a prediction method for multiscale simulations using the Heterogeneous Multiscale Method. A finite-volume scheme integrates the macro-scale conservation laws for elastodynamics, which are closed by momentum and energy fluxes evaluated at the micro-scale. In the original approach, molecular dynamics (MD) simulations are launched for every macro-scale volume element. Our adaptive sampling scheme replaces a large fraction of costly micro-scale MD simulations with fast table lookup and prediction. The cloud database Redis provides the plain table lookup, and with locality aware hashing we gather input data for our prediction scheme. For the latter we use kriging, which estimates an unknown value and its uncertainty (error) at a specific location in parameter space by using weighted averages of the neighboring points. We find that our adaptive scheme significantly improves simulation performance by a factor of 2.5 to 25, while retaining high accuracy for various choices of the algorithm parameters.