Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Component sizing”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Scalable multiscale modeling of platelets with 100 million particles

Here, we developed the core components of the AI-aided multiple time stepping algorithm for multiscale modeling of cell dynamics. This algorithm was implemented and analyzed on two supercomputer architectures with an application of simulating the aggregation of 250 platelets, or 102 million particles. To scale on these computers with complex memory and network architectures with GPUs, we devised a biomechanics-informed task mapping scheme to optimize load imbalance, communications, and memory utilization. Our simulations, scaling well up to 192 nodes on a Summit-like supercomputer with a peak speed of 11 petaflops, achieved a rate of 423 μs/day which is 500 times faster than the conventional algorithm using static time step and this has enabled studies of record size blood clots at record spatial–temporal resolutions. Additionally, we discovered the sensitive dependence of the scalability and execution time on the methods of decomposition, CPU–GPU coupling, and task mapping.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Processing and mechanical characterization of short carbon fiber-reinforced epoxy composites for material extrusion additive manufacturing

Fiber-reinforced polymer composites have been extensively utilized in recent years as feedstock materials for material extrusion additive manufacturing (AM) processes to improve strength, stiffness, and functionality of printed parts over unfilled printed polymers. However, the widespread adoption of AM of fiber-reinforced polymer composites requires a deeper understanding of the process-structure-property relationships in printed components, and such relationships are not well understood yet. Fiber length is critically important to the mechanical performance of short fiber composites, but very few studies to-date have focused on how the fiber length distribution (FLD) evolves during processing of composite feedstocks and how this evolution affects printing behavior and mechanical properties in 3D-printed composites. Here, FLD is measured for carbon fiber reinforced epoxy composites over a wide range of ink compositions and shear mixing times, and the distributions are fit with a Weibull-type distribution function. The effects of FLD on the tradeoff between ink processability, ink rheology, printing behavior and mechanical properties are investigated. Furthermore, the effects of printing parameters (nozzle size and print speed) on mechanical anisotropy and fiber orientation distribution (FOD) in printed composites are explored. Mechanical properties of printed composites are characterized via 3 pt-flexural testing, and microstructure is investigated using optical and scanning electron microscopy (SEM), and x-ray computed tomography. Finally, the fitted Weibull parameters are fed into a composite model that incorporates FLD and FOD, and model predictions are found to be in excellent agreement with experimental observations.

3D printing↗
DOE OSTI2021

Investigation of neutron noise in a micro-scale, natural circulation molten salt fission battery system

The Molten Salt Reactor (MSR) concept is a rapidly evolving Generation IV design that has recently attracted favorable attention due to the potential for reducing waste generation, realizing passive safety features, and seizing on the opportunity for cost effective economics. A specific novel micro-MSR, natural circulation, battery design concept involves placing all primary components within a single reactor vessel containment without the need of primary forced pumping. This small modular, integral design presents potential cost savings while producing safety, reliable, and transportable carbon-free power for decades. This investigation evaluates the neutron noise induced by density and flow fluctuations in a natural circulation MSR battery concept being developed at the University of Idaho. This study finds numerical solutions to the one-dimensional, one-group, coupled diffusion equations for a bare, homogeneous core, and evaluates neutron flux and delayed neutron precursor concentration noise due to core flow and fuel salt density perturbation sources. Noise analysis shows both point kinetic and space-dependent behavior is present despite the small size, low-flow, closely-coupled nature of the natural circulation MSR concept. Analysis shows that the low flow of the convective system results in noise behavior typically representative of a larger forced-circulation MSR system. Finally, noise techniques may also be useful in MSR core diagnostics as a non-invasive, low-cost radiological option to traditional monitoring schemes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Passive deconfinement of runaway electrons using an in-vessel helical coil

