Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytic continuation”

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 361 records · Page 20

Hyperstealth dark matter and long-lived particles

A new dark matter candidate is proposed that arises as the lightest baryon from a confining 𝑆⁢𝑈⁡(𝑁) gauge theory which equilibrates with the Standard Model only through electroweak interactions. Surprisingly, this candidate can be as light as a few GeV. The lower bound arises from the intersection of two competing requirements: (i) the equilibration sector of the model must be sufficiently heavy, at least several TeV, to avoid bounds from colliders, and (ii) the lightest dark meson (that may be the dark 𝜂′, 𝜎, or the lightest glueball) has suppressed interactions with the SM and must decay before big bang nucleosynthesis. The low-energy dark sector consists of one flavor that is electrically neutral and an almost electroweak singlet. The dark matter candidate is the lightest baryon consisting of 𝑁 of these light flavors leading to a highly suppressed elastic scattering rate with the Standard Model (SM). The equilibration sector consists of vectorlike dark quarks that transform under the electroweak group, ensuring that the dark sector can reach thermal equilibrium with the SM in the early Universe. The lightest dark meson lifetimes vary between 10 −3 ≲ 𝑐⁢𝜏 ≲ 10 7 m, providing an outstanding target for LHC production and experimental detection. We delineate the interplay between the lifetime of the light mesons, the suppressed direct detection cross section of the lightest baryon, and the scale of equilibration sector that can be probed at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Transient MOC with frequency transform and DSA on unstructured mesh

We present an implementation of the transient method of characteristics (MOC) with isotropic time derivatives, accelerated by diffusion synthetic acceleration (DSA). The fully implicit frequency transform method is used to solve the transient problem with analytic precursor integration. The code works on meshes composed of almost any of the commonly used non-curvilinear finite element types, and can handle the deformation of geometry in time-dependent transport calculations. We present results of a continuous Fourier analysis for the transient multigroup DSA problem, and representative benchmarking results are presented for the C5G7-TD benchmark in 2D showing reasonable performance and agreement compared to other codes. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Comparing the performance of density functionals in describing the adsorption of atoms and small molecules on Ni(111)

Reaction networks calculated using density functional theory (DFT) play a crucial role in understanding and predicting heterogeneous catalyst performance. However, results obtained using different density functionals can vary significantly. Accordingly, it is important to assess the accuracy of available functionals in capturing adsorption energetics derived from single-crystal adsorption microcalorimetry experiments. Here, we evaluate how different density functionals describe the molecular and dissociative adsorption of CH 3 I, the adsorption of CH 3 , I, and H, and the dissociative adsorption of CH 4 on Ni(111). First, we identify energetically-preferred adsorbed states for each adsorbate at three surface coverages and calculate the corresponding enthalpies of adsorption at 160 K. Then, we compare the results to experimental measurements published in the literature. For each adsorption process we find at least one functional that is quantitatively accurate within experimental error, but no functional is accurate (within experimental error) for all adsorption processes studied. Even when assuming an additional DFT error of ±20 kJ/mol on top of the experimental error, only the PBE-D3 and RPBE-D3 functionals are accurate for all considered adsorption systems. Overall, these results indicate that quantitative agreement between density functional theory calculations and experimental measurements is both system- and functional-dependent. Importantly, single-crystal adsorption microcalorimetry experiments will continue to play an important role for benchmarking density functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Compressive Creep of Polymer Electrolyte Membranes: A Case Study for Electrolyzers

For proton-exchange membrane (PEM) water electrolyzers to be commercially feasible, PEMs must perform over long lifetimes in liquid environments under compression while maintaining mechanical stability. A hydrated environment, while inherent for operation and conductivity, undermines PEM stability. Mechanical stability of PEMs is commonly characterized in tension, which is not applicable to electrolyzers, wherein PEMs could undergo high pressures. In this study, a compression creep procedure is developed using a custom-designed setup to monitor creep response of hydrated PEMs. Furthermore, our results show that PEMs exhibit continuous creep response under compression over 24 h, with a dependence on the applied pressure and hydration state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Trace Chemical Detection Using Intercalated MXenes as a Signal Enhancing Substrate in Optical Probes

