Search NASA⌕ Search

SEARCH · Search NASA

Results for “Distance Based Approximation”

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 217 records · Page 12

Enhanced Deuteron Coalescence Probability in Jets

The transverse-momentum ($p_T$) spectra and coalescence parameters $B_2$ of (anti)deuterons are measured in p-p collisions at $\sqrt{s}=13$ $\text{TeV}$ for the first time in and out of jets. In this measurement, the direction of the leading particle with the highest $p_T$ in the event ($p^{lead}_T > 5 GeV/c$) is used as an approximation for the jet axis. The event is consequently divided into three azimuthal regions, and the jet signal is obtained as the difference between the toward region, that contains jet fragmentation products in addition to the underlying event (UE), and the transverse region, which is dominated by the UE. The coalescence parameter in the jet is found to be approximately a factor of 10 larger than that in the underlying event. This experimental observation is consistent with the coalescence picture and can be attributed to the smaller average phase-space distance between nucleons in the jet cone as compared with the underlying event. The results presented in this Letter are compared to predictions from a simple nucleon coalescence model, where the phase-space distributions of nucleons are generated using pythia8 with the Monash 2013 tuning, and to predictions from a deuteron production model based on ordinary nuclear reactions with parametrized energy-dependent cross sections tuned on data. The latter model is implemented in pythia8.3. Both models reproduce the observed large difference between in-jet and out-of-jet coalescence parameters, although the almost flat trend of the $B^{Jet}_2$ is not reproduced by the models, which instead give a decreasing trend.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

eFESTA: Ensemble Feature Exploration with Surface Density Estimates

We propose surface density estimate (SDE) to model the spatial distribution of surface features-isosurfaces, ridge surfaces, and streamsurfaces-in 3D ensemble simulation data. The inputs of SDE computation are surface features represented as polygon meshes, and no field datasets are required (e.g., scalar fields or vector fields). The SDE is defined as the kernel density estimate of the infinite set of points on the input surfaces and is approximated by accumulating the surface densities of triangular patches. We also propose an algorithm to guide the selection of a proper kernel bandwidth for SDE computation. An ensemble Feature Exploration method based on Surface densiTy EstimAtes (eFESTA) is then proposed to extract and visualize the major trends of ensemble surface features. For an ensemble of surface features, each surface is first transformed into a density field based on its contribution to the SDE, and the resulting density fields are organized into a hierarchical representation based on the pairwise distances between them. The hierarchical representation is then used to guide visual exploration of the density fields as well as the underlying surface features. We demonstrate the application of our method using isosurface in ensemble scalar fields, Lagrangian coherent structures in uncertain unsteady flows, and streamsurfaces in ensemble fluid flows.

Computational modeling↗

Deformable image registration based automatic CT-to-CT contour propagation for head and neck adaptive radiotherapy in the routine clinical setting

Purpose: To evaluate the clinical potential of deformable image registration (DIR)-based automatic propagation of physician-drawn contours from a planning CT to midtreatment CT images for head and neck (H and N) adaptive radiotherapy. Methods: Ten H and N patients, each with a planning CT (CT1) and a subsequent CT (CT2) taken approximately 3–4 week into treatment, were considered retrospectively. Clinically relevant organs and targets were manually delineated by a radiation oncologist on both sets of images. Four commercial DIR algorithms, two B-spline-based and two Demons-based, were used to deform CT1 and the relevant contour sets onto corresponding CT2 images. Agreement of the propagated contours with manually drawn contours on CT2 was visually rated by four radiation oncologists in a scale from 1 to 5, the volume overlap was quantified using Dice coefficients, and a distance analysis was done using center of mass (CoM) displacements and Hausdorff distances (HDs). Performance of these four commercial algorithms was validated using a parameter-optimized Elastix DIR algorithm. Results: All algorithms attained Dice coefficients of >0.85 for organs with clear boundaries and those with volumes >9 cm{sup 3}. Organs with volumes <3 cm{sup 3} and/or those with poorly defined boundaries showed Dice coefficients of ∼0.5–0.6. For the propagation of small organs (<3 cm{sup 3}), the B-spline-based algorithms showed higher mean Dice values (Dice = 0.60) than the Demons-based algorithms (Dice = 0.54). For the gross and planning target volumes, the respective mean Dice coefficients were 0.8 and 0.9. There was no statistically significant difference in the Dice coefficients, CoM, or HD among investigated DIR algorithms. The mean radiation oncologist visual scores of the four algorithms ranged from 3.2 to 3.8, which indicated that the quality of transferred contours was “clinically acceptable with minor modification or major modification in a small number of contours.” Conclusions: Use of DIR-based contour propagation in the routine clinical setting is expected to increase the efficiency of H and N replanning, reducing the amount of time needed for manual target and organ delineations.