A helical coil designed to passively generate non-axisymmetric fields during a plasma disruption is shown (via electromagnetic analysis, linear MHD modeling, and relativistic drift orbit tracing) to be effective at deconfining runaway electrons (REs) on a time scale significantly faster than the plasma current quench. Magnetic equilibria from DIII-D RE-producing scenarios are used to calculate the toroidal electric field generated during the current quench phase of a disruption, which in turn drives current in the proposed n = 1 in-vessel helical coil, without the need for any external power supplies or disruption detection or prediction techniques. Simulations of the plasma evolution using the TokSys GS Evolve code predict the inductive coupling of coil currents up to 12% of the pre-disruption plasma current into the helical coil. The coil geometry is parametrically varied to maximize both the non-resonant and resonant components of the 3D magnetic perturbation, resulting in δB/B ≈ 10 –2 and a vacuumisland overlapwidth of up to 0.7ψ N . The REORBIT module of the MARS-F code is used to model the full non-axisymmetric magnetic field and trace RE drift orbits to determine the effect on RE deconfinement, with up to 70% of the RE orbits lost after 0.2 ms. A two-stage evolution of the RE orbit loss fraction is observed to be caused by resonant trapping between multiple magnetic island chains. Finally, electromagnetic and thermal stresses on the coil are calculated to be within operational limits for installation in DIII-D, and scale favorably to a reactor-size device. Furthermore, these findings motivate future experimental study of the helical coil concept in DIII-D or other tokamaks.

3D coil↗
DOE OSTI2024

Performance prediction applying different reduced turbulence models to the SMART tokamak

The SMall Aspect Ratio Tokamak (SMART) is currently being commissioned at the University of Seville and will be able to compare the performance of positive and negative triangularity plasmas at low aspect ratio. Predictive simulations have been performed for different machine scenarios and heating schemes using the TRANSP code. The objectives of these simulations are to predict the parameters expected in positive triangularity plasmas, to guide diagnostic development, and to validate transport models. Several reduced turbulence models have been used to predict electron and ion temperatures for the operational phase 2. All models provide similar results from approximately mid-radius to the separatrix but important discrepancies are found in the core region. These positive triangularity results are compared with experiments from a similar size machine like GLOBUS-M2. The multi-mode model (MMM) shows the best agreement. Simulations with different boundary conditions have been performed and no strong differences have been observed between them. The impact of neutral beam injection (NBI) on the predicted profiles has also been addressed. Rotation reduces turbulence levels so higher temperatures are achieved when included in the simulations. Studying the different contributions to the thermal diffusivities, it is observed that electron temperature gradient (ETG) turbulence dominates at the plasma core while micro-tearing modes (MTM) dominate at the edge in the electron channel. In the ion channel, the neoclassical contribution is dominant at the core and at the very edge while the Weiland component, which includes ion temperature gradient mode (ITG), trapped electron mode (TEM), kinetic ballooning mode (KBM), peeling mode (PM) and collisionless and collision dominated magnetohydrodynamic (MHD) modes governs the mid-radius region. For phase 3, two plasmas with different electron densities have been studied. The case with lower density matches well a specific discharge of GLOBUS-M2. The higher density plasma shows high performance with β N ≈ 3.8.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Conceptual Design of a Tension Leg Platform With 22.3 MW Vertical Axis Turbine

Here, this paper presents the conceptual design of a tension leg platform (TLP) for the ARCUS “towerless” vertical-axis wind turbine (VAWT). VAWTs are ideal for floating offshore sites and have several advantages over horizontal-axis wind turbines (HAWT) including reduced top mass, lower center of gravity, increased energy capture, and in turn lower cost. The towerless ARCUS VAWT drives these advantages further through increased structural efficiency and by enabling more optimized TLP designs with simplified installation procedures. For hull sizing, we have studied three turbine sizes with corresponding power ratings of 5.1 MW, 10.4 MW and 22.3 MW. The largest turbine was identified as having the greatest potential to reduce the levelized cost of energy (LCOE) and is the reference size used for the further detailed design process. The conceptual design of the VAWT TLP has been awarded with an ABS Approval in Principle Certificate. This paper contains brief analysis results and design findings for a TLP designed to house a VAWT, including the following topics: • Applicable Design Codes • Metocean Conditions • ARCUS Turbine Loads • Design Load Cases and Requirements - Pre-service TLP Stability - In-place TLP Global Performance • Platform Configurations, Hull Structure Scantling Design, Weight and CG Estimation, and General Arrangement Drawings • Hull Ballast Plan for both Pre-service and In-place Conditions • Pre-service Quayside Integration, Transportation and Wet Tow Stability Analysis • Global Performance Analysis for Motions and Tendon tensions • Summary of cost components and system levelized cost of energy

17 WIND ENERGY↗
DOE OSTI2022

Solid Phase Processing of Lead-Free Brass with Carbon Additives