MXenes are 2D materials composed of layered transition metal nitrides or carbides. These materials are synthesized by HF exfoliation from MAX phases (Ti{sub 3}AlC{sub 2}). The 2D nanomaterial was synthesized by the removal of the 'A' element, resulting in a Mxene product (Ti{sub 3}C{sub 2}). MXenes have the general formula M{sub n+1}X{sub n}T{sub x}, where M is an early transition metal, X is Carbon and/or Nitrogen, and T accounts for surface terminated functional groups such as Fluoride, hydroxyl, and oxygen. These materials have very unique properties, similar to graphene, that allows them to be applied in a variety of trace detection techniques including surface-enhanced Raman spectroscopy (SERS). MXenes have also been demonstrated to selectively uptake uranyl ion, UO{sub 2}{sup 2+}. If this property can be combined with SERS or fluorescence detection, it may be possible to use MXenes as the basis for an alternative method to kinetic phosphorescence analysis (KPA) for trace uranyl measurements. Objectives: To confirm that MXene Nano materials are suitable substrates for SERS and sensor development by enhancing Raman signaling. To determine if certain MXene preparation methods yield materials that are more suitable for trace sensing methods. To determine uranyl uptake properties of these MXene materials and test them for analytical signals. Sample Preparation: Preparation of Ti{sub 3}C{sub 2}MXene (at FSU). MXenes were prepared by etching Al from Ti{sub 3}AlC{sub 2} (MAX phase)material. Two etching techniques yield different MXene products: LiF/HCl: Milder reaction, larger MXene flakes. HF: Harsher reaction, smaller flakes, larger layer separation. Products washed to remove etchant, vacuum filtered, and dried. Dried MXene flakes are air-stable. Film preparation for sensor testing (at SRNL): Suspend powder in diH{sub 2}O, purge with Ar, sonicate for 30 min. Centrifuge and collect supernate with suspended particles. Observed LiF-etched Mxene yielded a higher density of particles and darker collected solution. Drop-cast (4 ml) supernate onto slides and dried with Ar. For Rhodamine B (RhB) testing, drop-cast 4 ml drops onto Mxene spots and dried with Ar. Scanning Electron Microscopy conditions: 10 kV Beam energy, high vacuum; Working distance of 8 mm; beam penetration depth appx. 4 microns, beam spot size appx. 2 nanometers. Results: Detection of aluminum correlates with bright spots on image. Presence of aluminum shows that LiF/HCl etching was less thorough than HF etching. Trace Cl detection in LiF images suggests incomplete rinsing. HF has smaller feature size, more layer structure, and increased homogeneity, consistent with expectations. Macroscopic Raman spectroscopy measurements: 532 nm excitation, ∼50 mW with a ∼100 micron spot size (InPhotonics RPB probe). Kaiser Optical Holospec f/1.8 spectrometer with cooled (-60 deg.C) Andor iDus OE420 CCD. LiF 1x supernate showed good signal for trace measurements of Rhodamine B. HF and 1/4x LiF supernates showed little Mxene or Rhodamine B signal. Low deposition densities led to excess background signal from glass slides. For LiF film, response is linear with Rhodamine B concentration over range tested. Will retest with Raman microscope (∼1 micron spot size) to characterize SERS of more dilute LiF and HF etched Mxenes. Conclusions: The LiF etched material was more suitable for macroscopic SERS measurements because it was more concentrated, resulting in a thicker film than the HF etched Mxene and diluted LiF sample. However, the other materials may give greater SERS enhancements, which we hope to determine from measurements with the Raman microscope. From characterization with SEM we concluded that the HF etched Mxene is more uniform/homogenous and has smaller particle size than the LiF etched Mxene. There is still aluminum present in both samples indicating that etching wasn't complete, but the removal of the aluminum was more efficient in the HF method. Path Forward: Observe SERS with Raman microscopy, to obtain better signals for the more diluted samples and be able to compare enhancement effects for the different MXenes. Characterize uranyl sorption into MXene films and test Raman and fluorescence signals. Revisit the etching conditions to improve removal of aluminum. FSU and SRNL will continue to collaborate to create and characterize different Mxene materials and test their usefulness for sensor applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Yes, No, Maybe So: Human Factors Considerations for Fostering Calibrated Trust in Foundation Models Under Uncertainty