60 APPLIED LIFE SCIENCES↗

Decay Length Estimation of Single- Two- and Three-Wire Systems Above Ground under HEMP Excitation

Here, we analytically model single-, two-, and three-wires above ground to determine the decay lengths of common and differential modes induced by an E1 high-altitude electromagnetic pulse (HEMP) excitation. Decay length information is pivotal to determine whether any two nodes in the power grid may be treated as uncoupled. We employ a frequency-domain method based on transmission line theory named ATLOG - Analytic Transmission Line Over Ground to model infinitely long and finite single wires, as well as solve the eigenvalue problem of a single-,two-, and three-wire system. Our calculations show that a single, semi-infinite power line can be approximated by a10 km section of line and that the second electrical reflection for all line lengths longer than the decay length are below half the rated operating voltage. Moreover, our findings show that the differential mode propagates longer distances than the common mode in two- and three-wire systems, and this should be taken into account when performing damage assessment from HEMP excitation. Furthermore, this analysis is a significant step toward simplifying the modeling of practical continental grid lengths, yet maintaining accuracy, a result of enormous impact.

42 ENGINEERING↗

Robust Spectral Unmixing of Sparse Multispectral Lidar Waveforms using Gamma Markov Random Fields

Here, this paper presents a new Bayesian spectral un-mixing algorithm to analyse remote scenes sensed via sparse multispectral Lidar measurements. To a first approximation, in the presence of a target, each Lidar waveform consists of a main peak, whose position depends on the target distance and whose amplitude depends on the wavelength of the laser source considered (i.e, on the target reflectivity). Besides, these temporal responses are usually assumed to be corrupted by Poisson noise in the low photon count regime. When considering multiple wavelengths, it becomes possible to use spectral information in order to identify and quantify the main materials in the scene, in addition to estimation of the Lidar-based range profiles. Due to its anomaly detection capability, the proposed hierarchical Bayesian model, coupled with an efficient Markov chain Monte Carlo algorithm, allows robust estimation of depth images together with abundance and outlier maps associated with the observed 3D scene. The proposed methodology is illustrated via experiments conducted with real multispectral Lidar data acquired in a controlled environment. The results demonstrate the possibility to unmix spectral responses constructed from extremely sparse photon counts (less than 10 photons per pixel and band).

97 MATHEMATICS AND COMPUTING↗

How to Turn Jets into Cylinders near Supermassive Black Holes in 3D General Relativistic Magnetohydrodynamic Simulations