Currently, “lead-free” brass alloys (like C27450/C27451/C6930), used extensively in drinking water fixtures and automotive, electrical, and electronic applications contain maximum 0.25% lead to maintain mechanical performance and machinability. Adding graphite to brass as an alternative to lead, using casting, powder metallurgy, and extrusion methods, has been explored previously. However, all these methods have proven to be energy-, time-, and resource-intensive, while not enabling performance equivalent to that of C36000 brass. In this project, we developed a one-step approach using friction extrusion and ShAPE to make lead-free brass/graphite components such as wires, rods and tubes with mechanical performance equivalent to commercial lead-free brass alloys. Manufacturing temperatures were maintained ~550-730 °C with feed rates ranging between 4 – 25 mm/min. Results show larger grains at the center of the rods and wires with smaller grains developing at the edges. Graphite particles were sheared in the direction of extrusion with higher strains observed towards the edges. Hardness of the brass/graphite samples was over 25% higher than that of the corresponding brass-only samples (rods and wires), also friction extruded. Our results show that the sub-micron graphite plays an important role in limiting process temperature and restraining grain growth during friction extrusion, thus reducing grain size in composites.

36 MATERIALS SCIENCE↗
DOE OSTI2024

An Operations and Maintenance Roadmap for U.S. Offshore Wind: Enabling a Cost-Effective and Sustainable U.S. Offshore Wind Energy Industry Through Innovative Operations and Maintenance

The United States is currently targeting 30GW of offshore wind to be installed by 2030, and 150GW by 2050. Even considering future turbine sizes, this represents thousands of new turbines installed in a diverse set of environments, each with their unique design, installation, and maintenance challenges. While much can be learned from European and Asian experience with offshore wind over the past two decades, it is important to understand the unique circumstances of the U.S. This document explores operations and maintenance of offshore wind energy, specific to the U.S. and attempts to lay out a roadmap for needed activities to ensure reliability of future installations. The roadmap was informed through dozens of interviews with a wide cross-section of the industry, including representatives from OEMs, owner/operators, service companies, certification agencies, service providers, and researchers. The roadmap first describes the problem by component - blades, drivetrain and nacelle, structures and foundations, and electrical systems - through a look at current practices and opportunities for improvement in the areas of Failure Mode Analysis and Mitigation; Monitoring, Sensing, and Inspection; and Maintenance Execution. Crosscutting areas of Digitalization, Robotics and Automation, Prognostics and Health Management and O&M Optimization, Experimentation and Demonstration, Standardization, and Design Optimization Considering Reliability and O&M are then discussed. Finally, the roadmap summarizes all of these topics with recommendations for short (1-3 years), medium (4-7 years), and long term (8-12 years) activities, with a description of needed public and private sector contributions.

17 WIND ENERGY↗
DOE OSTI2024

Multiphysics Time-Integration for Turbulent Combustion at the Exascale

Turbulent reacting flow systems are often modeled with coupled time-dependent partial differential equations (PDEs). Solving such equations can easily tax the world's largest supercomputers. One pragmatic strategy for attacking such problems is to split the PDEs into components that can more easily be solved in isolation. This generic operator-splitting strategy leads to a set of ordinary differential equations (ODEs) that need to be solved as part of an "outer-loop" time-stepping approach. In many combustion applications, the ODEs to be solved can be very stiff, exhibiting timescales that span many orders of magnitude. The SUNDIALS library provides a plethora of robust time integration algorithms for solving these ODEs on exascale-capable computing hardware, yet for many complex applications (such multicomponent fuels or emissions predictions), the chemical models remain too complex to solve using reasonable resources. The Quasi-Steady State Approximation (QSSA) can be an effective tool for reducing the size and stiffness of the simulations. In this talk, I will discuss the use of the SUDIALS library of ODE solvers together with automatic code generation tools to solve complex turbulent reacting flow problems using QSSA models.

chemistry↗
DOE OSTI2024

A novel Bayesian approach for decomposing the radio emission of quasars: I. Modelling the radio excess in red quasars

ABSTRACT Studies show that both radio jets from the active galactic nuclei (AGNs) and the star formation (SF) activity in quasar host galaxies contribute to the quasar radio emission; yet their relative contributions across the population remain unclear. Here, we present an improved parametric model that allows us to statistically separate the SF and AGN components in observed quasar radio flux density distributions, and investigate how their relative contributions evolve with AGN bolometric luminosity ($L_\mathrm{bol}$) and redshift (z) using a fully Bayesian method. Based on the newest data from LOw-Frequency ARray Two-metre Sky Survey data release 2, our model gives robust fitting results out to $z\sim 4$, showing a quasar host galaxy SF rate (SFR) evolution that increases with bolometric luminosity and with redshift out to $z\sim 4$. This differs from the global cosmic SFR density, perhaps due to the importance of galaxy mergers. The prevalence of radio AGN emissions increases with quasar luminosity, but has little dependence on redshift. Furthermore, our new methodology and large sample size allow us to subdivide our data set to investigate the role of other parameters. Specifically, in this paper, we explore quasar colour and demonstrate that the radio excess in red quasars is due to an enhancement in AGN-related emission, since the host galaxy SF contribution to the total radio emission is independent of quasar colour. We also find evidence that this radio enhancement occurs mostly in quasars with weak or intermediate radio power.