High-stakes analytical environments require analysts to evaluate evidence and generate conclusions to inform critical decisions often under conditions of uncertainty. Probabilistic decision-making based on incomplete or inaccurate information can reduce productivity, compromise national interests, and endanger public safety. Researchers are developing expert systems built on foundation models (FMs) to support analysts’ decision-making processes by enabling human-artificial intelligence (AI) teaming, in part through the quantification and expression of uncertainty information. As FMs continue to mature, it is imperative to correspondingly consider analysts’ needs for appropriately interpreting and using uncertainty information. However, prior research indicates that it remains unclear how analysts engage with FM-generated uncertainty information and the extent to which these interactions influence trust in, and reliance on, expert systems. We plan to review the state of the science and conduct an exploratory, qualitative study to (a) understand how properly communicated uncertainty can foster calibrated trust and appropriate reliance and (b) identify approaches for effectively conveying FM-generated uncertainty information during analytical workflows. We will administer semi-structured interviews with analysts from a specific high-stakes analytical environment to collect their current experiences with job-related uncertainty and their impressions when viewing FM-generated uncertainty information. During the interview protocol, participants will be presented with several different FM outputs and invited to discuss their thoughts and beliefs about the uncertainty information displayed. Participants may provide insights into how trust and reliance may be influenced by uncertainty. The results of this study will help us to better understand how analysts currently interpret and use uncertainty information. Our findings may inform human factors recommendations for effectively conveying uncertainty information to foster calibrated trust in, and appropriate reliance on, expert systems. Interaction designers and FM developers can use this knowledge to enhance human-AI teaming and ensure the responsible deployment of FM-based expert systems in analytical workflows.

97 MATHEMATICS AND COMPUTING↗

Sensor enabled data-driven predictive analytics for modeling and control with high penetration of DERs in distribution systems

The electric power grid is undergoing a tremendous transformation due to the increasing penetration of renewable energy resources beginning with wind and more recently with the distributed energy resources (DERs) such as solar and battery storage. DERs have dramatically changed the role of the distribution systems in the overall power grid, and they are expected to contribute a significant portion of power generation in the future. If current trends for DERs continue, system operation and control will need to change dramatically for improved grid reliability and resiliency. As renewable resources increase in penetration, new and challenging operational, planning, and design problems are expected to emerge. Some of the key challenges that arise in the planning and operation of the future grid are: 1) Quantifying the impact of high DER penetration in distribution systems on bulk grid behavior over multiple time scales. 2) Identifying whether a particular DER configuration/settings have a large impact on the overall grid behavior. These challenges can be addressed in an offline manner using detailed T&D grid models and they can also be addressed in an online manner using sensor measurements. In particular, the advancement and planned growth in sensor technology in power grid over various voltage levels provide us with a unique opportunity to tackle these challenges from a data analytic perspective without needing detailed T&D grid models. A few questions that naturally arise when addressing the challenges from DERs using sensor data are: 1) How can we use limited sensor measurements to monitor & control voltage stability and small signal stability of the bulk system? 2) How can we ensure that the developed data analytic methods are robust to data availability and quality issues? 3) How can we compute the developed analytics in a scalable manner using streaming measurements? In this project, we addressed the aforementioned challenges arising from DERs and answered the questions raised above on how to effectively use the sensor measurements to enhance the reliability and performance of the electric grid. Thus, the overarching goal of this project is to develop effective reduced/representative system models from data that make the computational complexity sufficiently manageable so as to be useful to simulate, analyze, and even control complex non-linear power systems dynamics with large penetrations of DERs. In order to achieve the objective, the project team established a four-fold technical approach 1) Formulated a combined transmission-distribution co-simulation framework for data generation and validation, 2) Derived reduced/representative models of power systems based on data-driven methods for efficient computation and appropriate representation of system behavior, 3) Developed data driven characterization of power system behavior based on transfer operator theory, machine learning and optimization for model estimation, 4) Incorporated a scalable data management and processing architecture using distributed Kafka streaming applications that coordinate input data streams to the developed data analytics. The key accomplishments of the project are: 1) Development of a scalable multi-timescale T&D co-simulation framework (both for steady state and for dynamic co-simulation) using commercial solvers (PSSE and GridLAB-D). The steady-state T&D co-simulation interface is shared with our industry partner (PJM). 2) A structured reduced order dynamic model of distribution systems that can represent partial motor stalling along with a systematic procedure to derive the model parameters. 3) A PMU based online method to monitor, localize and mitigate fault-induced delayed voltage recovery using DER reactive support and load control in distribution systems. 4) Development of linear operator based robust methodologies for dynamic state estimation, uncertainty quantification, system identification and trajectory prediction for power system dynamics. 5) An adaptive damping control for utilizing wind energy resources to provide oscillation damping and system stability. 6) Implementation of Kafka-based framework for efficient processing of streaming data using Linux-based local virtual environment.

