Results for “static models” — Page 43

DOE OSTI2014

Modeling CO 2 Sequestration in Saline Aquifer and Depleted Oil Reservoir To Evaluate Regional CO 2 Sequestration Potential of Ozark Plateau Aquifer System, South-Central Kansas

1. Drilled, cored, and logged three wells to the basement and collecting more than 2,700 ft of conventional core; obtained 20 m i2 of multicomponent 3D seismic imaging and merged and reprocessed more than 125 mi 2 of existing 3D seismic data for use in modeling CO 2 - EOR oil recovery and CO 2 storage in five oil fields in southern Kansas. 2. Determined the technical feasibility of injecting and sequestering CO 2 in a set of four depleted oil reservoirs in the Cutter, Pleasant Prairie South, Eubank, and Shuck fields in southwest Kansas; of concurrently recovering oil from those fields; and of quantifying the volumes of CO2 sequestered and oil recovered during the process. 3. Formed a consortium of six oil operating companies, five of which own and operate the four fields. The consortium became part of the Southwest Kansas CO 2 -EOR Initiative for the purpose of sharing data, knowledge, and interest in understanding the potential for CO 2 -EOR in Kansas. 4. Built a regional well database covering 30,000 mi 2 and containing stratigraphic tops from ~90,000 wells; correlated 30 major stratigraphic horizons; digitized key wells, including wireline logs and sample logs; and analyzed more than 3,000 drill stem tests to establish that fluid levels in deep aquifers below the Permian evaporites are not connected to the surface and therefore pressures are not hydrostatic. Connectivity with the surface aquifers is lacking because shale aquitards and impermeable evaporite layers consist of both halite and anhydrite. 5. Developed extensive web applications and an interactive mapping system that do the following: a. Facilitate access to a wide array of data obtained in the study, including core descriptions and analyses, sample logs, digital (LAS) well logs, seismic data, gravity and magnetics maps, structural and stratigraphic maps, inferred fault traces, earthquakes, Class I and II disposal wells, and surface lineaments. b. Provide real-time analysis of the project dataset, including automated integration and viewing of well logs, core, core analyses, brine chemistry, and stratigraphy using the Java Profile app. A cross-section app allows for the display of log data for up to four wells at a time. 6. Integrated interpretations from the project’s interactive web-based mapping system to gain insights to aid in assessing the efficacy of geologic CO 2 storage in Kansas and insights toward understanding recent seismicity to aid in evaluating induced vs. naturally occurring earthquakes. 7. Developed a digital type-log system, including web-based software to modify and refine stratigraphic nomenclature to provide stakeholders a common means for communication about the subsurface. 8. Contracted use of a nuclear magnetic resonance (NMR) log and ran it slowly to capture response and characterize larger pores common for carbonate reservoirs. Used NMR to extend core analyses to apply permeability, relative permeability to CO 2 , and capillary pressure to the major rock types, each uniquely expressed as a reservoir quality index (RQI), present in the Mississippian and Arbuckle rocks. 9. Characterized and evaluated the possible role of microbes in dense brines. Used microbes to compliment H/O stable isotopes to fingerprint brine systems. Used perforation/swabbing to obtain samples from multiple hydrostratigraphic units and confirmed equivalent results using less expensive drill stem tests (DST). 10. Used an integrated approach from whole core, logs, tests, and seismic to verify and quantify properties of vuggy, brecciated, and fractured carbonate intervals. 11. Used complex geocellular static and dynamic models to evaluate regional storage capacity using large parallel processing. 12. Carbonates are complex reservoirs and CO 2 -EOR needs to move to the next generation to increase effectiveness of CO 2 and efficiency and safety of the injection.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Modeling CO 2 Sequestration in Saline Aquifer and Depleted Oil Reservoir To Evaluate Regional CO 2 Sequestration Potential of Ozark Plateau Aquifer System, South-Central Kansas