Yue, B. -H (ORCID:0009000989352929)↗
DOE OSTI2024

Detecting rare neutral atomic-carbon absorbers with a deep neural network

ABSTRACT C i absorbers play an important role as indicators for exploring the presence of cold gas in the interstellar medium of galaxies. However, the current data base of C i absorbers is very limited due to their weak absorption feature and rarity. Here, we report results from a search of C i λλ1560, 1656 absorption lines using Mg ii absorbers as signposts with modified deep learning algorithms, which provides a very quick way to search for weak C i absorber candidates. A total of 107 C i absorbers were detected, which nearly doubles the size of previously known samples. In addition, we found 17 C i absorbers to be associated with 2175 Å dust absorbers (2DAs), i.e. about 16 per cent C i absorbers are associated with 2DAs. Comparing the average dust depletion patterns of C i absorbers with those of damped Lyman α absorbers (DLAs), Mg ii absorbers, Ca ii absorbers, and 2175 Å dust absorbers (2DAs) shows that C i absorbers generally have environments with more dust than DLAs, Mg ii, and Ca ii absorbers, but similar to dust in 2DAs. Similarity between the dust depletion pattern of C i absorbers to that of the warm disc in the Milky Way indicates that C i absorption clouds are possibly associated with disc components in distant galaxies. Therefore, C i absorbers are confirmed to be excellent probes to trace cold gas and dust in the Universe.

Ge, Jian↗
DOE OSTI2024

The Most Distant H i Galaxies Discovered by the 500 m Dish FAST

Abstract Neutral hydrogen (Hi) is the primary component of the cool interstellar medium (ISM) and is the reservoir of fuel for star formation. Owing to the sensitivity of existing radio telescopes, our understanding of the evolution of the ISM in galaxies remains limited, as it is based on only a few hundred galaxies detected in Hibeyond the local Universe. With the high sensitivity of the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we carried out a blind Hisearch, the FAST Ultra-Deep Survey, which extends to redshifts up to 0.42 and a sensitivity of 50μJy beam −1 . Here, we report the first discovery of six galaxies in Hi atz> 0.38. For these galaxies, the FAST angular resolution of ∼4′ corresponds to a mean linear size of ∼ 1.3 h 70 − 1 Mpc. These galaxies are among the most distant Hiemission detections known, with one having the most massive Hicontent ( 10 10.93 ± 0.04 h 70 − 2 M ⊙ ). Using recent data from the DESI survey and new observations with the Hale, Big Telescope Alt-azimuth, and Keck telescopes, optical counterparts are detected for all galaxies within the 3σpositional uncertainty ( 0.5 h 70 − 1 Mpc) and 200 km s −1 in recession velocity. Assuming that the dominant source of Hiis the identified optical counterpart, we find evidence of evolution in the Hicontent of galaxies over the last 4.2 Gyr. Our new high-redshift Higalaxy sample provides the opportunity to better investigate the evolution of cool gas in galaxies. A larger sample size in the future will allow us to refine our knowledge of the formation and evolution of galaxies.

Astronomy & Astrophysics↗
DOE OSTI2021

Impact of dry intrusion events on the composition and mixing state of particles during the winter Aerosol and Cloud Experiment in the Eastern North Atlantic (ACE-ENA)