Abstract Accreting supermassive black holes (SMBHs) produce highly magnetized relativistic jets that tend to collimate gradually as they propagate outward. However, recent radio interferometric observations of the 3C 84 galaxy reveal a stunning, cylindrical jet already at several hundred SMBH gravitational radii,r≳ 350r g . We explore how such extreme collimation emerges via a suite of 3D general relativistic magnetohydrodynamic simulations. We consider an SMBH surrounded by a magnetized torus immersed in a constant-density ambient medium that starts at the edge of the SMBH sphere of influence, chosen to be much larger than the SMBH gravitational radius,r B = 10 3 r g . We find that radiatively inefficient accretion flows (e.g., M87) produce winds that collimate the jets into parabolas near the black hole. After the disk winds stop collimating the jets atr≲r B , they turn conical. Once outsider B , the jets run into the ambient medium and form backflows that collimate the jets into cylinders some distance beyondr B . Interestingly, for radiatively efficient accretion, as in 3C 84, the radiative cooling saps the energy out of the disk winds; at early times, they cannot efficiently collimate the jets, which skip the initial parabolic collimation stage, start out conical near the SMBH, and turn into cylinders already atr≃ 300r g , as observed in 3C 84. Over time, the jet power remains approximately constant, whereas the mass accretion rate increases; the winds grow in strength and start to collimate the jets, which become quasi-parabolic near the base, and the transition point to a nearly cylindrical jet profile moves outward while remaining insider B .

Astronomy & Astrophysics↗

Hexagonal Distributed Embedded Energy Converters (HexDEECs)

Distributed Embedded Energy Converter Technologies (DEEC-Tec) is a new domain for marine renewable energy research that utilizes a conglomeration of small distributed embedded energy converters (DEECs) that, in aggregate, form larger metamaterial frameworks. These resulting DEEC-Tec metamaterials can then, in turn, be used to construct flexible ocean wave energy converters called flexWECs. DEEC-Tec enables flexWECs: (i) to be inherently broad-banded ocean wave frequency energy converters and (ii) to have an inherent lack of highly loaded rigid bodies. The DEEC-Tec domain also benefits the marine renewable energy domain by inherently availing ways that marine energy can be harvested and converted that heretofore has not yet been considered possible: real-time execution of transforming topologies (e.g., actively changing a flexWEC's shape and form) and morphologies (e.g., actively changing a flexWEC's stiffness and damping throughout its entire structure). Presented, is one specific type of DEEC, a HexDEEC, that shows promise in aiding the adoption and further development of the DEEC-Tec domain - it is a small energy transducer being developed by the United States National Renewable Energy Laboratory. The HexDEEC is a small (characteristic length approximating a centimeter) energy transducer that converts the dynamic deformations of an elastomer into electricity through a charging-discharging cycle of a capacitor whose capacitance is varied by those elastic deformations. The HexDEEC is composed of a hyperelastic hexagonal housing (nominally silicon rubber) with six electrodes on its inner faces. The upper three electrodes share the same charge while the lower three electrodes oppose the upper electrode charges. Externally, the HexDEEC has two arms extending away from the middle vertices of the hexagon. Via principles governing the relationship between electrical capacitance and electrical potential (voltage and charge), electricity is generated when the HexDEEC's arms are dynamically pulled or released under tensile loading as doing so causes the distance between the upper and lower sets of electrodes to change - varying the energy converter's overall capacitance. Analytical and numerical modeling have already been used to estimate the electrical energy produced by a HexDEEC. The cursory models approximate the HexDEEC as a parallel plate variable capacitor - simplifying from six to two opposing plates with a constant dielectric volume between those two plates. To account for the elastic HexDEEC material properties, software such as SolidWorks and STAR-CCM+ have been used to generate hyperelastic models; notably, Mooney-Rivlin based models. Individual HexDEECs have been fabricated by drawing uncured liquid silicon rubber into molds via vacuum pressure. To simplify manufacturing, HexDEEC sub-components - e.g., electrodes, wires - can be placed within those molds such that they are directly embedded into the hexagonal housing during the curing process. Furthermore, DEEC-Tec metamaterials can be created by interweaving or sequentially layering multiple HexDEEC strands together. The HexDEEC based metamaterial could then generate electricity through its gross deformations. Ultimately, HexDEECs represent a specific type of energy transducer that can be leveraged, by the DEEC-Tec domain, to create metamaterials used to construct novel flexWECs.

