Energy transfer in molecular collisions

Two general, numerically exact, quantum mechanical methods have been developed for the calculation of energy transfer in molecular collisions. The methods do not treat electronic transitions because of the exchange symmetry of the electrons. All interactions between the atoms in the system are written as potential energies.

The first method is a matrix generalization of the invariant imbedding procedure, ^{17, 20} adapted for multi-channel collision processes. The second method is based on a direct integration of the matrix Schrödinger equation, with a re-orthogonalization transform applied during the integration.

Both methods have been applied to a collinear collision model for two diatoms, interacting via a repulsive exponential potential. Two major studies were performed. The first was to determine the energy dependence of the transition probabilities for an H₂ on the H₂ model system. Transitions are possible between translational energy and vibrational energy, and from vibrational modes of one H₂ to the other H₂. The second study was to determine the variation of vibrational energy transfer probability with differences in natural frequency of two diatoms similar to N₂.

Comparisons were made to previous approximate analytical solutions of this same problem. For translational to vibrational energy transfer, the previous approximations were not adequate. For vibrational to vibrational energy transfer of one vibrational quantum, the approximations were quite good.

Stationary Eddies and Zonal Variations of the Global Hydrological Cycle in a Changing Climate

This thesis advances our physical understanding of the sensitivity of the hydrological cycle to global warming. Specifically, it focuses on changes in the longitudinal (zonal) variation of precipitation minus evaporation (P - E), which is predominantly controlled by planetary-scale stationary eddies. By studying idealized general circulation model (GCM) experiments with zonally varying boundary conditions, this thesis examines the mechanisms controlling the strength of stationary-eddy circulations and their role in the hydrological cycle. The overarching goal of this research is to understand the cause of changes in regional P - E with global warming. An understanding of such changes can be useful for impact studies focusing on water availability, ecosystem management, and flood risk.

Based on a moisture-budget analysis of ERA-Interim data, we establish an approximation for zonally anomalous P - E in terms of surface moisture content and stationary-eddy vertical motion in the lower troposphere. Part of the success of this approximation comes from our finding that transient-eddy moisture fluxes partially cancel the effect of stationary-eddy moisture advection, allowing divergent circulations to dominate the moisture budget. The lower-tropospheric vertical motion is related to horizontal motion in stationary eddies by Sverdrup and Ekman balance. These moisture- and vorticity-budget balances also hold in idealized and comprehensive GCM simulations across a range of climates.

By examining climate changes in the idealized and comprehensive GCM simulations, we are able to show the utility of the vertical motion P - E approximation for splitting changes in zonally anomalous P - E into thermodynamic and dynamic components. Shifts in divergent stationary-eddy circulations dominate changes in zonally anomalous P - E. This limits the local utility of the "wet gets wetter, dry gets drier” idea, where existing P - E patterns are amplified with warming by the increase in atmospheric moisture content, with atmospheric circulations held fixed. The increase in atmospheric moisture content manifests instead in an increase in the amplitude of the zonally anomalous hydrological cycle as measured by the zonal variance of P - E. However, dynamic changes, particularly the slowdown of divergent stationary-eddy circulations, limit the strengthening of the zonally anomalous hydrological cycle. In certain idealized cases, dynamic changes are even strong enough to reverse the tendency towards "wet gets wetter, dry gets drier” with warming.

Motivated by the importance of stationary-eddy vertical velocities in the moisture budget analysis, we examine controls on the amplitude of stationary eddies across a wide range of climates in an idealized GCM with simple topographic and ocean-heating zonal asymmetries. An analysis of the thermodynamic equation in the vicinity of topographic forcing reveals the importance of on-slope surface winds, the midlatitude isentropic slope, and latent heating in setting the amplitude of stationary waves. The response of stationary eddies to climate change is determined primarily by the strength of zonal surface winds hitting the mountain. The sensitivity of stationary-eddies to this surface forcing increases with climate change as the slope of midlatitude isentropes decreases. However, latent heating also plays an important role in damping the stationary-eddy response, and this damping becomes stronger with warming as the atmospheric moisture content increases. We find that the response of tropical overturning circulations forced by ocean heat-flux convergence is described by changes in the vertical structure of moist static energy and deep convection. This is used to derive simple scalings for the Walker circulation strength that capture the monotonic decrease with warming found in our idealized simulations.

Through the work of this thesis, the advances made in understanding the amplitude of stationary-waves in a changing climate can be directly applied to better understand and predict changes in the zonally anomalous hydrological cycle.