Abstract. Long-range transport of continental emissions has a far-reaching influence over remote regions, resulting in substantial change in the size, morphology, and composition of the local aerosol population and cloud condensation nuclei (CCN) budget. Here, we investigate the physicochemical properties of atmospheric particles collected on board a research aircraft flown over the Azores during the winter 2018 Aerosol and Cloud Experiment in the Eastern North Atlantic (ACE-ENA) campaign. Particles were collected within the marine boundary layer (MBL) and free troposphere (FT) after long-range atmospheric transport episodes facilitated by dry intrusion (DI) events. Chemical and physical properties of individual particles were investigated using complementary capabilities of computer-controlled scanning electron microscopy and X-ray spectromicroscopy to probe particle external and internal mixing state characteristics. Furthermore, real-time measurements of aerosol size distribution, cloud condensation nuclei (CCN) concentration, and back-trajectory calculations were utilized to help bring into context the findings from offline spectromicroscopy analysis. While carbonaceous particles were found to be the dominant particle type in the region, changes in the percent contribution of organics across the particle population (i.e., external mixing) shifted from 68 % to 43 % in the MBL and from 92 % to 46 % in FT samples during DI events. This change in carbonaceous contribution is counterbalanced by the increase in inorganics from 32 % to 57 % in the MBL and 8 % to 55 % in FT. The quantification of the organic volume fraction (OVF) of individual particles derived from X-ray spectromicroscopy, which relates to the multi-component internal composition of individual particles, showed a factor of 2.06 ± 0.16 and 1.11 ± 0.04 increase in the MBL and FT, respectively, among DI samples. We show that supplying particle OVF into the κ-Köhler equation can be used as a good approximation of field-measured in situ CCN concentrations. We also report changes in the κ values in the MBL from κMBL, non-DI=0.48 to κMBL, DI=0.41, while changes in the FT result in κFT, non-DI=0.36 to κFT, DI=0.33, which is consistent with enhancements in OVF followed by the DI episodes. Our observations suggest that entrainment of particles from long-range continental sources alters the mixing state population and CCN properties of aerosol in the region. The work presented here provides field observation data that can inform atmospheric models that simulate sources and particle composition in the eastern North Atlantic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Hardenability and microstructural evolution of a precipitation strengthened Ni 50 Ti 21 Hf 25 Al 4 alloy

NiTi-based quaternary alloys are used in a variety of mechanical components, such as bearings, actuators, and dampers, owing to their good hardenability, wear resistance, and corrosion resistance. Additionally, one of the most notable characteristics of NiTi-based alloys is their shape memory effect and pseudoelastic properties. Connecting the macroscopic processing parameters employed in the design of new intermetallic alloys to the nanoscale structural characteristics dictating their behavior is crucial for improving their mechanical properties and expanding the spectrum of potential applications. Here, in this work, an arc melted Ni 50 Ti 21 Hf 25 Al 4 (at%) alloy was solution treated at 1050 °C followed by quenching and aging at 600 °C to investigate the effect of aging time on the microstructure and mechanical properties. Two types of nano-sized precipitates were observed and determined as face-centered orthorhombic H-phase (TiHf)Ni and L2 1 Heusler precipitates Ni 2 TiAl. The morphology and orientation of the H-phase were investigated using scanning and transmission electron microscopy (SEM and TEM), elucidating the coarsening kinetics and strengthening contribution of that phase to the intermetallic mechanical behavior. Following coarsening, the presence of Heusler nanoprecipitates was detected under overaged conditions through TEM imaging and nanobeam electron diffraction patterns. A peak hardness condition of 756 HV was achieved after 70 h of aging, indicating that the co-precipitation of H-phase and Heusler precipitates through a well-designed aging treatment can lead to optimal mechanical performance, thus elevating the alloy’s potential as a viable material for industrial applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A diffusion–deformation model with damage for polymer undergoing rapid decompression failure

Green hydrogen is emerging as one of the most promising alternatives to replace fossil fuels. While hydrogen gas has a good energy density by weight, its poor energy density by volume requires it to be stored under high pressure for commercial use. The hydrogen infrastructure developed to handle this high pressure hydrogen contains multiple components made of polymers as sealing agents. Although polymers do not react chemically with hydrogen gas, they undergo a mechanical failure when high pressure hydrogen gas is suddenly depressurized. This phenomenon known as 'rapid decompression failure' occurs due to the diffusion of hydrogen through polymer and getting trapped inside pre-existing cavities or voids. Here, in this paper, a continuum mechanics-based coupled diffusion-deformation-damage model was developed to predict the hydrogen distribution, stress distribution, and damage propagation inside the polymer while it undergoes rapid decompression failure. The polymer was modeled as a hyperelastic material because it represents the nonlinear material response observed in uniaxial tensile tests perfectly. The effects of hydrogen diffusivity, pre-existing cavity size, cavity location, applied hydrogen pressure, and depressurization rate on damage initiation were studied. It was found that the coefficient of diffusion plays an important role in damage initiation and damage was mostly concentrated in the inside areas rather than near the surface. Experiments were conducted with EPDM polymer which agreed well with the predicted trends using the given model. The effect of adding carbon black and silica filler particles and plasticizer to the pure EPDM polymer was also studied. It was found that damage during RDF decreases with the addition of fillers, but increases with the addition of the plasticizer. Finally, the damage evolution in the presence of two cavities was also studied, and was found that the interaction of stress fields around the cavities alters the damage occurring during RDF.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Involving the new generations in Fermilab endeavors

