Search NASA⌕ Search

SEARCH · Search NASA

Results for “Function Estimation”

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 703 records · Page 39

Local Stability and Global Instability in Iron-opaque Disks

The thermal stability of accretion disks and the possibility of seeing a limit-cycle behavior strongly depends on the ability of the disk plasma to cool down. Various processes connected with radiation–matter interaction appearing in hot accretion disk plasma contribute to opacity. For the case of geometrically thin and optically thick accretion disks, we can estimate the influence of several different components of function κ, given by the Roseland mean. In the case of high temperatures of ∼10{sup 7} K, the electron Thomson scattering is dominant. At lower temperatures, atomic processes become important. The slope dlogκ/dlogT can have a locally stabilizing or destabilizing effect on the disk. Although the local MHD simulation postulates the stabilizing influence of the atomic processes, only the global time-dependent model can reveal the global disk stability range estimation. This is due to the global diffusive nature of those processes. In this paper, using the previously tested GLADIS code with a modified prescription of the viscous dissipation, we examine the stabilizing effect of the iron opacity bump.

79 ASTRONOMY AND ASTROPHYSICS↗

Assessment of electronic structure methods for the determination of the ground spin states of Fe( ii ), Fe( iii ) and Fe( iv ) complexes

Our ability to understand and simulate the reactions catalyzed by iron depends strongly on our ability to predict the relative energetics of spin states. In this work, we studied the electronic structures of Fe 2+ ion, gaseous FeO and 14 iron complexes using Kohn–Sham density functional theory with particular focus on determining the ground spin state of these species as well as the magnitudes of relevant spin-state energy splittings. The 14 iron complexes investigated in this work have hexacoordinate geometries of which seven are Fe(II), five are Fe(III) and two are Fe(IV) complexes. These are calculated using 20 exchange–correlation functionals. In particular, we use a local spin density approximation (LSDA) – GVWN5, four generalized gradient approximations (GGAs) – BLYP, PBE, OPBE and OLYP, two non-separable gradient approximations (NGAs) – GAM and N12, two meta-GGAs – M06-L and M11-L, a meta-NGA – MN15-L, five hybrid GGAs – B3LYP, B3LYP*, PBE0, B97-3 and SOGGA11-X, four hybrid meta-GGAs – M06, PW6B95, MPW1B95 and M08-SO and a hybrid meta-NGA – MN15. The density functional results are compared to reference data, which include experimental results as well as the results of diffusion Monte Carlo (DMC) calculations and ligand field theory estimates from the literature. For the Fe 2+ ion, all functionals except M11-L correctly predict the ground spin state to be quintet. However, quantitatively, most of the functionals are not close to the experimentally determined spin-state splitting energies. For FeO all functionals predict quintet to be the ground spin state. For the 14 iron complexes, the hybrid functionals B3LYP, MPW1B95 and MN15 correctly predict the ground spin state of 13 out of 14 complexes and PW6B95 gets all the 14 complexes right. Here, the local functionals, OPBE, OLYP and M06-L, predict the correct ground spin state for 12 out of 14 complexes. Two of the tested functionals are not recommended to be used for this type of study, in particular M08-SO and M11-L, because M08-SO systematically overstabilizes the high spin state, and M11-L systematically overstabilizes the low spin state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The 2.4 μm Galaxy Luminosity Function As Measured Using WISE. I. Measurement Techniques

The astronomy community has at its disposal a large back catalog of public spectroscopic galaxy redshift surveys that can be used for the measurement of luminosity functions (LFs). Utilizing the back catalog with new photometric surveys to maximum efficiency requires modeling the color selection bias imposed on the selection of target galaxies by flux limits at multiple wavelengths. The likelihood derived herein can address, in principle, all possible color selection biases through the use of a generalization of the LF, Φ(L), over the space of all spectra: the spectro-luminosity functional, Ψ[L{sub ν}]. It is, therefore, the first estimator capable of simultaneously analyzing multiple redshift surveys in a consistent way. We also propose a new way of parametrizing the evolution of the classic Schechter function parameters, L {sub ⋆} and ϕ {sub ⋆}, that improves both the physical realism and statistical performance of the model. The techniques derived in this paper are used in a companion paper by Lake et al. to measure the LF of galaxies at the rest-frame wavelength of 2.4 μm using the Widefield Infrared Survey Explorer (WISE).