I. Baryon-antibaryon phase transition at high temperature. II. Inclusive virtual photon-hadron reactions in the parton model

Part I

Present experimental data on nucleon-antinucleon scattering allow a study of the possibility of a phase transition in a nucleon-antinucleon gas at high temperature. Estimates can be made of the general behavior of the elastic phase shifts without resorting to theoretical derivation. A phase transition which separates nucleons from antinucleons is found at about 280 MeV in the approximation of the second virial coefficient to the free energy of the gas.

Part II

The parton model is used to derive scaling laws for the hadrons observed in deep inelastic electron-nucleon scattering which lie in the fragmentation region of the virtual photon. Scaling relations are obtained in the Bjorken and Regge regions. It is proposed that the distribution functions become independent of both q² and ν where the Bjorken and Regge regions overlap. The quark density functions are discussed in the limit x→1 for the nucleon octet and the pseudoscalar mesons. Under certain plausible assumptions it is found that only one or two quarks of the six types of quarks and antiquarks have an appreciable density function in the limit x→1. This has implications for the quark fragmentation functions near the large momentum boundary of their fragmentation region. These results are used to propose a method of measuring the proton and neutron quark density functions for all x by making measurements on inclusively produced hadrons in electroproduction only. Implications are also discussed for the hadrons produced in electron-positron annihilation.

The statistical bootstrap model

A review is presented of the statistical bootstrap model of Hagedorn and Frautschi. This model is an attempt to apply the methods of statistical mechanics in high-energy physics, while treating all hadron states (stable or unstable) on an equal footing. A statistical calculation of the resonance spectrum on this basis leads to an exponentially rising level density ρ(m) ~ cm^-3 e^βom at high masses.

In the present work, explicit formulae are given for the asymptotic dependence of the level density on quantum numbers, in various cases. Hamer and Frautschi's model for a realistic hadron spectrum is described.

A statistical model for hadron reactions is then put forward, analogous to the Bohr compound nucleus model in nuclear physics, which makes use of this level density. Some general features of resonance decay are predicted. The model is applied to the process of NN annihilation at rest with overall success, and explains the high final state pion multiplicity, together with the low individual branching ratios into two-body final states, which are characteristic of the process. For more general reactions, the model needs modification to take account of correlation effects. Nevertheless it is capable of explaining the phenomenon of limited transverse momenta, and the exponential decrease in the production frequency of heavy particles with their mass, as shown by Hagedorn. Frautschi's results on "Ericson fluctuations" in hadron physics are outlined briefly. The value of β_o required in all these applications is consistently around [120 MeV]^-1 corresponding to a "resonance volume" whose radius is very close to ƛ_π. The construction of a "multiperipheral cluster model" for high-energy collisions is advocated.

On the Dynamics of the Adenylate Energy System: Homeorhesis vs Homeostasis

Biochemical energy is the fundamental element that maintains both the adequate turnover of the biomolecular structures and the functional metabolic viability of unicellular organisms. The levels of ATP, ADP and AMP reflect roughly the energetic status of the cell, and a precise ratio relating them was proposed by Atkinson as the adenylate energy charge (AEC). Under growth-phase conditions, cells maintain the AEC within narrow physiological values, despite extremely large fluctuations in the adenine nucleotides concentration. Intensive experimental studies have shown that these AEC values are preserved in a wide variety of organisms, both eukaryotes and prokaryotes. Here, to understand some of the functional elements involved in the cellular energy status, we present a computational model conformed by some key essential parts of the adenylate energy system. Specifically, we have considered (I) the main synthesis process of ATP from ADP, (II) the main catalyzed phosphotransfer reaction for interconversion of ATP, ADP and AMP, (III) the enzymatic hydrolysis of ATP yielding ADP, and (IV) the enzymatic hydrolysis of ATP providing AMP. This leads to a dynamic metabolic model (with the form of a delayed differential system) in which the enzymatic rate equations and all the physiological kinetic parameters have been explicitly considered and experimentally tested in vitro. Our central hypothesis is that cells are characterized by changing energy dynamics (homeorhesis). The results show that the AEC presents stable transitions between steady states and periodic oscillations and, in agreement with experimental data these oscillations range within the narrow AEC window. Furthermore, the model shows sustained oscillations in the Gibbs free energy and in the total nucleotide pool. The present study provides a step forward towards the understanding of the fundamental principles and quantitative laws governing the adenylate energy system, which is a fundamental element for unveiling the dynamics of cellular life.