1. Drilled, cored, and logged three wells to the basement and collecting more than 2,700 ft of conventional core; obtained 20 m i2 of multicomponent 3D seismic imaging and merged and reprocessed more than 125 mi 2 of existing 3D seismic data for use in modeling CO 2 - EOR oil recovery and CO 2 storage in five oil fields in southern Kansas. 2. Determined the technical feasibility of injecting and sequestering CO 2 in a set of four depleted oil reservoirs in the Cutter, Pleasant Prairie South, Eubank, and Shuck fields in southwest Kansas; of concurrently recovering oil from those fields; and of quantifying the volumes of CO2 sequestered and oil recovered during the process. 3. Formed a consortium of six oil operating companies, five of which own and operate the four fields. The consortium became part of the Southwest Kansas CO 2 -EOR Initiative for the purpose of sharing data, knowledge, and interest in understanding the potential for CO 2 -EOR in Kansas. 4. Built a regional well database covering 30,000 mi 2 and containing stratigraphic tops from ~90,000 wells; correlated 30 major stratigraphic horizons; digitized key wells, including wireline logs and sample logs; and analyzed more than 3,000 drill stem tests to establish that fluid levels in deep aquifers below the Permian evaporites are not connected to the surface and therefore pressures are not hydrostatic. Connectivity with the surface aquifers is lacking because shale aquitards and impermeable evaporite layers consist of both halite and anhydrite. 5. Developed extensive web applications and an interactive mapping system that do the following: a. Facilitate access to a wide array of data obtained in the study, including core descriptions and analyses, sample logs, digital (LAS) well logs, seismic data, gravity and magnetics maps, structural and stratigraphic maps, inferred fault traces, earthquakes, Class I and II disposal wells, and surface lineaments. b. Provide real-time analysis of the project dataset, including automated integration and viewing of well logs, core, core analyses, brine chemistry, and stratigraphy using the Java Profile app. A cross-section app allows for the display of log data for up to four wells at a time. 6. Integrated interpretations from the project’s interactive web-based mapping system to gain insights to aid in assessing the efficacy of geologic CO 2 storage in Kansas and insights toward understanding recent seismicity to aid in evaluating induced vs. naturally occurring earthquakes. 7. Developed a digital type-log system, including web-based software to modify and refine stratigraphic nomenclature to provide stakeholders a common means for communication about the subsurface. 8. Contracted use of a nuclear magnetic resonance (NMR) log and ran it slowly to capture response and characterize larger pores common for carbonate reservoirs. Used NMR to extend core analyses to apply permeability, relative permeability to CO 2 , and capillary pressure to the major rock types, each uniquely expressed as a reservoir quality index (RQI), present in the Mississippian and Arbuckle rocks. 9. Characterized and evaluated the possible role of microbes in dense brines. Used microbes to compliment H/O stable isotopes to fingerprint brine systems. Used perforation/swabbing to obtain samples from multiple hydrostratigraphic units and confirmed equivalent results using less expensive drill stem tests (DST). 10. Used an integrated approach from whole core, logs, tests, and seismic to verify and quantify properties of vuggy, brecciated, and fractured carbonate intervals. 11. Used complex geocellular static and dynamic models to evaluate regional storage capacity using large parallel processing. 12. Carbonates are complex reservoirs and CO 2 -EOR needs to move to the next generation to increase effectiveness of CO 2 and efficiency and safety of the injection.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Light hydrocarbon and noble gas migration as an analogue for potential CO 2 leakage: numerical simulations and field data from three hydrocarbon systems

Abstract Sites where anthropogenic CO 2 captured from industrial sources is stored in deep geological formations for climate change mitigation are required to show secure retention of the injected CO 2 . Monitoring, reporting, and verification (MRV) plans are needed to indicate that no CO 2 release has occurred. We explored the degree to which direct comparison between a surface anomaly and reservoir geochemistry using various geochemical parameters can be used for attribution. We used data collected on light hydrocarbons and noble gases throughout the sedimentary column at three CO 2 enhanced oil recovery (EOR) sites to understand the processes that may cause fluid evolution. Light hydrocarbon and noble gases were sampled from reservoirs, gas‐bearing intervals above reservoirs, and groundwater. Vertical profiles indicated that lighter components are relatively enriched during migration (i.e. chromatographic separation). Static and numerical models were designed to simulate episodic gas migration and geochemical alteration of these geochemical parameters from solubility and sorption. The effects of hydrocarbon solubility were minimal (Bernard ratio changes within 5.2%) although field data were within the range of expected alterations from sorption. Forward models of CO 2 migration and noble gas interactions showed that CO 2 stripping causes an enrichment of crustal noble gases. In areas where natural fluxes of CO 2 from depth are non‐existent, the occurrence of crustal noble gas signature may distinguish fugitive CO 2 from the reservoir from natural near‐surface sources, and could be considered to explain apparent fluid anomalies. © 2019 Society of Chemical Industry and John Wiley & Sons, Ltd.