79 ASTRONOMY AND ASTROPHYSICS↗

Variational control forces for enhanced sampling of nonequilibrium molecular dynamics simulations

We introduce a variational algorithm to estimate the likelihood of a rare event within a nonequilibrium molecular dynamics simulation through the evaluation of an optimal control force. Optimization of a control force within a chosen basis is made possible by explicit forms for the gradients of a cost function in terms of the susceptibility of driven trajectories to changes in variational parameters. We consider probabilities of time-integrated dynamical observables as characterized by their large deviation functions and find that in many cases, the variational estimate is quantitatively accurate. Additionally, we provide expressions to exactly correct the variational estimate that can be evaluated directly. We benchmark this algorithm against the numerically exact solution of a model of a driven particle in a periodic potential, where the control force can be represented with a complete basis. Additionally, we then demonstrate the utility of the algorithm in a model of repulsive particles on a line, which undergo a dynamical phase transition, resulting in singular changes to the form of the optimal control force. In both systems, we find fast convergence and are able to evaluate large deviation functions with significant increases in statistical efficiency over alternative Monte Carlo approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-Layer Predictive Controller for V2G and G2V Services Using on Wireless Power Transfer Technology

This paper presents a comprehensive solution that allows electric vehicles (EV) to autonomously charge and discharge their batteries wirelessly during long term parking and/or the transient stops. A two-layer power-flow controller for bidirectional wireless power transfer system (BWPTS) in EVs applications is proposed. The proposed controller can manage the bidirectional power-flow between EV and surrounding infrastructures such as power grid, home microgrid, building micro-grid, road or another vehicle. It consists of two levels of control; the first is responsible for communicating with the surrounding infrastructures and gathering information from driver, charging station, power grid and battery management system and then, based on these information, it estimates the EV's psychological price as a function of its battery's state-of-charge (SOC) and compares it with the energy price to decide whether to charge, discharge or abstain, and how much the charging or discharging rate. The second layer receives the reference signal from the first one and generates the control parameters for two synchronized resonant converters (one is on the vehicle side and the other is on the grid side) to provide the requited power-flow. The proposed controller is adaptively estimating the system parameters to consider the misalignment conditions effects on the system performance. The parameter estimation is achieved using only one voltage sensor. The second layer control is designed based on a new analytical modeling for the power flow in the system. For verification purposes, a prototype for BWPTS was built and driven by the proposed controller, which was implemented using a field-programmable gate array (FPGA) integrated circuit. The proposed controller provides very fast and stable response during both the transient and steady state operation in comparison with the conventional PI controller.

33 ADVANCED PROPULSION SYSTEMS↗

Selective electron-phonon coupling strength from nonequilibrium optical spectroscopy: The case of MgB 2

The coupling between quasiparticles and bosonic excitations rules the energy transfer pathways in condensed matter systems. The possibility of inferring the strength of specific coupling channels from their characteristic timescales measured in nonequilibrium experiments is still an open question. Here, in this study, we investigate MgB 2 , in which conventional superconductivity at temperatures as high as 39 K is mediated by the strong coupling between the conduction electrons and the 𝐸 2⁢𝑔 phonon mode. By means of broadband time-resolved optical spectroscopy, we show that this selective electron-phonon coupling dictates the nonequilibrium optical response of MgB 2 at early times (<100 fs) after photoexcitation. Furthermore, based on an effective temperature model analysis, we estimate its contribution to the total electron-boson coupling function extracted from complementary equilibrium spectroscopy approaches, namely, optical reflectivity and angle-resolved photoemission spectroscopy. The coupling strength with the 𝐸 2⁢𝑔 phonon modes is thus estimated to be 𝜆 ≃ 0.56, which is approximately half of the total coupling constant, in agreement with 𝑎𝑏 initio calculations from the literature. As a benchmark, broadband time-resolved optical spectroscopy is performed also on the isostructural and nonsuperconducting compound AlB 2 , showing that the nonequilibrium optical response relaxes on a slower timescale due to the lack of strongly coupled phonon modes. Our findings demonstrate the possibility to resolve and quantify selective electron-phonon coupling from nonequilibrium optical spectroscopy.