Quantum ricochets: surface capture, release and energy loss of fast ions hitting a polar surface at grazing incidence

A diffraction mechanism is proposed for the capture, multiple bouncing and final escape of a fast ion (keV) impinging on the surface of a polarizable material at grazing incidence. Capture and escape are effected by elastic quantum diffraction consisting of the exchange of a parallel surface wave vector G= 2p/ a between the ion parallel momentum and the surface periodic potential of period a. Diffraction- assisted capture becomes possible for glancing angles F smaller than a critical value given by Fc 2- 2./ a-| Vim|/ E, where E is the kinetic energy of the ion,. = h/ Mv its de Broglie wavelength and Vim its average electronic image potential at the distance from the surface where diffraction takes place. For F< Fc, the ion can fall into a selected capture state in the quasi- continuous spectrum of its image potential and execute one or several ricochets before being released by the time reversed diffraction process. The capture, ricochet and escape are accompanied by a large, periodic energy loss of several tens of eV in the forward motion caused by the coherent emission of a giant number of quanta h. of Fuchs- Kliewer surface phonons characteristic of the polar material. An analytical calculation of the energy loss spectrum, based on the proposed diffraction process and using a model ion-phonon coupling developed earlier (Lucas et al 2013 J. Phys.: Condens. Matter 25 355009), is presented, which fully explains the experimental spectrum of Villette et al (2000 Phys. Rev. Lett. 85 3137) for Ne+ ions ricocheting on a LiF(001) surface.

Modeling Innovative Power Take-Off Based on Double-Acting Hydraulic Cylinders Array for Wave Energy Conversion

One of the key systems of a Wave Energy Converter for extraction of wave energy is the Power Take-Off (PTO) device. This device transforms the mechanical energy of a moving body into electrical energy. This paper describes the model of an innovative PTO based on independently activated double-acting hydraulic cylinders array. The model has been developed using a simulation tool, based on a port-based approach to model hydraulics systems. The components and subsystems used in the model have been parameterized as real components and their values experimentally obtained from an existing prototype. In fact, the model takes into account most of the hydraulic losses of each component. The simulations show the flexibility to apply different restraining torques to the input movement depending on the geometrical configuration and the hydraulic cylinders on duty, easily modified by a control law. The combination of these two actions allows suitable flexibility to adapt the device to different sea states whilst optimizing the energy extraction. The model has been validated using a real test bench showing good correlations between simulation and experimental tests.

Valuation of Wind Energy Projects: A Real Options Approach

We address the valuation of an operating wind farm and the finite-lived option to invest in it under different reward/support schemes: a constant feed-in tariff, a premium on top of the electricity market price (either a fixed premium or a variable subsidy such as a renewable obligation certificate or ROC), and a transitory subsidy, among others. Futures contracts on electricity with ever longer maturities enable market-based valuations to be undertaken. The model considers up to three sources of uncertainty: the electricity price, the level of wind generation, and the certificate (ROC) price where appropriate. When analytical solutions are lacking, we resort to a trinomial lattice combined with Monte Carlo simulation; we also use a two-dimensional binomial lattice when uncertainty in the ROC price is considered. Our data set refers to the UK. The numerical results show the impact of several factors involved in the decision to invest: the subsidy per MWh generated, the initial lump-sum subsidy, the maturity of the investment option, and electricity price volatility. Different combinations of variables can help bring forward investments in wind generation. One-off policies, e.g., a transitory initial subsidy, seem to have a stronger effect than a fixed premium per MWh produced.

Energy demands of diverse spiking cells from the neocortex, hippocampus, and thalamus

It has long been known that neurons in the brain are not physiologically homogeneous. In response to current stimulus, they can fire several distinct patterns of action potentials that are associated with different physiological classes ranging from regular-spiking cells, fast-spiking cells, intrinsically bursting cells, and low-threshold cells. In this work we show that the high degree of variability in firing characteristics of action potentials among these cells is accompanied with a significant variability in the energy demands required to restore the concentration gradients after an action potential. The values of the metabolic energy were calculated for a wide range of cell temperatures and stimulus intensities following two different approaches. The first one is based on the amount of Na+ load crossing the membrane during a single action potential, while the second one focuses on the electrochemical energy functions deduced from the dynamics of the computational neuron models. The results show that the thalamocortical relay neuron is the most energy-efficient cell consuming between 7 and 18 nJ/cm(2) for each spike generated, while both the regular and fast spiking cells from somatosensory cortex and the intrinsically-bursting cell from a cat visual cortex are the least energy-efficient, and can consume up to 100 nJ/cm(2) per spike. The lowest values of these energy demands were achieved at higher temperatures and high external stimuli.