DEEC-Tec↗

ON SUN-TO-EARTH PROPAGATION OF CORONAL MASS EJECTIONS

We investigate how coronal mass ejections (CMEs) propagate through, and interact with, the inner heliosphere between the Sun and Earth, a key question in CME research and space weather forecasting. CME Sun-to-Earth kinematics are constrained by combining wide-angle heliospheric imaging observations, interplanetary radio type II bursts, and in situ measurements from multiple vantage points. We select three events for this study, the 2012 January 19, 23, and March 7 CMEs. Different from previous event studies, this work attempts to create a general picture for CME Sun-to-Earth propagation and compare different techniques for determining CME interplanetary kinematics. Key results are obtained concerning CME Sun-to-Earth propagation: (1) the Sun-to-Earth propagation of fast CMEs can be approximately formulated into three phases: an impulsive acceleration, then a rapid deceleration, and finally a nearly constant speed propagation (or gradual deceleration); (2) the CMEs studied here are still accelerating even after the flare maximum, so energy must be continuously fed into the CME even after the time of the maximum heating and radiation has elapsed in the corona; (3) the rapid deceleration, presumably due to interactions with the ambient medium, mainly occurs over a relatively short timescale following the acceleration phase; and (4) CME-CME interactions seem a common phenomenon close to solar maximum. Our comparison between different techniques (and data sets) has important implications for CME observations and their interpretations: (1) for the current cases, triangulation assuming a compact CME geometry is more reliable than triangulation assuming a spherical front attached to the Sun for distances below 50-70 solar radii from the Sun, but beyond about 100 solar radii we would trust the latter more; (2) a proper treatment of CME geometry must be performed in determining CME Sun-to-Earth kinematics, especially when the CME propagation direction is far away from the observer; and (3) our approach to comparing wide-angle heliospheric imaging observations with interplanetary radio type II bursts provides a novel tool in investigating CME propagation characteristics. Future CME observations and space weather forecasting are discussed based on these results.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigation of the Effect of Framework Flexibility on CO 2 Adsorption in SIFSIX-3-Cu Using a Machine-Learned Force Field

Metal–organic frameworks (MOFs) offer promise as selective CO 2 sorbents, but successful MOF sorbent materials need high CO 2 binding affinity and selectivity for CO 2 over water. This work focuses on the use of machine-learned force fields (MLFFs) to model CO 2 adsorption in flexible MOFs, with a focus on SIFSIX-3-Cu, an anion-pillared MOF known for its high CO 2 affinity. A preliminary high-throughput screening of over 900 anion-pillared MOFs was performed using rigid UFF+DDEC6 force fields to predict zero-loading heats of adsorption for CO 2 and H 2 O. SIFSIX-3-Cu was selected for further computational study due to its predicted CO 2 heat of adsorption and experimental relevance. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies and forces, with an iterative sampling scheme combining molecular dynamics, geometry optimization, random geometric insertion, and NVT Monte Carlo-based configuration generation to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included, contrasting with previous models that approximated the MOF as rigid. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations with the MLFF produced CO 2 adsorption isotherms in good agreement with experimental data at direct air capture (DAC) pressures (e.g., 40 Pa), in contrast to previous overestimations of CO 2 sorption by models with rigid structures. Bond and angle histogram analysis showed that MOF flexibility increased the variance of fluorine–fluorine diagonal distances at adsorption sites, resulting in a lower predicted sorption for flexible, asymmetric SIFSIX-3-Cu pore geometries compared to the rigid, symmetric DFT-optimized SIFSIX-3-Cu pore geometry. A detailed description of flexibility afforded by the MLFF resulted in an accurately predicted CO 2 uptake (0.88 mmol/g) at low pressure (40 Pa) compared to the experimentally measured value (1.24 mmol/g). In conclusion, these results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low-pressure applications.

adsorption↗