Mor, Selene [Università Cattolica del Sacro Cuore,↗

Benchmarking the GW Approximation and Bethe–Salpeter Equation for Groups IB and IIB Atoms and Monoxides

Energies from the GW approximation and the Bethe–Salpeter equation (BSE) are benchmarked against the excitation energies of transition-metal (Cu, Zn, Ag, and Cd) single atoms and monoxide anions. We demonstrate that best estimates of GW quasiparticle energies at the complete basis set limit should be obtained via extrapolation or closure relations, while numerically converged GW -BSE eigenvalues can be obtained on a finite basis set. Calculations using real-space wave functions and pseudopotentials are shown to give best-estimate GW energies that agree (up to the extrapolation error) with calculations using all-electron Gaussian basis sets. We benchmark the effects of a vertex approximation (ΓLDA) and the mean-field starting point in GW and the BSE, performing computations using a real-space, transition-space basis and scalar-relativistic pseudopotentials. Here, while no variant of GW improves on perturbative G 0 W 0 at predicting ionization energies, G 0 W 0 Γ LDA -BSE computations give excellent agreement with experimental absorption spectra as long as off-diagonal self-energy terms are included. We also present G 0 W 0 quasiparticle energies for the CuO – , ZnO – , AgO – , and CdO – anions, in comparison to available anion photoelectron spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Reservoir Temperature Estimator (RTEst): A multicomponent geothermometry tool

The Reservoir Temperature Estimator (RTEst) is a multicomponent geothermometry tool for estimating reservoir geochemical parameters including reservoir temperature, CO 2 fugacity, mass of water lost or gained, and a reaction factor. It estimates these parameters and their associated uncertainties by minimizing an objective function that is the weighted sum of squares of the saturation indices of a user-selected set of minerals believed to be equilibrated with the reservoir fluid. RTEst accomplishes these estimates by combining the geochemical modeling capabilities of The Geochemist’s Workbench® with the optimization/parameter estimation resources of PEST®. An included interface aids the user in selecting plausible mineral phases to comprise the objective function and calculates their weighting factors. The working principles of RTEst are described and its efficacy is illustrated by presenting results of its application to various geothermal fields with known conditions. These examples show RTEst can account for the alteration of ascending reservoir fluid by mineral (calcite) re-equilibration with changes in temperature, reconstruct waters with CO 2 loss, correct for the deficit of water and other volatiles (CO 2,gas , H 2 S gas ) from boiling, and determine the amount of mixing of thermal and non-thermal waters. RTEst can use data with basis species below detection limit, missing, or unreliable either by assuming equilibrium with a controlling mineral (fixed-analyte method) or by treating the analyte concentration as an optimization parameter. The inverse of variance weighting method included in RTEst provides more representative results than either the normalization or unit weighting methods. Finally, the ability of RTEst to calculate reservoir temperatures, gas fugacity, and mixing fractions demonstrates its usefulness as a tool for evaluating geothermal systems.

15 GEOTHERMAL ENERGY↗

Constraining accuracy of the pairwise velocities in N -body simulations using scale-free models

We present a continuation of an analysis that aims to quantify resolution of N-body simulations by exploiting large (up to N = 4096 3 ) simulations of scale-free cosmologies run using abacus. Here, we focus on radial pairwise velocities of the matter field, both by direct estimation and through the cumulative two-point correlation function (using the pair conservation equation). We find that convergence at the 1 per cent level of the mean relative pairwise velocity can be demonstrated over a range of scales, evolving from a few times the grid spacing at early times to slightly below this scale at late times. We show the analysis of two different box sizes as well as from averaging results from the smaller boxes, and compare the power of the two aforementioned estimators in constraining accuracy at each scale. Down to scales of the order of the smoothing parameter, convergence is obtained at ∼5 per cent precision, and shows a behaviour indicating asymptotic stable clustering. In conclusion, we also infer for LCDM simulations conservative estimates on the evolution of the lower cut-off to resolution (at 1 and 5 per cent precision) as a function of redshift.

79 ASTRONOMY AND ASTROPHYSICS↗

ON THE ROTATION SPEED OF THE MILKY WAY DETERMINED FROM H i EMISSION

The circular rotation speed of the Milky Way at the solar radius, Θ{sub 0}, has been estimated to be 220 km s{sup −1} by fitting the maximum velocity of H i emission as a function of Galactic longitude. This result is in tension with a recent estimate of Θ{sub 0} = 240 km s{sup −1}, based on Very Long Baseline Interferometry (VLBI) parallaxes and proper motions from the BeSSeL and VERA surveys for large numbers of high-mass star-forming regions across the Milky Way. We find that the rotation curve best fitted to the VLBI data is slightly curved, and that this curvature results in a biased estimate of Θ{sub 0} from the H i data when a flat rotation curve is assumed. This relieves the tension between the methods and favors Θ{sub 0} = 240 km s{sup −1}.

79 ASTRONOMY AND ASTROPHYSICS↗

Efficient Nighttime Biogenic SOA Formation in a Polluted Residual Layer

Abstract Organic nitrates formed from nighttime reaction between anthropogenic nitrate radicals (NO3) and biogenic volatile organic compounds (BVOCs) are potentially an important but highly uncertain source of secondary organic aerosol (SOA). Here we report on the enhanced nighttime biogenic SOA formation observed in a polluted residual layer over Sacramento, California, the morning of 15 June 2010 during the Carbonaceous Aerosols and Radiative Effects Study (CARES). Trajectory analysis showed that the residual layer air, containing trace amounts of isoprene and monoterpenes left over after nighttime oxidation, was influenced by the San Francisco Bay Area emissions the previous evening. The residual layer aerosol was also enriched in nitrate, with about 64% of it estimated to be in the form of organic nitrates. The nitrate:organic mass ratio of the SOA was about 0.47±0.044, which corresponds to the range typically found in isoprene mononitrates. Assuming the SOA was composed of organic mononitrates, its nominal molecular weight was estimated at 186±11 g mol-1, consistent with the highly functionalized isoprene hydroxynitrates that have been observed in the particle phase in the southeast United States. Overall, our findings show that the efficiency of nighttime biogenic SOA formation, expressed as the change in organic aerosol mass relative to carbon monoxide (?OA/?CO), was about 100 µg m-3 ppmv, which falls in the range 77-157 µg m-3 ppmv previously estimated for enhanced daytime SOA formation from mixed anthropogenic and biogenic emissions during CARES.

Zaveri, Rahul A.↗

From molecular to macroscopic: predicting liquid–liquid phase equilibria and small-angle scattering of mixtures of organic liquids from atomistic simulation using Kirkwood–Buff theory

Macroscopic phase equilibria between solutions define the functionality of many biological and industrial processes, yet they are challenging to predict due to the inherent complexity of liquids containing large molecules. This work introduces an approach for the purely predictive calculation of such phase equilibria in temperature-composition space from molecular dynamics (MD) simulations at one temperature in the single-phase region. We use an approach developed previously to obtain the entropic and enthalpic contributions to the free energy of mixing from the atomic-scale information given by MD simulations via Kirkwood–Buff theory. This allows us to accurately estimate the free energy of mixing as a function of temperature, and thus obtain liquid–liquid phase equilibria, including liquid–liquid critical points, associated binodal and spinodal lines, and composition fluctuations across a region of temperature and composition. Results for binary malonamide–alkane systems are validated by comparison to a direct experimental probe of the fluctuations: the small angle X-ray scattering intensity near zero wavenumber. The MDKB → Phase method demonstrated here provides a significant improvement in predicting liquid–liquid equilibria and free energy as a function of temperature for our systems of interest compared to conventional thermodynamic models. The accurate performance of this purely predictive approach lies in its preservation of atomistic details when determining thermodynamic properties. Furthermore, its inherent extensibility to multi-component systems will likely make the MDKB → Phase approach a valuable general tool for connecting molecular interactions to macroscopic phase equilibria and for the computational screening of materials for targeted thermodynamic behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Long-term trends in urban NO 2 concentrations and associated paediatric asthma incidence: estimates from global datasets

Combustion-related nitrogen dioxide (NO 2 ) air pollution is associated with paediatric asthma incidence. We aimed to estimate global surface NO 2 concentrations consistent with the Global Burden of Disease study for 1990-2019 at a 1 km resolution, and the concentrations and attributable paediatric asthma incidence trends in 13 189 cities from 2000 to 2019. We scaled an existing annual average NO 2 concentration dataset for 2010-12 from a land use regression model (based on 5220 NO 2 monitors in 58 countries and land use variables) to other years using NO 2 column densities from satellite and reanalysis datasets. We applied these concentrations in an epidemiologically derived concentration-response function with population and baseline asthma rates to estimate NO 2 -attributable paediatric asthma incidence. We estimated that 1.85 million (95% uncertainty interval [UI] 0.93-2.80 million) new paediatric asthma cases were attributable to NO 2 globally in 2019, two thirds of which occurred in urban areas (1.22 million cases; 95% UI 0.60-1.8 million). The proportion of paediatric asthma incidence that is attributable to NO 2 in urban areas declined from 19.8% (1.22 million attributable cases of 6.14 million total cases) in 2000 to 16.0% (1.24 million attributable cases of 7.73 million total cases) in 2019. Urban attributable fractions dropped in high-income countries (-41%), Latin America and the Caribbean (-16%), central Europe, eastern Europe, and central Asia (-13%), and southeast Asia, east Asia, and Oceania (-6%), and rose in south Asia (+23%), sub-Saharan Africa (+11%), and north Africa and the Middle East (+5%). The contribution of NO 2 concentrations, paediatric population size, and asthma incidence rates to the change in NO 2 -attributable paediatric asthma incidence differed regionally. Despite improvements in some regions, combustion-related NO 2 pollution continues to be an important contributor to paediatric asthma incidence globally, particularly in cities. Mitigating air pollution should be a crucial element of public health strategies for children.

60 APPLIED LIFE SCIENCES↗

Extending the SAGA Survey (xSAGA). I. Satellite Radial Profiles as a Function of Host-galaxy Properties

We present "Extending the Satellites Around Galactic Analogs Survey" (xSAGA), a method for identifying low-z galaxies on the basis of optical imaging and results on the spatial distributions of xSAGA satellites around host galaxies. Using spectroscopic redshift catalogs from the SAGA Survey as a training data set, we have optimized a convolutional neural network (CNN) to identify z < 0.03 galaxies from more-distant objects using image cutouts from the DESI Legacy Imaging Surveys. From the sample of >100,000 CNN-selected low-z galaxies, we identify >20,000 probable satellites located between 36–300 projected kpc from NASA-Sloan Atlas central galaxies in the stellar-mass range $9.5\lt \mathrm{log}({M}_{\star }/{M}_{\odot })\lt 11$. We characterize the incompleteness and contamination for CNN-selected samples and apply corrections in order to estimate the true number of satellites as a function of projected radial distance from their hosts. Satellite richness depends strongly on host stellar mass, such that more-massive host galaxies have more satellites, and on host morphology, such that elliptical hosts have more satellites than disky hosts with comparable stellar masses. We also find a strong inverse correlation between satellite richness and the magnitude gap between a host and its brightest satellite. The normalized satellite radial distribution between 36–300 kpc does not depend on host stellar mass, morphology, or magnitude gap. The satellite abundances and radial distributions we measure are in reasonable agreement with predictions from hydrodynamic simulations. Our results deliver unprecedented statistical power for studying satellite galaxy populations and highlight the promise of using machine-learning for extending galaxy samples of wide-area surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Phase Imaging: A Compressive Sensing Approach

Since Wolfgang Pauli posed the question in 1933, whether the probability densities |Ψ(r)|² (real-space image) and |Ψ(q)|² (reciprocal space image) uniquely determine the wave function Ψ(r) [1], the so called Pauli Problem sparked numerous methods in all fields of microscopy [2, 3]. Reconstructing the complete wave function Ψ(r) = a(r)e-iφ(r) with the amplitude a(r) and the phase φ(r) from the recorded intensity enables the possibility to directly study the electric and magnetic properties of the sample through the phase. In transmission electron microscopy (TEM), electron holography is by far the most established method for phase reconstruction [4]. Requiring a high stability of the microscope, next to the installation of a biprism in the TEM, holography cannot be applied to any microscope straightforwardly. Recently, a phase retrieval approach was proposed using conventional TEM electron diffractive imaging (EDI). Using the SAD aperture as reciprocal-space constraint, a localized sample structure can be reconstructed from its diffraction pattern and a real-space image using the hybrid input-output algorithm [5]. We present an alternative approach using compressive phase-retrieval [6]. Our approach does not require a real-space image. Instead, random complimentary pairs of checkerboard masks are cut into a 200 nm Pt foil covering a conventional TEM aperture (cf. Figure 1). Used as SAD aperture, subsequently diffraction patterns are recorded from the same sample area. Hereby every mask blocks different parts of gold particles on a carbon support (cf. Figure 2). The compressive sensing problem has the following formulation. First, we note that the complex-valued reciprocal-space wave-function is the Fourier transform of the (also complex-valued) real-space wave-function, Ψ(q) = F[Ψ(r)], and subsequently the diffraction pattern image is given by |Ψ(q)|2 = |F[Ψ(r)]|2. We want to find Ψ(r) given a few differently coded diffraction pattern measurements yn = |F[HnΨ(r)]|2, where the matrices Hn encode the mask structure of the aperture. This is a nonlinear inverse problem, but has been shown to be solvable even in the underdetermined case [6]. Since each diffraction pattern yn contains diffraction information from selected regions of the same sample, the differences in each pattern contain local phase information, which can be combined to form a full estimate of the real-space wave-function[7]. References: [1] W. Pauli in “Die allgemeinen Prinzipien der Wellenmechanik“, ed. H Geiger and W Scheel, (Julius Springer, Berlin). [2] A. Tonomura, Rev. Mod. Phys. 59 (1987), p. 639. [3] J. Miao et al, Nature 400 (1999), p. 342. [4] H. Lichte et al, Annu. Rev. Mater. Res. 37 (2007), p. 539. [5] J. Yamasaki et al, Appl. Phys. Lett. 101 (2012), 234105. [6] P Schniter and S Rangan. Signal Proc., IEEE Trans. on. 64(4), (2015), pp. 1043. [7] Supported by the Chemical Imaging, Signature Discovery, and Analytics in Motion initiatives at PNNL. PNNL is operated by Battelle Memorial Inst. for the US DOE; contract DE-AC05-76RL01830.

Schneider, Sebastian↗

Transverse Λ polarization in unpolarized $pp$ → jet Λ ↑ X

In this Letter, we investigate the spontaneous transverse polarization of Λ hyperons produced in unpolarized pp collisions inside a jet, by adopting a TMD approach where transverse momentum effects are included only in the fragmentation process. We will present predictions based on the parametrizations of the Λ polarizing fragmentation function as extracted from fits to Belle e + e – data. These estimates will be compared against preliminary STAR data. We will then be able to explore the universality properties of the quark polarizing fragmentation function and, for the first time, the role of its gluon counterpart.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hyper Open-Shell Excited Spin States of Transition-Metal Compounds: FeF 2 , FeF 2 ···Ethane, and FeF 2 ···Ethylene

Spin-state energetics are important for understanding properties that involve more than one spin state, for example, catalysis occurring on two or more potential energy surfaces corresponding to different electronic spins. Very often, multiple-spin processes involve transition-metal compounds, and therefore, it is important to understand the electronic structure and energetics of such compounds in different spin states. In this work, we benchmark relative spin-state energies of FeF 2 with respect to the quintet ground spin state using both singleconfigurational and multi-configurational methods, and we examine how they are affected by the binding of ethane and ethylene to the iron center. We also benchmark the binding energies of the complexes. The single-configurational methods used in this work are the Hartree–Fock method, 32 exchange–correlation functionals, and the CCSD(T) coupled-cluster method in both restricted and unrestricted formalisms. The multi-configurational methods that have been used are CASSCF, CASPT2, CASPT3, MRCI, MRCI+Q, and MR-ACPF. The spin-state splitting energies depend on the functional chosen and of the 32 exchange-correlation functionals investigated here we find that for the septet and spin-projected triplet states of FeF 2 , the M06 functional is the best when compared to our best estimates from multi-reference calculations. If all the nine excitation energies are considered, where there are three excited spin states (singlet, triplet, and septet) for each of the three systems (FeF 2 , FeF 2 ···ethane, and FeF 2 ···ethylene), the three best performing functionals are HLE16, SOGGA11-X, and M06-2X. We find that the binding of ethane perturbs the relative spin-state energy of FeF 2 by only a small amount, but the stronger binding of ethylene has a larger effect. For the spin-state splitting energies of FeF 2 using single-reference CCSD(T), we found the predicted results depend very strongly on precisely how the calculations are done, in particular, on the spin-restricted or spin-unrestricted character of the SCF reference state, which can differ even by around 50 kcal/mol for the SCF reference state and the subsequent CCSD(T) calculations. Furthermore, on analyzing the wave functions of both the spinrestricted or spin-unrestricted formalisms, we find that the lowest-energy singlet and triplet states of the complexes, just like FeF 2 in isolation, can have more unpaired electrons than they are usually assumed to have, i.e. they can be hyper open-shell electronic configurations, and that this can significantly lower the energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A three-dimensional model of terrain-induced updrafts for movement ecology studies

Abstract Background Spatially explicit simulation models of animal movements through the atmosphere necessarily require a representation of the spatial and temporal variation of atmospheric conditions. In particular, for movements of soaring birds that rely extensively on vertical updrafts to avoid flapping flight, accurate and reliable estimation of the vertical component of wind is critical. The interaction between wind and complex terrain shapes both the horizontal and vertical wind fields, highlighting the need to model the coupling between local terrain features and atmospheric conditions at scales relevant to animal movement. Methods In this work, we propose a new empirical model for estimating the orographic updraft field. The model is developed using computational fluid dynamics simulations of canonical atmospheric conditions over moderately complex terrain. To isolate buoyancy and thermal effects, and focus on terrain-induced effects, we use only simulations of a neutrally stratified atmosphere to develop the model. The model, which we name Engineering Vertical Velocity Estimator (EVVE), is simple to implement and is a function of the underlying terrain elevation map, the desired height above ground level (AGL), and wind conditions at a reference height (80 m). We validate the model with data from the Alaiz mountain (Spain) field campaign. Results Compared to observations, the proposed improved model estimates the updrafts at 120 m AGL with a mean error of 0.11 m/s ( $$\sigma =0.28$$ σ = 0.28 m/s), compared to 0.85 m/s ( $$\sigma =0.58$$ σ = 0.58 m/s) for its baseline. For typical land-based wind turbine hub heights of 80 m AGL, the proposed model has a mean error of 0.04 m/s ( $$\sigma =0.25$$ σ = 0.25 m/s), compared to baseline 0.54 m/s ( $$\sigma =0.45$$ σ = 0.45 m/s) estimations. We illustrate an application of the model in movement ecology by comparing simulated tracks and presence maps of golden eagles ( Aquila chrysaetos ) moving across two distinct landscapes. The tracks and presence maps are obtained using a simple heuristic-based movement model, with the updraft field given by the proposed model and a wind vector-based estimation approach that is currently in wide use in movement ecology studies of raptors and other soaring birds. Conclusions We highlight that movement model results can be sensitive to the underlying orographic updraft model, especially in studies of fine-scale movements in regions of complex topography. We suggest adopting the proposed model rather than the wind vector estimation method for studies of soaring bird movements.

17 WIND ENERGY↗