Anderson, Jacob↗

DOE OSTI2023

A dynamic model of wind turbine yaw for active farm control

This paper presents a graph-based dynamic yaw model to predict the dynamic response of the hub-height velocities and the power of a wind farm to a change in yaw. The model builds on previous work where the turbines define the nodes of the graph and the edges represent the interactions between turbines. Advances associated with the dynamic yaw model include a novel analytical description of the deformation of wind turbine wakes under yaw to represent the velocity deficits and a more accurate representation of the interturbine travel time of wakes. The accuracy of the model is improved by coupling it with time- and space-dependent estimates of the wind farm inflow based on real-time data from the wind farm. The model is validated both statically and dynamically using large-eddy simulations. An application of the model is presented that incorporates the model into an optimal control loop to control the farm power output.

17 WIND ENERGY↗

DOE OSTI2018

An asymptotically compatible meshfree quadrature rule for nonlocal problems with applications to peridynamics

In this paper, we present a meshfree quadrature rule for compactly supported nonlocal integro-differential equations (IDEs) with radial kernels. We apply this rule to develop a meshfree discretization of a peridynamic solid mechanics model that requires no background mesh. Existing discretizations of peridynamic models have been shown to exhibit a lack of asymptotic compatibility to the corresponding linearly elastic local solution. By posing the quadrature rule as an equality constrained least squares problem, we obtain asymptotically compatible convergence by introducing polynomial reproduction constraints. Our approach naturally handles traction-free conditions, surface effects, and damage modeling for both static and dynamic problems. We demonstrate high-order convergence to the local theory by comparing to manufactured solutions and to cases with crack singularities for which an analytic solution is available. We verify the applicability of the approach to realistic problems by reproducing high-velocity impact results from the Kalthoff–Winkler experiments.

42 ENGINEERING↗

DOE OSTI2021

A workflow to assess the efficacy of brine extraction for managing injection-induced seismicity potential using data from a CO 2 injection site near Decatur, Illinois

Injection of CO 2 for storage in a geologic formation increases pore pressure and alters in situ stresses. Depending on the orientation of any existing fault and fracture planes, such as critically stressed planes, this stress alteration will modify normal stresses acting on planes and could result in frictional sliding and release stored energy in the form of seismicity. Brine extraction (BE) is a technique that can be applied prior to, or during, CO 2 injection to reduce pore pressure for increasing storage capacity and, potentially, for reducing the likelihood of frictional sliding. Here a workflow is described to assess the efficacy of BE for mitigating frictional sliding (i.e., seismicity) during injection and entails: site characterization, stress calculations and failure assessment, static and dynamic modeling, and BE operational planning. Site characterization describes the stress field used to calculate the Coulomb Failure Function (CFF) that constrains allowable pore pressure changes and injection rates in the numerical simulation of CO 2 injection scenarios. The inclusion of BE in the workflow allows for determination of the potential need for pressure reduction, and evaluation of the effectiveness of this operation. Example application of the workflow using an injection field dataset near Decatur, IL, provides insight on fracture planes and stresses at the site, formation properties and the impact of variable CO 2 injection-rate targets on whether BE plans are required. The study workflow indicates that BE could enhance CO 2 injection rate by 39% and correspondingly reduce the potential for injection-induced seismicity as indicated by a reduction in CFF.

58 GEOSCIENCES↗

DOE OSTI2021

A conservative level set method for N-phase flows with a free-energy-based surface tension model

Not Available

Howard, Amanda A. (ORCID:0000000264116198)↗

DOE OSTI2019

Response of deep aquifers to climate variability