Critical Height of the Torus Instability in Two-ribbon Solar Flares

We studied the background field for 60 two-ribbon flares of M-and-above classes during 2011–2015. These flares are categorized into two groups, i.e., eruptive and confined flares, based on whether a flare is associated with a coronal mass ejection or not. The background field of source active regions is approximated by a potential field extrapolated from the B {sub z} component of vector magnetograms provided by the Helioseismic and Magnetic Imager. We calculated the decay index n of the background field above the flaring polarity inversion line, and defined a critical height h {sub crit} corresponding to the theoretical threshold ( n {sub crit} = 1.5) of the torus instability. We found that h {sub crit} is approximately half of the distance between the centroids of opposite polarities in active regions and that the distribution of h {sub crit} is bimodal: it is significantly higher for confined flares than for eruptive ones. The decay index increases monotonously with increasing height for 86% (84%) of the eruptive (confined) flares but displays a saddle-like profile for the rest, 14% (16%), which are found exclusively in active regions of multipolar field configuration. Moreover, n at the saddle bottom is significantly smaller in confined flares than that in eruptive ones. These results highlight the critical role of background field in regulating the eruptive behavior of two-ribbon flares.

79 ASTRONOMY AND ASTROPHYSICS↗

Feasibility of Multimodal Deformable Registration for Head and Neck Tumor Treatment Planning

Purpose: To investigate the feasibility of using deformable registration in clinical practice to fuse MR and CT images of the head and neck for treatment planning. Method and Materials: A state-of-the-art deformable registration algorithm was optimized, evaluated, and compared with rigid registration. The evaluation was based on manually annotated anatomic landmarks and regions of interest in both modalities. We also developed a multiparametric registration approach, which simultaneously aligns T1- and T2-weighted MR sequences to CT. This was evaluated and compared with single-parametric approaches. Results: Our results show that deformable registration yielded a better accuracy than rigid registration, without introducing unrealistic deformations. For deformable registration, an average landmark alignment of approximatively 1.7 mm was obtained. For all the regions of interest excluding the cerebellum and the parotids, deformable registration provided a median modified Hausdorff distance of approximatively 1 mm. Similar accuracies were obtained for the single-parameter and multiparameter approaches. Conclusions: This study demonstrates that deformable registration of head-and-neck CT and MR images is feasible, with overall a significanlty higher accuracy than for rigid registration.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Electronic structure, pore size distribution, and sorption characterization of an unusual MOF, {[Ni(dpbz)][Ni(CN) 4 ]} n , dpbz = 1,4-bis(4-pyridyl)benzene

