Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface shape fitting”

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 55 records · Page 3

A method to predict texture effect on ion beam channeling analysis of polycrystals and the application to study the mosaic spreading effect in highly oriented pyrolytic graphite

We propose a method to convert the channeling Rutherford backscattering spectrum yield map of a single crystal to a polycrystal through a matrix rotation technique. The rotation matrix is determined by the deviation of the crystal axial direction from the original z axis. The final yield map is created after averaging the rotated yields using the texture function as the weight factor. For highly oriented pyrolytic graphite (HOPG) exhibiting mosaic spread, the method leads to a Gaussian kernel averaging of the map obtained from a single crystal. The yield map of a single crystal is obtained by a simulation of ion trajectories in a potential field described by Moliere screened Coulomb potentials. Yield maps are calculated under various σ values (standard deviations of mosaic spread). The simulated results are compared with experimental results obtained using 1.2 MeV alpha particle. σ is extracted through the best fitting, demonstrating that the method can be used to obtain texture details. The effects of mosaic spread on minimum yield χ min and the half-width at half maximum of angular scans ψ ½ are systematically modeled and compared with previous theoretical equations. The study also shows that previous theoretical equations are valid only at small σ values. The proposed method can be applied to any type of polycrystal and is not limited to HOPG. It provides near-surface mosaic spread and crystallography information with a longitudinal depth resolution of tens of nanometers and is not influenced by grain shapes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Enhanced coercivity in Fe5C2/SiO2 core/shell nanocrystals

Rod-shaped Fe5C2 and core/shell Fe5C2/SiO2 nanocrystals were synthesized via a solution-based chemical method. Structural analysis confirmed the monoclinic phase of Fe5C2 with space group C2/c. Zero-field-cooling (ZFC) and field-cooling (FC) magnetization curves revealed distinct magnetic behaviors: uncoated Fe5C2 exhibited a low-temperature FC plateau indicative of strong dipolar interactions, while Fe5C2/SiO2 showed a monotonic increase in FC magnetization, suggesting reduced dipolar interactions due to SiO2 surface passivation. Isothermal remanent magnetization (IRM) and DC demagnetization (DCD) measurements supported this trend, with δM plots confirming weaker dipolar interactions in the coated sample. Bloch’s law fitting of temperature-dependent saturation magnetization showed a smaller Bloch’s constant for pure Fe5C2 and a larger value for Fe5C2/SiO2, reflecting enhanced surface disorder and reduced exchange coupling in the latter. Notably, Fe5C2/SiO2 demonstrated increased coercivity, attributed to decreased dipolar interaction and elevated surface anisotropy. Kneller’s law fitting yielded higher blocking temperatures for Fe5C2 (476 K) than Fe5C2/SiO2 (456 K), highlighting the impact of dipolar interactions on magnetic relaxation. These findings illustrate how SiO2 coatings effectively modulate dipolar interactions and enhance coercivity in Fe5C2 nanocrystals.