There is a general agreement that deep aquifers experience significant lag time in their response to climatic variations. Analysis of Temporal Gravity Recovery and Climate Experiment (GRACE), Soil Moisture and Ocean Salinity mission (SMOS), satellite altimetry, stable isotopic composition of groundwater, and precipitation and static global geopotential models over the Nubian Sandstone Aquifer System (NSAS) revealed rapid aquifer response to climate variability. Findings include: (1) The recharge areas of the NSAS (Northern Sudan Platform subbasin) witnessed a dry period (2002-2012), where average annual precipitation (AAP) was modest (85 mm) followed by a wet period (2013-2016; AAP: 107 mm), and during both periods the APP remained negligible (<10 mm) over the northern parts of the NSAS (Dakhla subbasin); (2) the secular trends in terrestrial water storage (TWS) over the Dakhla subbasin were estimated at -3.8 ± 1.3 mm/yr and + 7.8 ± 1 mm/yr for the dry and wet periods, respectively; (3) spatial variations in TWS values and phase are consistent with rapid groundwater flow from the Northern Sudan Platform subbasin and Lake Nasser towards the Dakhla subbasin during the wet period and from the lake during the dry period; and (4) networks of densely fractured and karstified bedrocks provide preferential pathways for groundwater flow. The proposed model is supported by (1) rapid response in groundwater levels in distant wells (>280 km from source areas) and in soil moisture content in areas with shallow (<2 m) groundwater levels to fluctuations in Lake Nasser surface water, and (2) the isotopic composition (O, H) of groundwater along the preferred pathways, consistent with mixing of enriched (Lake Nasser water or precipitation over Sudan) and depleted (NSAS fossil water) end-members. Here, findings provide new insights into the response of large, deep aquifers to climate variability and address the sustainability of the NSAS and similar fossil aquifers worldwide.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate

Characterizing the properties of groundwater aquifers is essential for predicting aquifer response and managing groundwater resources. In this work, we develop a high-dimensional scalable Bayesian inversion framework governed by a three-dimensional quasi-static linear poroelastic model to characterize lateral permeability variations in groundwater aquifers. In this work, we determine the maximum a posteriori (MAP) point of the posterior permeability distribution from centimeter-level surface deformation measurements obtained from Interferometric Synthetic Aperture Radar (InSAR). The scalability of our method to high parameter dimension is achieved through the use of adjoint-based derivatives, inexact Newton methods to determine the MAP point, and a Mat´ern class sparse prior precision operator. Together, these guarantee that the MAP point is found at a cost, measured in number of forward/adjoint poroelasticity solves, that is independent of the parameter dimension. We apply our methodology to a test case for a municipal well in Mesquite, Nevada, in which InSAR and GPS surface deformation data are available. We solve problems with up to 320,824 state variable degrees of freedom (DOFs) and 16,896 parameter DOFs. A consistent treatment of noise level is employed so that the aquifer characterization result does not depend on the pixel spacing of surface deformation data. Our results show that the use of InSAR data significantly improves characterization of lateral aquifer heterogeneity, and the InSAR-based aquifer characterization recovers complex lateral displacement trends observed by independent daily GPS measurements.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Prediction of chronic kidney disease progression using recurrent neural network and electronic health records

Chronic kidney disease (CKD) is a progressive loss in kidney function. Early detection of patients who will progress to late-stage CKD is of paramount importance for patient care. To address this, we develop a pipeline to process longitudinal electronic heath records (EHRs) and construct recurrent neural network (RNN) models to predict CKD progression from stages II/III to stages IV/V. The RNN model generates predictions based on time-series records of patients, including repeated lab tests and other clinical variables. Our investigation reveals that using a single variable, the recorded estimated glomerular filtration rate (eGFR) over time, the RNN model achieves an average area under the receiver operating characteristic curve (AUROC) of 0.957 for predicting future CKD progression. When additional clinical variables, such as demographics, vital information, lab test results, and health behaviors, are incorporated, the average AUROC increases to 0.967. In both scenarios, the standard deviation of the AUROC across cross-validation trials is less than 0.01, indicating a stable and high prediction accuracy. Our analysis results demonstrate the proposed RNN model outperforms existing standard approaches, including static and dynamic Cox proportional hazards models, random forest, and LightGBM. The utilization of the RNN model and the time-series data of previous eGFR measurements underscores its potential as a straightforward and effective tool for assessing the clinical risk of CKD patients concerning their disease progression.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2013

Design, analysis and testing of a parallel-kinematic high-bandwidth XY nanopositioning stage