Since 1984 the Italian groups of the Istituto Nazionale di Fisica Nucleare (INFN) and Italian Universities, collaborating with the DOE laboratory of Fermilab (US) have been running a two-month summer training program for Italian university students. While in the first year the program involved only four physics students of the University of Pisa, in the following years it was extended to engineering students. This extension was very successful and the engineering students have been since then extremely well accepted by the Fermilab Technical, Accelerator, and Scientific Computing Division groups. Over the many years of its existence, this program has proven to be the most effective way to engage new students in Fermilab endeavors. Many students have extended their collaboration with Fermilab with their Master’s Thesis and PhD. Since 2004 the program has been supported in part by DOE in the frame of an exchange agreement with INFN. Over its almost 40 years of history, the program has grown in scope and size and has involved more than 550 Italian students from more than 20 Italian Universities, Several Institutes of Research, including ASI and INAF in Italy, and the ISSNAF Foundation in the US, have provided additional financial support. Since the program does not exclude appropriately selected non-Italian students, a handful of students from European and non-European Universities were also accepted over the years. Each intern is supervised by a Fermilab Mentor responsible for performing the training program. Training programs spanned from Tevatron, CMS, Muon (g-2), Mu2e, and Short Baseline Neutrino Experiments and DUNE design and experimental data analysis, development of particle detectors (silicon trackers, calorimeters, drift chambers, neutrino and dark matter detectors), design of electronic and accelerator components, development of infrastructures and software for exascale data handling, research on superconductive elements and on accelerating cavities, and theory of particle accelerators. Since 2010, within an extended program supported by the Italian Space Agency and the Italian National Institute of Astrophysics, a total of 30 students in physics, astrophysics, and engineering have been hosted for two months in the summer at US space science Research Institutes and laboratories. In 2015 the University of Pisa included these programs within its educational programs. Accordingly, Summer School students are enrolled at the University of Pisa for the duration of the internship and are identified and ensured as such. At the end of the internship, the students are required to write summary reports on their achievements. After positive evaluation by a University Examining Board, interns are acknowledged credits for their Diploma Supplement. The program was canceled in 2020 and 2021 due to the pandemic but restarted successfully in 2022. We believe this program can be taken as a model and easily adopted by interested institutions.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Evaluation of the excitation spectra with diffusion Monte Carlo on an auxiliary bosonic ground state

We aim to improve upon the variational Monte Carlo (VMC) approach for excitations replacing the Jastrow factor by an auxiliary bosonic (AB) ground state and multiplying it by a fermionic component factor. The instantaneous change in imaginary time of an arbitrary excitation in the original interacting fermionic system is obtained by measuring observables via the ground-state distribution of walkers of an AB system that is subject to an auxiliary effective potential. The effective potential is used to (i) drive the AB system’s ground-state configuration space toward the configuration space of the excitations of the original fermionic system and (ii) subtract from a diffusion Monte Carlo (DMC) calculation contributions that can be included in conventional approximations, such as mean-field and configuration interaction (CI) methods. In this novel approach, the AB ground state is treated statistically in DMC, whereas the fermionic component of the original system is expanded in a basis. The excitation energies of the fermionic eigenstates are obtained by sampling a fermion–boson coupling term on the AB ground state. We show that this approach can take advantage of and correct for approximate eigenstates obtained via mean-field calculations or truncated interactions. We demonstrate that the AB ground-state factor incorporates the correlations missed by standard Jastrow factors, further reducing basis truncation errors. Relevant parts of the theory have been tested in soluble model systems and exhibit excellent agreement with exact analytical data and CI and VMC approaches. In particular, for limited basis set expansions and sufficient statistics, AB approaches outperform CI and VMC in terms of basis size for the same systems. The implementation of this method in current codes, despite being demanding, will be facilitated by reusing procedures already developed for calculating ground-state properties with DMC and excitations with VMC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