Joshi, Pramanand [Department of Physics, Universit↗

Disentangling the Impacts of Microtopography and Shrub Distribution on Snow Depth in a Subarctic Watershed: Toward a Predictive Understanding of Snow Spatial Variability

Snow plays a critical role in carbon cycling, vegetation dynamics, and permafrost hydrology at high latitudes by influencing surface energy exchange. Predicting snow distribution patterns is essential for understanding the evolution of Arctic ecosystems, yet scaling process-level knowledge to landscape predictions remains challenging. Here, we analyze snow depth (2019 and 2022), terrain elevation, and vegetation height from a watershed on the Seward Peninsula, Alaska, to examine how topography and shrubs shape snow redistribution across spatial scales. We find that snow depth is strongly coupled to terrain at scales below ∼60 m but becomes increasingly decoupled at larger scales. The topographic model of snow depth variation, which transforms terrain data to align with these scale-dependent snow patterns, is well correlated with local snow depth variations (linear fit R 2 > 0.5 for 85% of 100-m patches). A machine learning reconstruction of shrub canopy snow trapping reveals a simple exponential relationship between canopy structure and snow accumulation ( R 2 = 0.59), highlighting the combined influence of topography and vegetation on snow distribution. Together, these empirical relationships capture much of the observed snow variability in the watershed ( R 2 = 0.49, root mean square error (RMSE) = 30 cm), though systematic limitations persist in areas of strong scour and at coarser scales where wind-terrain interactions are more complex. These findings provide a framework for more efficient snow depth prediction and offer insights to improve snow-vegetation feedback representation in Earth System Models.

54 ENVIRONMENTAL SCIENCES↗

BISON fuel performance modeling optimization for experiment X447 and X447A using axial swelling and cladding strain measurements

With the recent need to qualify new reactor designs such as the Versatile Test Reactor (VTR), fuel performance calculations need to be performed to determine safety criteria of the proposed designs. In order to validate the fuel performance results obtained by a fuel performance code, BISON, for new reactor designs, legacy fuel from EBR-II and FFTF MFF with Post -Irradiation Examination (PIE) data need to be used as validation cases to benchmark models. Here in this work, BISON has been paired with the Fuels Irradiation & Physics Database (FIPD) and IFR Materials Information System (IMIS) to supply PIE data for comparison with simulations of EBR-II experiments X447/X447A. X447/X447A were assessed by implementing models for Fuel Cladding Chemical Interaction (FCCI) within BISON and optimizing the friction coefficient between the fuel surface and the cladding, the anisotropic swelling factor, and the HT9 first thermal creep scalar (which scales the first term in the HT9 creep equation) to best match the PIE axial fuel swelling height and cladding profilometry for all pins in X447/X447A. The optimal values were found using a generic algorithm developed to select different values for the three parameters until end criteria was met and error couldn’t be reduced further. The BISON-simulated cladding profilometry was evaluated using Standard Error of the Estimate (SEE) to account for the profile shape of the cladding profilometry. Optimal values for the friction coefficient, anisotropic fuel swelling factor, and HT9 first thermal creep scalar were found to best fit the BISON simulation results to the PIE measurements found in IMIS and FIPD. Improvements to current models are suggested to account for the underprediction of fuel swelling at low burnups and the overprediction of fuel swelling at higher burnups observed for the axial fuel swelling height. Although two pins in EBR-II X447/X447A (DP70 and DP75) were known to fail due to FCCI, none of the pins simulated in BISON reached a cumulative damage fraction (CDF) above 0.008 with FCCI correlations coupled in the BISON simulations. The error estimate generated for all pins in X447/X447A using optimal values was 209 µm, which is deemed acceptable.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Carbon nanotube-based mixed-matrix membranes with supramolecularly engineered interface for enhanced gas separation performance

Here in this study, novel nanocomposite membranes were fabricated from single-walled carbon nanotubes (SWNTs) and a triptycene-containing polyimide. The desirable interfacial morphology and homogenous dispersion of SWNTs were achieved by strong π-π stacking and supramolecular shape-fitting interactions between the nanotube and the paddlewheel-shaped triptycene moieties in polymer matrix. Mixed-matrix membranes (MMMs) containing 2–15 wt% of as-purchased SWNTs (AP-SWNTs) were prepared and tested for their gas transport properties. While all MMMs showed greatly enhanced permeabilities with well-maintained selectivities, a non-linear dependence of permeability on filler content was observed. Another series of MMMs with purified SWNTs (P-SWNTs) and acid-treated SWNT (A-SWNTs) were also prepared. The functionalized A-SWNT MMM exhibited the best separation performance at 2 wt% filler content among all MMMs, which can be attributed to the improved interfacial affinity with carboxylic acid functionalized surface of A-SWNTs.

42 ENGINEERING↗

DESI DR2 Galaxy luminosity functions

We present galaxy luminosity functions (LFs) for the Dark Energy Spectroscopic Instrument (DESI) DR2 Bright Galaxy Survey (BGS) in the g, r, z, and $w1$ bands over $0.002\lt z\lt 0.6$. Our analysis uses updated k-corrections and evolutionary corrections, including new polynomial kcorrection fits derived from BGS Year 1 data that supersede earlier GAMA-based prescriptions. Exploiting the statistical power of DESI, we measure LFs to very faint magnitudes, reaching $^{0.1}M_r-5\log h\sim -10$. Independent measurements from the North and South survey regions show excellent agreement around the LF knee, but the very small statistical uncertainties reveal that simple analytic forms fail to capture the full LF shape. The bright end departs from a pure exponential decline, while the faint end exhibits complex, non-powerlaw behaviour, including a pronounced upturn at $^{0.1}M_r-5\log h\gtrsim -15$, which is stronger for red galaxies than for blue. We show that our LFs are largely complete for galaxies with surface brightness $\mu _{50}\lt 25$, and that an apparent steepening fainter than $-13$ is driven primarily by local overdensity and fragmentation of large galaxies. A systematic North–South offset at the brightest magnitudes is traced to red galaxies and may reflect shallower North photometry underestimating extended earlytype profiles, although this remains inconclusive. We therefore also provide LFs based on model Petrosian magnitudes. Redshift splitting reveals small but significant residuals, indicating limitations of a simple global evolutionary model. Using the redshift limits of J. Loveday et al. (2012), we find excellent agreement with GAMA, with substantially reduced statistical errors. These measurements provide a precise reference for studies of environmental and population-dependent LFs and for testing galaxy formation models.

79 ASTRONOMY AND ASTROPHYSICS↗

Theoretical study of the CO 2 –O 2 van der Waals complex: potential energy surface and applications

A four-dimensional-potential energy surface (4D-PES) of the atmospherically relevant carbon dioxide–oxygen molecule (CO 2 –O 2 ) van der Waals complex is mapped using the ab initio explicitly correlated coupled cluster method with single, double, and perturbative triple excitations (UCCSD(T)-F12b), and extrapolation to the complete basis set (CBS) limit using the cc-pVTZ-F12/cc-pVQZ-F12 bases and the l –3 formula. An analytic representation of the 4D-PES was fitted using the method of interpolating moving least squares (IMLS). These calculations predict that the most stable configuration of CO 2 –O 2 complex corresponds to a planar slipped-parallel structure with a binding energy of V ~ –243 cm –1 . Another isomer is found on the PES, corresponding to a non-planar cross-shaped structure, with V ~ –218 cm –1 . The transition structure connecting the two minima is found at V ~ –211 cm –1 . We also performed comparisons with some CO 2 –X van der Waals complexes. Moreover, we provide a SAPT analysis of this molecular system. Then, we discuss the complexation induced shifts of CO 2 and O 2 . Afterwards, this new 4D-PES is employed to compute the second virial coefficient including temperature dependence. A comparison between quantities obtained in our calculations and those from experiments found close agreement attesting to the high quality of the PES and to the importance of considering a full description of the anisotropic potential for the derivation of thermophysical properties of CO 2 –O 2 mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Density-Functional Tight-Binding Parameters for Bulk Zirconium: A Case Study for Repulsive Potentials

Density-functional tight-binding (DFTB) parameters are presented for the simulation of the bulk phases of zirconium. Electronic parameters were obtained using a band structure fitting strategy, while two-center repulsive potentials were created by particle swarm optimization. As objective functions for the repulsive potential fitting, we employed the Birch–Murnaghan equations of state for hexagonal close-packed (HCP), body-centered cubic (BCC) and ω phases of Zr from density-functional theory (DFT). When fractional atomic coordinates are not allowed to change in the generation of the equation-of-state curves, long-range repulsive DFTB potentials are able to almost perfectly reproduce equilibrium structures, relative DFT energies of the bulk phases, and bulk moduli. However, the same potentials lead to artifacts in the DFTB potential energy surfaces when atom positions in the unit cell are allowed to fully relax during the change of unit cell parameters. Finally, conventional short-range repulsive DFTB potentials, while inferior in their ability to reproduce DFT bulk energetics, are able to correctly reproduce the qualitative shape of the DFT potential energy surfaces, including the location of global minima, and can therefore be considered more transferable.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Weak Gravitational Lensing of Low Surface Brightness Galaxies in the Dark Energy Year 3 Catalog

We present galaxy-galaxy lensing measurements of a sample of low surface brightness galaxies (LSBGs) drawn from the Dark Energy Survey Year 3 (Y3) data. LSBGs are diffuse galaxies with a surface brightness dimmer than the ambient night sky. Given their faintness, the use of standard observational techniques proves challenging. Weak gravitational lensing probes both the baryonic and dark matter content of galaxies, rendering it a powerful technique to estimate LSBG masses. The LSBG lens sample consists of 23, 790 total extended galaxies separated into red and blue color types at \(g-i\ge 0.60\) and \(g-i\le 0.60\), respectively. We use the Y3 \sc{metacalibration} shape catalog as the source sample, with a number density of \(5.59 \mathrm{gal}/\mathrm{arcmin}^2\). We measure the tangential shear around the lens galaxies across angular scales of \(0.25-400\) arcmin and find a signal-to-noise of 6.67 for red galaxies, 2.17 for blue galaxies, and 5.30 for the total sample. We fit a model built from two NFW profiles corresponding to the LSBG dark matter subhalo and host halo to the red LSBG shear measurements with an MCMC. We estimate the host halo mass at \(7.3 ^{+2.0}_{-1.6}\times 10^{12} M_\mathrm{\odot}\). We place a 95\% upper bound on the subhalo mass at \(3\times 10^{11} M_\mathrm{\odot}\). We utilize the lens sample photometry to obtain an estimate of the red LSBG stellar mass distribution. We compare the ratio between the stellar mass and the subhalo mass to the parameterized, satellite-specific SHMR. This work represents the first example of an attempted constraint on the masses of LSBGs using weak gravitational lensing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Weak Gravitational Lensing of Low Surface Brightness Galaxies in the DES Year 3 Catalog

We present galaxy-galaxy lensing measurements of a sample of low surface brightness galaxies (LSBGs) drawn from the Dark Energy Survey Year 3 (Y3) data. LSBGs are diffuse galaxies with a surface brightness dimmer than the ambient night sky. Given their faintness, the use of standard observational techniques proves challenging. Weak gravitational lensing probes both the baryonic and dark matter content of galaxies, rendering it a powerful technique to estimate LSBG masses. The LSBG lens sample consists of 23, 790 total extended galaxies separated into red and blue color types at \(g-i\ge 0.60\) and \(g-i\le 0.60\), respectively. We use the Y3 \metacal\ shape catalog as the source sample, with a number density of {\color{red}\(5.59 \mathrm{gal}/\mathrm{arcmin}^2\)}. We measure the tangential shear around the lens galaxies across angular scales of \(0.25-400\) arcmin and find a signal-to-noise of 6.67 for red galaxies, 2.17 for blue galaxies, and 5.30 for the total sample. W e fit a model built from two NFW profiles corresponding to the LSBG dark matter subhalo and host halo to the red LSBG shear measurements with an MCMC. We estimate the host halo mass at \(7.3 ^{+2.0}_{-1.6}\times 10^{12} M_\mathrm{\odot}\). We place a 95\% upper bound on the subhalo mass at \(2.8\times 10^{11} M_\mathrm{\odot}\). We utilize the lens sample photometry to obtain an estimate of the red LSBG stellar mass distribution. We compare the ratio between the stellar mass and the subhalo mass to the parameterized, satellite-specific SHMR \citep{Moster2010}. This work represents the first example of an attempted constraint on the masses of LSBGs using weak gravitational lensing.

79 ASTRONOMY AND ASTROPHYSICS↗

Disentangling top-down drivers of mortality underlying diel population dynamics of Prochlorococcus in the North Pacific Subtropical Gyre

Photosynthesis fuels primary production at the base of marine food webs. Yet, in many surface ocean ecosystems, diel-driven primary production is tightly coupled to daily loss. This tight coupling raises the question: which top-down drivers predominate in maintaining persistently stable picocyanobacterial populations over longer time scales? Motivated by high-frequency surface water measurements taken in the North Pacific Subtropical Gyre (NPSG), we developed multitrophic models to investigate bottom-up and top-down mechanisms underlying the balanced control of Prochlorococcus populations. Here, we find that incorporating photosynthetic growth with viral- and predator-induced mortality is sufficient to recapitulate daily oscillations of Prochlorococcus abundances with baseline community abundances. In doing so, we infer that grazers in this environment function as the predominant top-down factor despite high standing viral particle densities. The model-data fits also reveal the ecological relevance of light-dependent viral traits and non-canonical factors to cellular loss. Finally, we leverage sensitivity analyses to demonstrate how variation in life history traits across distinct oceanic contexts, including variation in viral adsorption and grazer clearance rates, can transform the quantitative and even qualitative importance of top-down controls in shaping Prochlorococcus population dynamics.

54 ENVIRONMENTAL SCIENCES↗

MaNGA galaxy properties – II. A detailed comparison of observed and simulated spiral galaxy scaling relations

ABSTRACT We present a catalogue of dynamical properties for 2368 late-type galaxies from the MaNGA survey. The latter complements the catalogue of photometric properties for the same sample based on deep optical dark energy sky instrument legacy imaging survey (DESI) photometry processed with autoprof . Rotation curves (RCs), extracted by model-fitting H α velocity maps from the MaNGA Data Analysis Pipeline, extend out to 1.4 (1.9) Re for the primary (secondary) MaNGA samples, respectively. The RCs and ancillary MaNGA Pipe3D data products were used to construct various fundamental galaxy scaling relations that are also compared uniformly with similar relations from numerical investigation of a hundred astrophysical object (NIHAO) zoom-in simulations. Simulated NIHAO galaxies were found to broadly reproduce the observed MaNGA galaxy population for $\log (M_*/{\rm {\rm M}_{\odot })\, \gt\, 8.5}$. Some discrepancies remain, such as those pertaining to central stellar densities and the diversity of RCs due to strong feedback schemes. Also presented are spatially resolved scatters for the velocity-size-stellar mass (VRM*) structural relations using MaNGA and NIHAO samples. The scatter for these relations in the galaxian interiors is a consequence of the diversity of inner RC shapes, while scatter in the outskirts is dictated by the large range of stellar surface densities, which is itself driven by sporadic star formation. The detailed, spatially resolved scatter analysis highlights the complex interplay between local and global astrophysical processes and provides a strong constraint to numerical simulations.

Astronomy & Astrophysics↗

Combining PCA and nonlinear fitting of peak models to re-evaluate C 1s XPS spectrum of cellulose

Cellulose is an example of a material that responds to XPS by the creation of new chemistry not present in the as-received sample. While improvements in instrumentation may be seen in general as beneficial to surface science, recent studies have shown that the consequences for some materials are detrimental. Here, in this work, these problems are illustrated through an analysis of cellulose spectra obtained during a degradation study. C 1s spectra are decomposed into two well-formed component curves that are open to chemical interpretation. In particular, a component-curve representative of pure cellulose is obtained as well as a second component curve that implies cellulose is degraded through the creation of carbon chemistry involving C—O, C$=$O and O—C$=$O. Since cellulose is a crystalline material, formed through the alignment of molecules under the influence of hydrogen bonds, the analysis and findings presented in this paper are relevant to any material analyzed by XPS whose properties are dependent on hydrogen bonds. The analysis techniques are based on an informed vectorial approach, which extracts directly from data spectral shapes that are used to monitor sample degradation via linear least squares optimization. Related mathematics of Principal Component Analysis and linear analysis are presented.

36 MATERIALS SCIENCE↗

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗

Wavelength modulation laser-induced fluorescence for plasma characterization

Laser-Induced Fluorescence (LIF) spectroscopy is an essential tool for probing ion and atom velocity distribution functions (VDFs) in complex plasmas. VDFs carry information about the kinetic properties of species that is critical for plasma characterization. Accurate interpretation of these functions is challenging due to factors such as multicomponent distributions, broadening effects, and background emissions. Our research investigates the use of Wavelength Modulation (WM) LIF to enhance the sensitivity of VDF measurements. Unlike standard Amplitude Modulation (AM) methods, WM–LIF measures the derivative of the LIF signal. This approach makes variations in VDF shape more pronounced. VDF measurements with WM–LIF were investigated with both numerical modeling and experimental measurements. The developed model enables the generation of both WM and AM signals, facilitating comparative analysis of fitting outcomes. Experiments were conducted in a weakly collisional argon plasma with magnetized electrons and non-magnetized ions. Measurements of the argon ion VDFs employed a narrow-band tunable diode laser, which scanned the 4p 4 D 7/2 –3d 4 F 9/2 transition centered at 664.553 nm in vacuum. A lock-in amplifier detected the second harmonic WM signal, which was generated by modulating the laser wavelength with an externally controlled piezo-driven mirror of the diode laser. Finally, our findings indicate that the WM–LIF signal is more sensitive to fitting parameters, allowing for better identification of VDF parameters such as the number of distribution components, their temperatures, and velocities. In addition, WM–LIF can serve as an independent method to verify AM measurements and is particularly beneficial in environments with substantial light noise or background emissions, such as those involving thermionic cathodes and reflective surfaces.

47 OTHER INSTRUMENTATION↗

Refinements to the Graves–Pitarka Kinematic Rupture Generator, Including a Dynamically Consistent Slip-Rate Function, Applied to the 2019 Mw 7.1 Ridgecrest Earthquake

The main objective of this study is to develop physics–based constraints on the spatiotemporal variation of the slip–rate function using a simplified dynamic rupture model. First, we performed dynamic rupture modeling of the 2019 M w 7.1 Ridgecrest, California, earthquake, to analyze the effects of depth–dependent stress and material friction on slip rate. Then, we used our modeling results to guide refinements to the slip–rate function that were implemented in the Graves–Pitarka kinematic rupture generation technique. The dynamic ruptures were computed on a surface–rupturing, planar strike–slip fault that includes a weak (negative to low–stress–drop) zone in the upper 4 km of the crust. Below the weak zone, we placed high–stress–drop patches designed to mirror the large–slip areas seen in various rupture model inversions of the event. The locations of the high–stress–drop patches and the hypocenter were varied in multiple realizations to investigate how changing the dynamic conditions affected the resulting rupture kinematics, in particular, the slip rate. From these simulations, we observed a systematic change in the shape of the slip–rate function from Kostrov type below the weak zone to a predominantly symmetric shape within the weak zone, along with a depth–dependent reduction of peak slip rate. We generalized these shallow rupture features into a depth–dependent parametric variation of the slip–rate function and implemented it in the Graves–Pitarka kinematic rupture model generator. The performance of the updated kinematic approach was then verified in 0–4 Hz simulations of the M w 7.1 Ridgecrest earthquake, which showed that incorporating the depth–dependent variation in the shape of the slip–rate function improves the fit to the observed near–fault ground motions in the 0.5–3 s period range.

58 GEOSCIENCES↗

A White Dwarf with Transiting Circumstellar Material Far outside the Roche Limit

We report the discovery of a white dwarf exhibiting deep, irregularly shaped transits, indicative of circumstellar planetary debris. Using Zwicky Transient Facility DR2 photometry of ZTF J013906.17+524536.89 and follow-up observations from the Las Cumbres Observatory, we identify multiple transit events that recur every ≈107.2 days, much longer than the 4.5–4.9 hr orbital periods observed in WD 1145+017, the only other white dwarf known with transiting planetary debris. The transits vary in both depth and duration, lasting 15–25 days and reaching 20%–45% dips in flux. Optical spectra reveal strong Balmer lines, identifying the white dwarf as a DA with ${T}_{\mathrm{eff}}=10,530\pm 140\,{\rm{K}}$ and $\mathrm{log}(g)=7.86\pm 0.06$. A Ca ii K absorption feature is present in all spectra both in and out of transit. Spectra obtained during one night at roughly 15% transit depth show increased Ca ii K absorption with a model atmospheric fit suggesting [Ca/H] = –4.6 ± 0.3, whereas spectra taken on three nights out of transit have [Ca/H] of –5.5, –5.3, and –4.9 with similar uncertainties. While the Ca ii K line strength varies by only 2σ, we consider a predominantly interstellar origin for Ca absorption unlikely. Here, we suggest a larger column density of circumstellar metallic gas along the line of site or increased accretion of material onto the white dwarf's surface are responsible for the Ca absorption, but further spectroscopic studies are required. In addition, high-speed time series photometry out of transit reveals variability with periods of 900 and 1030 s, consistent with ZZ Ceti pulsations.

79 ASTRONOMY AND ASTROPHYSICS↗

First-Principles Treatment of Vibrational Broadening in X-Ray Excited Valence Band Spectra for n-SrTiO 3 (001)

The valence band maximum (VBM) is an important quantity for semiconductors as it locates the Fermi level relative to the band edge. Accurate measurement of this quantity in near-surface regions of semiconductors by photoemission is a first step toward determining the electronic properties of heterostructures involving these materials. While extrapolating the leading edge of the valence band to the energy axis in photoemission spectra is a widely used way to find the VBM, this method can be ambiguous if the leading edge exhibits multiple slopes. Another way to determine the VBM is to fit the leading edge to an appropriately broadened, cross-section modulated theoretical density of states (DOS). Three kinds of broadening that should be included for maximum accuracy are those due to: (1) finite instrumental resolution, (2) valence hole lifetime, and (3) vibrational excitations. While steps (1) and (2) are straightforward to implement, (3) is more difficult because the appropriate amount of broadening is not known a priori. Here, we demonstrate that explicit inclusion of vibrational broadening using ab initio molecular dynamics facilitates accurate VBM determination for n-SrTiO 3 (001). The total DOS is constructed by summing time-averaged projections at elevated temperature onto s - , p - , and d orbitals for the constituent atoms and modulating with the associated photoemission cross sections. Subsequent convolutions of the total DOS, first with a Gaussian of width equal to the experimental energy resolution and second with a Lorentzian to simulate valence hole lifetime effects, yield line shapes that reproduce the experimental leading edges rather well. The VBM is then given by the energy at which the vibrationally broadened total DOS (prior to the convolutions) goes to zero. The VBMs generated by this method quantitatively agree with those resulting from extrapolating from the middle of the measured leading edge for SrTiO 3 .

36 MATERIALS SCIENCE↗