This paper presents the design, analysis, and testing of a parallel-kinematic high-bandwidth XY nanopositioning stage driven by piezoelectric stack actuators. The stage is designed with two kinematic chains. In each kinematic chain, the end-effector of the stage is connected to the base by two symmetrically distributed flexure modules, respectively. Each flexure module comprises a fixed-fixed beam and a parallelogram flexure serving as two orthogonal prismatic joints. With the purpose to achieve high resonance frequencies of the stage, a novel center-thickened beam which has large stiffness is proposed to act as the fixed-fixed beam. The center-thickened beam also contributes to reducing cross-coupling and restricting parasitic motion. To decouple the motion in two axes totally, a symmetric configuration is adopted for the parallelogram flexures. Based on the analytical models established in static and dynamic analysis, the dimensions of the stage are optimized in order to maximize the first resonance frequency. Then finite element analysis is utilized to validate the design and a prototype of the stage is fabricated for performance tests. According to the results of static and dynamic tests, the resonance frequencies of the developed stage are over 13.6 kHz and the workspace is 11.2 μm × 11.6 μm with the cross-coupling between two axes less than 0.52%. It is clearly demonstrated that the developed stage has high resonance frequencies, a relatively large travel range, and nearly decoupled performance between two axes. For high-speed tracking performance tests, an inversion-based feedforward controller is implemented for the stage to compensate for the positioning errors caused by mechanical vibration. The experimental results show that good tracking performance at high speed is achieved, which validates the effectiveness of the developed stage.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2014

Accurate modeling of switched reluctance machine based on hybrid trained WNN

According to the strong nonlinear electromagnetic characteristics of switched reluctance machine (SRM), a novel accurate modeling method is proposed based on hybrid trained wavelet neural network (WNN) which combines improved genetic algorithm (GA) with gradient descent (GD) method to train the network. In the novel method, WNN is trained by GD method based on the initial weights obtained per improved GA optimization, and the global parallel searching capability of stochastic algorithm and local convergence speed of deterministic algorithm are combined to enhance the training accuracy, stability and speed. Based on the measured electromagnetic characteristics of a 3-phase 12/8-pole SRM, the nonlinear simulation model is built by hybrid trained WNN in Matlab. The phase current and mechanical characteristics from simulation under different working conditions meet well with those from experiments, which indicates the accuracy of the model for dynamic and static performance evaluation of SRM and verifies the effectiveness of the proposed modeling method.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2015

Wetting properties of molecularly rough surfaces

We employ molecular dynamics simulations to study the wettability of nanoscale rough surfaces in systems governed by Lennard-Jones (LJ) interactions. We consider both smooth and molecularly rough planar surfaces. Solid substrates are modeled as a static collection of LJ particles arranged in a face-centered cubic lattice with the (100) surface exposed to the LJ fluid. Molecularly rough solid surfaces are prepared by removing several strips of LJ atoms from the external layers of the substrate, i.e., forming parallel nanogrooves on the surface. We vary the solid-fluid interactions to investigate strongly and weakly wettable surfaces. We determine the wetting properties by measuring the equilibrium droplet profiles that are in turn used to evaluate the contact angles. Macroscopic arguments, such as those leading to Wenzel’s law, suggest that surface roughness always amplifies the wetting properties of a lyophilic surface. However, our results indicate the opposite effect from roughness for microscopically corrugated surfaces, i.e., surface roughness deteriorates the substrate wettability. Adding the roughness to a strongly wettable surface shrinks the surface area wet with the liquid, and it either increases or only marginally affects the contact angle, depending on the degree of liquid adsorption into the nanogrooves. For a weakly wettable surface, the roughness changes the surface character from lyophilic to lyophobic due to a weakening of the solid-fluid interactions by the presence of the nanogrooves and the weaker adsorption of the liquid into the nanogrooves.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Ultrafast excited-state deactivation of 9-methylhypoxanthine in aqueous solution: A QM/MM MD study

Photoinduced ultrafast non-adiabatic decay of 9-methylhypoxanthine (9MHPX) in aqueous solution was investigated by ab initio surface-hopping dynamics calculations using a combined quantum mechanical/molecular mechanical approach. The absorption spectra of 9MHPX in aqueous solution were also explored by the hybrid cluster-continuum model at the level of time-dependent density functional theory along with the polarizable continuum model (PCM). The static electronic-structure calculations indicate that the absorption spectra of 9MHPX simulated by TD-B3LYP/PCM and TD-X3LYP/PCM can reproduce very well the experimental findings, with the accuracy of about 0.20 eV. According to dynamics simulations, irradiation of 9MHPX populates the bright excited singlet S{sub 1} state, which may undergo an ultrafast non-radiative deactivation to the S{sub 0} state. The lifetime of the S{sub 1} state of 9MHPX in aqueous solution is predicted to be 115.6 fs, slightly longer than that in the gas phase (88.8 fs), suggesting that the solvent water has no significant influence on the excited-state lifetime of 9MHPX. Such a behavior in 9MHPX is distinctly different from its parent hypoxanthine keto-N9H tautomer in which the excited-state lifetime of the latter in water solution was remarkably enhanced as compared to the gas phase. The significant difference of the photodynamical behaviors between 9MHPX and keto-N9H can be ascribed to their different hydrogen bond environment in aqueous solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Explosive dispersal of particles in high speed environments