A note on the von Bertalanffy growth function concerning the allocation of surplus energy to reproduction

We propose an extended form of the von Bertalanffy growth function (VBGF), where the allocation of surplus energy to reproduction is considered. Any function can be used in our model to describe the ratio of energy allocation for reproduction to that for somatic growth. As an example, two models for energy allocation were derived: a step-function and a logistic function. The extended model can jointly describe growth in adult and juvenile stages. The change in growth rate between the two stages can be either gradual or steep; the latter gives a biphasic VBGF. The results of curve fitting indicated that a consideration of reproductive energy is meaningful for model extension. By controlling parameter values, our comprehensive model gives various growth curve shapes ranging from indeterminate to determinate growth. An increase in the number of parameters is unavoidable in practical applications of this new model. Additional information on reproduction will improve the reliability of model estimates.

Energy retrofits in social housing. Analysis of its thermal behaviour.

468 p.

Hartle's model within the general theory of perturbatiye matchings: the change in mass

Hartle's model provides the most widely used analytic framework to describe isolated compact bodies rotating slowly in equilibrium up to second order in perturbations in the context of General Relativity. Apart from some explicit assumptions, there are some implicit, like the "continuity" of the functions in the perturbed metric across the surface of the body. In this work we sketch the basics for the analysis of the second order problem using the modern theory of perturbed matchings. In particular, the result we present is that when the energy density of the fluid in the static configuration does not vanish at the boundary, one of the functions of the second order perturbation in the setting of the original work by Hartle is not continuous. This discrepancy affects the calculation of the change in mass of the rotating star with respect to the static configuration needed to keep the central energy density unchanged.

Low Effort L-i Nuclear Fusion Plasma Control Using Model Predictive Control Laws

One of the main problems of fusion energy is to achieve longer pulse duration by avoiding the premature reaction decay due to plasma instabilities. The control of the plasma inductance arises as an essential tool for the successful operation of tokamak fusion reactors in order to overcome stability issues as well as the new challenges specific to advanced scenarios operation. In this sense, given that advanced tokamaks will suffer from limited power available from noninductive current drive actuators, the transformer primary coil could assist in reducing the power requirements of the noninductive current drive sources needed for current profile control. Therefore, tokamak operation may benefit from advanced control laws beyond the traditionally used PID schemes by reducing instabilities while guaranteeing the tokamak integrity. In this paper, a novel model predictive control (MPC) scheme has been developed and successfully employed to optimize both current and internal inductance of the plasma, which influences the L-H transition timing, the density peaking, and pedestal pressure. Results show that the internal inductance and current profiles can be adequately controlled while maintaining the minimal control action required in tokamak operation.

Effects of El Niño events on energy demand and egg production of rockfish (Scorpaenidae: Sebastes): a bioenergetics approach

Fish bioenergetics models estimate relationships between energy budgets and environmental and physiological variables. This study presents a generic rockfish (Sebastes) bioenergetics model and estimates energy consumption by northern California blue rockf ish (S. mystinus) under average (baseline) and El Niño conditions. Compared to males, female S. mystinus required more energy because they were larger and had greater reproductive costs. When El Niño conditions (warmer temperatures; lower growth, condition, and fecundity) were experienced every 3−7 years, energy consumption decreased on an individual and a per-recruit basis in relation to baseline conditions, but the decrease was minor (<4% at the individual scale, <7% at the per-recruit scale) compared to decreases in female egg production (12−19% at the individual scale, 15−23% at the per-recruit scale). When mortality in per-recruit models was increased by adding fishing, energy consumption in El Niño models grew more similar to that seen in the baseline model. However, egg production decreased significantly — an effect exacerbated by the frequency of El Niño events. Sensitivity analyses showed that energy consumption estimates were most sensitive to respiration parameters, energy density, and female fecundity, and that estimated consumption increased as parameter uncertainty increased. This model provides a means of understanding rockfish trophic ecology in the context of community structure and environmental change by synthesizing metabolic, demographic, and environmental information. Future research should focus on acquiring such information so that models like the bioenergetics model can be used to estimate the effect of climate change, community shifts, and different harvesting strategies on rockfish energy demands.