The monoclinic (Ni(L)[Ni(CN) 4 ] (L= 1,4-Bis(4-pyridyl) benzene) compound (defined as Ni-dpbz) is a flexible metal organic framework which assumes a pillared structure with layers defined by 2D Ni[Ni(CN) 4 ] n nets and dpbz ligands as pillars. The structure features an entrapped dpbz ligand that links between the open ends of four-fold Ni sites from two neighboring chains. This arrangement results in an unusual 5-fold pseudo square-pyramid environment for Ni and a significantly long Ni-N distance of 2.369(4) Å. Herein, using Density Functional Theory calculations, the different bonding characteristics between the 5-fold and 6-fold Ni's were determined. We found that there is weak covalent bonding between the 5-fold Ni and N in the entrapped ligand, and the 6-fold Ni-N bonds provide effective electronic conduction. The disordered dimethyl sulfoxide (DMSO) solvent molecules are not bonded to the framework. The material has a single pore with a diameter of 4.1 Å. This pore includes approximately 55% of the total free volume (based on a zero-diameter probe). The accessible pore surface area and pore volume were calculated to be 507 m 2 /g and 6.99 cm 3 /kg, respectively. The maximum amount of CO 2 that can be accommodated in the pores after DMSO is removed was found to be 204 mg/g, agreeing with the results of adsorption/desorption experiments of about 220 mg/g.

36 MATERIALS SCIENCE↗

Chiral effective field theory calculations of weak transitions in light nuclei

In this work, we report quantum Monte Carlo calculations of weak transitions in A≤10 nuclei, based on the Norfolk two- and three-nucleon chiral interactions, and associated one- and two-body axial currents. Furthermore, we find that the contribution from two-body currents is at the 2–3% level, with the exception of matrix elements entering the rates of 8Li, 8B, and 8He β decay. These matrix elements are suppressed in impulse approximation based on the (leading order) Gamow Teller transition operator alone; two-body currents provide a 20–30% correction, which is, however, insufficient to bring theory in agreement with experimental data. For the other transitions, the agreement with the data is satisfactory, and the results exhibit a negligible to mild model dependence when different combinations of Norfolk interactions are utilized to construct the nuclear wave functions. We report a complete study of two-body weak transition densities which reveals the expected universal behavior of two-body currents at short distances throughout the range of A=3 to A=10 systems considered here.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes

The present contribution focuses on the accuracy of reflection-type boundary conditions in the Stokes–Brinkman–Darcy modeling of porous flows solved with the lattice Boltzmann method (LBM), which we operate with the two-relaxation-time (TRT) collision and the Brinkman-force based scheme (BF), called BF-TRT scheme. In parallel, we compare it with the Stokes–Brinkman–Darcy linear finite element method (FEM) where the Dirichlet boundary conditions are enforced on grid vertices. In bulk, both BF-TRT and FEM share the same defect: in their discretization a correction to the modeled Brinkman equation appears, given by the discrete Laplacian of the velocity-proportional resistance force. This correction modifies the effective Brinkman viscosity, playing a crucial role in the triggering of spurious oscillations in the bulk solution. While the exact form of this defect is available in lattice-aligned, straight or diagonal, flows; in arbitrary flow/lattice orientations its approximation is constructed. At boundaries, we verify that such a Brinkman viscosity correction has an even more harmful impact. Already at the first order, it shifts the location of the no-slip wall condition supported by traditional LBM boundary schemes, such as the bounce-back rule. For that reason, this work develops a new class of boundary schemes to prescribe the Dirichlet velocity condition at an arbitrary wall/boundary-node distance and that supports a higher order accuracy in the accommodation of the TRT-Brinkman solutions. For their modeling, we consider the standard BF scheme and its improved version, called IBF; this latter is generalized in this work to suppress or to reduce the viscosity correction in arbitrarily oriented flows. Our framework extends the one- and two-point families of linear and parabolic link-wise boundary schemes, respectively called B-LI and B-MLI, which avoid the interference of the Brinkman viscosity correction in their closure relations. The performance of LBM and FEM is thoroughly evaluated in three benchmark tests, which are run throughout three distinctive permeability regimes. The first configuration is a horizontal porous channel, studied with a symbolic approach, where we construct the exact solutions of FEM and BF/IBF with different boundary schemes. The second problem refers to an inclined porous channel flow, which brings in as new challenge the formation of spurious boundary layers in LBM; that is, numerical artefacts that arise due to a deficient accommodation of the bulk solution by the low-accurate boundary scheme. The third problem considers a porous flow past a periodic square array of solid cylinders, which intensifies the previous two tests with the simulation of a more complex flow pattern. The ensemble of numerical tests provides guidelines on the effect of grid resolution and the TRT free collision parameter over the accuracy and the quality of the velocity field, spanning from Stokes to Darcy permeability regimes. It is shown that, with the use of the high-order accurate boundary schemes, the simple, uniform-mesh-based TRT-LBM formulation can even surpass the accuracy of FEM employing hardworking body-fitted meshes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On the use of NMR distance measurements for assessing surface site homogeneity

The past few decades have seen tremendous growth in the area of single-site heterogeneous catalysis, which aims to combine the best aspects of homogeneous and heterogeneous catalysis, namely molecular-level site control and ease of separation/recycling. Despite this, we still do not have a means of assessing site homogeneity and whether the produced catalyst is indeed a “single-site”. Recent developments have enabled the use of NMR-based distance measurements to determine the conformations and configurations of surface sites, leading to the question whether such measurements can be used to distinguish materials containing either single or multiple surface sites with otherwise indistinguishable NMR properties. Here, we describe a Monte Carlo-based multi-structure search algorithm and its application to the determination of multi-site structures from supported metal complexes. The sensitivity of REDOR data to the existence of multiple sites is assessed using synthetic data and prior literature examples are revisited to determine whether the single-site approximation was indeed appropriate. We lastly apply this new methodology to differentiate the configurations of zirconocene complexes grafted onto alumina supports that were thermally treated at different temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermal imaging of high power diode lasers subject to back-irradiance

In this study, CCD-based thermoreflectance imaging and finite element modeling are used to study the two-dimensional (2D) temperature profile of a junction-down broad-area diode laser facet subject to back-irradiance. By determining the temperature rise in the active region (ΔΤAR) at different diode laser optical powers, back-irradiance reflectance levels, and back-irradiance spot locations, we find that ΔΤAR increases by nearly a factor of three when the back-irradiance spot is centered in the absorbing substrate approximately 5 μm away from the active region, a distance roughly equal to half of the back-irradiance spot FWHM (9 μm). This corroborates prior work studying the relationship between the back-irradiance spot location and catastrophic optical damage, suggesting a strong thermal basis for reduced laser lifetime in the presence of back-irradiance for diode lasers fabricated on absorbing substrates.

42 ENGINEERING↗

Full Law Analysis Scattering System Hub (FLASSH)

In neutronic systems such as thermal reactors, the thermal neutron plays a fundamental role in the system behavior. Therefore, high fidelity thermal neutron scattering cross section evaluations are required in the design of such systems. Because the de Broglie wavelength of thermal neutrons is comparable to the atomic separation distance in matter, the neutron scattering behavior in the energy region of atomic thermal motion is fundamentally different from high energy regions. In recent years, algorithms and theories of thermal neutron scattering have been developed and verified. Based on that effort, a capability for full analysis of the scattering law in an atomistic system (FLASSH) has been developed. FLASSH not only reproduces the evaluations of conventional thermal neutron scattering codes, such as the LEAPR and THERMR modules of the NJOY code, but also improves the accuracy of evaluations by removing or relaxing numerous approximations and by implementing advanced algorithm. In this work, α-SiO{sub 2} and crystalline beryllium metal are taken as examples to demonstrate the evaluation of the inelastic thermal scattering cross section and the coherent elastic scattering cross section using FLASSH.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Excitonic couplings between molecular crystal pairs by a multistate approximation

In this paper, we present a diabatization scheme to compute the excitonic couplings between an arbitrary number of states in molecular pairs. The method is based on an algebraic procedure to find the diabatic states with a desired property as close as possible to that of some reference states. In common with other diabatization schemes, this method captures the physics of the important short-range contributions (exchange, overlap, and charge-transfer mediated terms) but it becomes particularly suitable in presence of more than two states of interest. The method is formulated to be usable with any level of electronic structure calculations and to diabatize different types of states by selecting different molecular properties. These features make the diabatization scheme presented here especially appropriate in the context of organic crystals, where several excitons localized on the same molecular pair may be found close in energy. In this paper, the method is validated on the tetracene crystal dimer, a well characterized case where the charge transfer (CT) states are closer in energy to the Frenkel excitons (FE). The test system was studied as a function of an external electric field (to explore the effect of changing the relative energy of the CT excited state) and as a function of different intermolecular distances (to probe the strength of the coupling between FE and CT states). Additionally, we illustrate how the approximation can be used to include the environment polarization effect.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