In this paper, we present the results of the explosive dispersal of particles in high-speed environments. We carry out Euler–Lagrange numerical simulations of a source at quiescent ambient conditions as well as moving at Mach numbers of 3 and 6. Particle volume fractions of 0%, 1%, and 4.5% are presented. The detonation profile is computed with the Jones–Wilkins–Lee equation of state using a reactive burn model. Non-static cases provide a framework to consider the effect of a bow shock and pre-existing high-speed flow conditions on the dispersal process. We also compute averages of both static and dynamic pressures, as well as impulse density histories on virtual probe planes to characterize the momentum of the flow and particles that would deposit on a target. Results suggest that the presence of the particles can have a substantial effect on the pressure average of the virtual target planes.

Durant, B. A.↗

DOE OSTI2013

SUPERMASSIVE BLACK HOLE BINARY EVOLUTION IN AXISYMMETRIC GALAXIES: THE FINAL PARSEC PROBLEM IS NOT A PROBLEM

During a galaxy merger, the supermassive black hole (SMBH) in each galaxy is thought to sink to the center of the potential and form an SMBH binary; this binary can eject stars via three-body scattering, bringing the SMBHs ever closer. In a static spherical galaxy model, the binary stalls at a separation of about a parsec after ejecting all the stars in its loss cone-this is the well-known final parsec problem. Earlier work has shown that the centrophilic orbits in triaxial galaxy models are key in refilling the loss cone at a high enough rate to prevent the black holes from stalling. However, the evolution of binary SMBHs has never been explored in axisymmetric galaxies, so it is not clear if the final parsec problem persists in these systems. Here we use a suite of direct N-body simulations to follow SMBH binary evolution in galaxy models with a range of ellipticity. For the first time, we show that mere axisymmetry can solve the final parsec problem; we find the SMBH evolution is independent of N for an axis ratio of c/a = 0.8, and that the SMBH binary separation reaches the gravitational radiation regime for c/a = 0.75.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Importance of dynamic lattice effects for crystal field excitations in the quantum spin ice candidate Pr 2 Zr 2 O 7

We explore dynamic interactions between the crystal lattice and magnetic degrees of freedom in a frustrated magnetic system using the example of a pyrochlore quantum spin-ice candidate Pr 2 Zr 2 O 7 . Using Raman scattering spectroscopy we demonstrate that crystal electric field excitations of Pr3+, which define the magnetic properties of Pr 2 Zr 2 O 7 , cannot be understood within a model of a static lattice. We identify vibronic interactions with a phonon which lead to a splitting of a doublet crystal field excitation at around 55 meV. We also observe an unconventional behavior of a splitting of the non-Kramers ground state doublet of Pr 3+ , revealed by observing excitations to the first excited singlet state $E$$^{0}_{g}$ → A 1 g at around 10 meV. The splitting has a strong temperature dependence, where the doublet structure is most prominent between 50 and 100 K, and the weight of one of the components strongly decreases on cooling contrary to simple thermal population tendency. Furthermore, we suggest a static or dynamic deviation of Pr 3+ from the position in the ideal crystal structure can be the origin of the effect, with the deviation strongly decreasing at low temperatures.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Cooperative mechanisms behind nanoscale smoothing on metal surfaces: From adatom diffusion to step nucleation

Many technologies require stable or metastable surface morphology. Here, in this paper we study the factors that control the metastability of a common feature of rough surfaces: “hillocks.” We use low energy electron microscopy to follow the evolution of the individual atomic steps in hillocks on Pd(111). We show that the uppermost island in the stack often adopts a static, metastable configuration. Modeling this result shows that the degree of the metastability depends on the configuration of steps dozens of atomic layers lower. Our model allows us to link surface metastability to the atomic processes of surface evolution.

Fernandes Cauduro, Andre Luis↗