DER integration↗

Bifurcation Analysis Reveals Solution Structures of Phase Field Models

The phase field method is playing an increasingly important role in understanding and predicting morphological evolution in materials and biological systems. Here, in this study, we develop a new analytical approach based on the bifurcation analysis to explore the mathematical solution structure of phase field models. Revealing such solution structures not only is of great mathematical interest but also may provide guidance to experimentally or computationally uncover new morphological evolution phenomena in materials undergoing electronic and structural phase transitions. To elucidate the idea, we apply this analytical approach to three representative phase field equations: the Allen-Cahn equation, the Cahn-Hilliard equation, and the Allen-Cahn-Ohta-Kawasaki system. The solution structures of these three phase field equations are also verified numerically by the homotopy continuation method.

97 MATHEMATICS AND COMPUTING↗

A-Site Cation Chemistry in Halide Perovskites

Metal halide perovskites are an important class of semiconductors now being implemented as photovoltaic absorbers and explored for light emission, among other device applications. The semiconducting properties of halide perovskites are deeply intertwined with their composition and structure. Specifically the symmetry, tilting, and distortions of the metal halide octahedra impact the band structure and other optoelectronic properties. In this review, we examine the various compositions of monovalent A-site cations in threedimensional (3D) halide perovskites AMX 3 (M = divalent metal; X = halide). We focus on how the A-site cation templates the inorganic metal-halide perovskite framework, resulting in changes in the crystal structure symmetry, as well as M–X bonding parameters, summarized in a comprehensive table of AMX 3 structures. The A-site cation motion, effects of alloying, and 2D Ruddlesden–Popper perovskite structures with unique A-site cations are further overviewed. Correlations are shown between these A-site cation dominated structural parameters and the resulting optoelectronic properties such as band gap. This review should serve as a reference for the A-site cation structural chemistry of metal halide perovskites and inspire continued research into less explored metal halide perovskite compositions and structures.

14 SOLAR ENERGY↗

Electrochemical direct air capture of CO 2 using neutral red as reversible redox-active material

Direct air capture of carbon dioxide is a viable option for the mitigation of CO 2 emissions and their impact on global climate change. Conventional processes for carbon capture from ambient air require 230 to 800 kJ thermal per mole of CO 2 , which accounts for most of the total cost of capture. Here, we demonstrate electrochemical direct air capture using neutral red as a redox-active material in an aqueous solution enabled by the inclusion of nicotinamide as a hydrotropic solubilizing agent. The electrochemical system demonstrates a high electron utilization of 0.71 in a continuous flow cell with an estimated minimum work of 35 kJ e per mole of CO 2 from 15% CO 2 . Further exploration using ambient air (410 ppm CO 2 in the presence of 20% oxygen) as a feed gas shows electron utilization of 0.38 in a continuous flow cell to provide an estimated minimum work of 65 kJ e per mole of CO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Correlation-aware binning for small-angle neutron scattering via Gaussian-process inference

Binning in small-angle neutron scattering (SANS) is typically performed empirically, with fixed parameters chosen for convenience rather than statistical optimality. Such practices often fail to balance statistical precision and spatial resolution, leading to inconsistencies across instruments and datasets. Here we establish a correlation-aware framework that determines the optimal bin width from first principles by extending the classical Freedman–Diaconis (FD) rule to account for inter-bin correlations with a Gaussian process. In this formulation, the scattering intensity is treated as a smooth stochastic field whose statistical coherence is described by a covariance matrix. Analytical expressions of errors derived from this model yield closed-form criteria that separate the total deviation into contributions from counting noise, aliasing distortion and curvature-dependent correlation effects. Expressed in reduced variables, the resulting dimensionless error surface reveals a continuous transition from the uncorrelated FD regime to the correlation-dominated limit, providing a unified description of noise suppression and resolution control. Because the formulation depends only on the profile characteristics of scattering intensity I(Q), specifically its average intensity and first- and second-order derivatives, it applies generally to any SANS measurement regardless of sample, instrument or geometry. Experimental validation using small- and ultra-small-angle neutron scattering data confirms the predicted scaling behavior, demonstrating that correlation-aware inference systematically reduces mean-squared error and enables information-efficient reproducible data reduction across materials and instruments.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗

Performance Assessments of Demand Flexibility from Grid-Interactive Efficient Buildings: Issues and Considerations

This SEE Action Network report explains basic concepts and fundamental considerations for assessing the actual demand flexibility performance of buildings participating in demand flexibility programs and responding to time-varying retail rates. Demand flexibility is the capability of distributed energy resources (DERs) to adjust a building’s load profile across different timescales. Assessments determine the timing, location, quantity, and quality of grid services provided. The results can be used for financial settlements and to improve performance of demand flexibility, support its consideration in resource potential studies and electricity system planning, and contribute to cost-effectiveness evaluations. While practitioners and regulators regularly find opportunities to improve performance assessments of demand-related services, to a large degree current best practices are sufficient for basic service offerings as demand flexibility is implemented today. However, advances in assessment practices will be required in a future with grid-interactive efficient buildings (GEBs) that provide continuous demand flexibility by integrating multiple DERs and flexibility modes (load shed, load shift, modulate, and generate). Example practices include using new baseline constructs—including whether baselines are even required for certain building flexibility modes or configurations, deploying more advanced metering and analytics, further developing cybersecurity standards, improving communication standards for increased interoperability, and establishing performance metrics and assessment procedures for load modulation as it is more fully defined and implemented. This report provides information on prioritizing and designing performance assessment elements for demand flexibility programs and time-varying retail rates, assessment protocols, and research and development needs. Additional research sponsored by DOE’s Building Technologies Office offers more detailed technical information on assessing performance of demand flexibility, including definition of performance metrics.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Characterization Results for the September 2023 Tank Farm 3H Evaporator Overhead Sample

On an annual basis, Savannah River Mission Completion (SRMC) provides 2H and 3H evaporator overhead samples to Savannah River National Lab (SRNL) to be analyzed per Section 5.2 of the Effluent Treatment Project (ETP) Waste Compliance Plan (WCP) and the Waste Acceptance Criteria (WAC). This report presents the average characterization results for the September 2023 3H evaporator overhead sample. The sample was clear and colorless with no visible solids. The results provide measurements for cesium-137 ( 137 Cs), strontium-90 ( 90 Sr), and iodine-129 ( 129 I) with the radionuclide concentration limits specified by the WAC. These analyses were performed in duplicate, and all three measured radionuclide concentrations were within ETP WAC limits. A summary of the analytical results for this 3H evaporator overhead sample includes the following: The measured cesium-137 activity in the 3H evaporator overhead sample averaged 5.70E+01 dpm/mL, which is below the ETP WAC limit of 1.30E+03 dpm/mL. As seen through historical results, Cs-137 continues to present a decreasing activity since 2020. The strontium-90 activity in the 3H evaporator overhead sample averaged <2.16E00 dpm/mL, which is below the ETP WAC limit of 1.76E+02 dpm/mL. The iodine-129 activity in the 3H evaporator overhead sample averaged <9.08E-02 dpm/mL, which is below the ETP WAC limit of 1.00E+00 dpm/mL

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Spatially quasi-periodic bifurcations from periodic traveling water waves and a method for detecting bifurcations using signed singular values

We present a method of detecting bifurcations by locating zeros of a signed version of the smallest singular value of the Jacobian. This enables the use of quadratically convergent root-bracketing techniques or Chebyshev interpolation to locate bifurcation points. Only positive singular values have to be computed, though the method relies on the existence of an analytic or smooth singular value decomposition (SVD). The sign of the determinant of the Jacobian, computed as part of the bidiagonal reduction in the SVD algorithm, eliminates slope discontinuities at the zeros of the smallest singular value. We use the method to search for spatially quasi-periodic traveling water waves that bifurcate from large-amplitude periodic waves. The water wave equations are formulated in a conformal mapping framework to facilitate the computation of the quasi-periodic Dirichlet-Neumann operator. We find examples of pure gravity waves with zero surface tension and overhanging gravity-capillary waves. In both cases, the waves have two spatial quasi-periods whose ratio is irrational. We follow the secondary branches via numerical continuation beyond the realm of linearization about solutions on the primary branch to obtain traveling water waves that extend over the real line with no two crests or troughs of exactly the same shape. The pure gravity wave problem is of relevance to ocean waves, where capillary effects can be neglected. Such waves can only exist through secondary bifurcation as they do not persist to zero amplitude. The gravity-capillary wave problem demonstrates the effectiveness of using the signed smallest singular value as a test function for multi-parameter bifurcation problems. This test function becomes mesh independent once the mesh is fine enough.

97 MATHEMATICS AND COMPUTING↗

Orientation reversal and the Chern-Simons natural boundary

We show that the fundamental property of preservation of relations, underlying resurgent analysis, provides a new perspective on crossing a natural boundary, an important general problem in theoretical and mathematical physics. This reveals a deeper rigidity aspect of resurgence in a quantum field theory path integral. The physical context here is the non-perturbative completion of complex Chern-Simons theory that associates to a 3-manifold a collection of q-series invariants labeled by Spinc structures, for which crossing the natural boundary corresponds to orientation reversal of the 3-manifold. Our new resurgent perspective leads to a practical numerical algorithm that generates q-series which are dual to unary q-series composed of false theta functions. Until recently, these duals were only known in a limited number of cases, essentially based on Ramanujan’s mock theta functions, and the common belief was that the duals might not even exist in the general case. Resurgence analysis identifies as primary objects Mordell integrals: up to changes of variables, they are Laplace transforms of resurgent functions. Their unique Borel summed transseries decomposition on either side of the Stokes line is simply the unique decomposition into real and imaginary parts. In turn, the latter are combinations of unary q-series in terms of q and its modular counterpart $\overset{\sim }{q}$ , and are resurgent by construction. The Mordell integral is analytic across the natural boundary of the q and $\overset{\sim }{q}$ series, and uniqueness of a similar decomposition which preserves algebraic relations on the other side of the boundary defines the unique boundary crossing of the q series. We demonstrate that this continuation can be efficiently implemented numerically. In the cases where unique mock modular identities are known, they are found by this numerical procedure, but the procedure can go well beyond the known list of identities. A particularly interesting feature of the resurgent approach is that it reveals new aspects, and is very different from other known approaches based on indefinite theta series, Appell-Lerch sums, and representation theory of logarithmic vertex operator algebras.

Chern-Simons theories↗

Tungsten oxide-coated copper gallium selenide sustains long-term solar hydrogen evolution

This work demonstrates that ultrathin (4 nm) tungsten oxide (WO 3 ) coatings on copper gallium selenide (CuGa 3 Se 5 ) photocathodes have the potential for long-term solar hydrogen evolution. We describe how with a combination of a robust 1.84 eV CuGa 3 Se 5 absorber layer, a WO 3 protective coating, and a Pt catalyst, we obtain a new durability milestone for any non-silicon photoelectrochemical hydrogen-producing device by passing 21,490 C cm -2 of charge across six weeks of continuously-illuminated chronoamperometric testing under applied bias.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

How good are polarizable and flexible models for water: Insights from a many-body perspective

We present a systematic analysis of state-of-the-art polarizable and flexible water models from a many-body perspective, with a specific focus on their ability to represent the Born–Oppenheimer potential energy surface of water from the gas to the liquid phase. Using coupled cluster data in the completed basis set limit as a reference, we examine the accuracy of the polarizable models in reproducing individual many-body contributions to interaction energies and harmonic frequencies of water clusters and compare their performance with that of MB-pol, an explicit many-body model that has been shown to correctly predict the properties of water across the entire phase diagram. Based on these comparisons, we use MB-pol as a reference to analyze the ability of the polarizable models to reproduce the energy landscape of liquid water under ambient conditions. We find that, while correctly reproducing the energetics of minimum-energy structures, the polarizable models examined in this study suffer from inadequate representations of many-body effects for distorted configurations. To investigate the role played by geometry-dependent representations of 1-body charge distributions in reproducing coupled cluster data for both interaction and many-body energies, we introduce a simplified version of MB-pol that adopts fixed atomic charges and demonstrate that the new model retains the same accuracy as the original MB-pol model. Based on the analyses presented in this study, we believe that future developments of both polarizable and explicit many-body models should continue in parallel and would benefit from synergistic efforts aimed at integrating the best aspects of the two theoretical/computational frameworks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Communication—Electrodeposition of Indium Directly on Silicon

Direct electrodeposition of indium onto silicon paves the way for advances in microelectronics, photovoltaics, and optoelectronics. Indium is generally electrodeposited onto silicon utilizing a physically or thermally deposited metallic seed layer. Eliminating this layer poses benefits in microelectronics by reducing resistive interfaces and in vapor-liquid-solid conversion to III-V material by allowing direct contact to the single-crystal silicon substrate for epitaxial conversion. We investigated conditions to directly electrodeposit indium onto n-type Si(100). We illustrate that a two-step galvanostatic plating at low temperatures can consistently produce smooth, continuous films of indium over large areas, in bump morphologies, and conformally into inverted pyramids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