Crystal Coulomb energies: I. Organic donor-acceptor complexes--a review. II. Classical Ewald calculations of the Coulomb binding energy of four donor-acceptor crystals and Wurster's blue perchlorate. III. A rapid-convergence quantum-mechanical formalism for crystal electronic energies

Chapter I

Theories for organic donor-acceptor (DA) complexes in solution and in the solid state are reviewed, and compared with the available experimental data. As shown by McConnell et al. (Proc. Natl. Acad. Sci. U.S., 53, 46-50 (1965)), the DA crystals fall into two classes, the holoionic class with a fully or almost fully ionic ground state, and the nonionic class with little or no ionic character. If the total lattice binding energy 2ε₁ (per DA pair) gained in ionizing a DA lattice exceeds the cost 2ε_o of ionizing each DA pair, ε₁ + ε_o less than 0, then the lattice is holoionic. The charge-transfer (CT) band in crystals and in solution can be explained, following Mulliken, by a second-order mixing of states, or by any theory that makes the CT transition strongly allowed, and yet due to a small change in the ground state of the non-interacting components D and A (or D⁺ and A^-). The magnetic properties of the DA crystals are discussed.

Chapter II

A computer program, EWALD, was written to calculate by the Ewald fast-convergence method the crystal Coulomb binding energy E_C due to classical monopole-monopole interactions for crystals of any symmetry. The precision of E_C values obtained is high: the uncertainties, estimated by the effect on E_C of changing the Ewald convergence parameter η, ranged from ± 0.00002 eV to ± 0.01 eV in the worst case. The charge distribution for organic ions was idealized as fractional point charges localized at the crystallographic atomic positions: these charges were chosen from available theoretical and experimental estimates. The uncertainty in E_C due to different charge distribution models is typically ± 0.1 eV (± 3%): thus, even the simple Hückel model can give decent results.

E_C for Wurster's Blue Perchl orate is -4.1 eV/molecule: the crystal is stable under the binding provided by direct Coulomb interactions. E_C for N-Methylphenazinium Tetracyanoquino- dimethanide is 0.1 eV: exchange Coulomb interactions, which cannot be estimated classically, must provide the necessary binding.

EWALD was also used to test the McConnell classification of DA crystals. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine: 7,7,8,8-Tetracyanoquinodimethan) E_C = -4.0 eV while 2ε_o = 4.6₅ eV: clearly, exchange forces must provide the balance. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine:para-Chloranil) E_C = -4.4 eV, while 2ε_o = 5.0 eV: again EC falls short of 2ε₁. As a Gedankenexperiment, two nonionic crystals were assumed to be ionized: for (1:1)-(Hexamethyl- benzene:para-Chloranil) E_C = -4.5 eV, 2ε_o = 6.6 eV; for (1:1)- (Napthalene:Tetracyanoethylene) E_C = -4.3 eV, 2ε_o = 6.5 eV. Thus, exchange energies in these nonionic crystals must not exceed 1 eV.

Chapter III

A rapid-convergence quantum-mechanical formalism is derived to calculate the electronic energy of an arbitrary molecular (or molecular-ion) crystal: this provides estimates of crystal binding energies which include the exchange Coulomb inter- actions. Previously obtained LCAO-MO wavefunctions for the isolated molecule(s) ("unit cell spin-orbitals") provide the starting-point. Bloch's theorem is used to construct "crystal spin-orbitals". Overlap between the unit cell orbitals localized in different unit cells is neglected, or is eliminated by Löwdin orthogonalization. Then simple formulas for the total kinetic energy Q^(XT)_λ, nuclear attraction [λ/λ]^XT, direct Coulomb [λλ/λ'λ']^XT and exchange Coulomb [λλ'/λ'λ]^XT integrals are obtained, and direct-space brute-force expansions in atomic wavefunctions are given. Fourier series are obtained for [λ/λ]^XT, [λλ/λ'λ']^XT, and [λλ/λ'λ]^XT with the help of the convolution theorem; the Fourier coefficients require the evaluation of Silverstone's two-center Fourier transform integrals. If the short-range interactions are calculated by brute-force integrations in direct space, and the long-range effects are summed in Fourier space, then rapid convergence is possible for [λ/λ]^XT, [λλ/λ'λ']^XT and [λλ'/λ'λ]^XT. This is achieved, as in the Ewald method, by modifying each atomic wavefunction by a "Gaussian convergence acceleration factor", and evaluating separately in direct and in Fourier space appropriate portions of [λ/λ]^XT, etc., where some of the portions contain the Gaussian factor.

Divergences in the Long Run Trends in the Price of Energy and of Energy Services

32 p.

Sustainable IT and IT for sustainability

Energy and sustainability have become one of the most critical issues of our generation. While the abundant potential of renewable energy such as solar and wind provides a real opportunity for sustainability, their intermittency and uncertainty present a daunting operating challenge. This thesis aims to develop analytical models, deployable algorithms, and real systems to enable efficient integration of renewable energy into complex distributed systems with limited information.

The first thrust of the thesis is to make IT systems more sustainable by facilitating the integration of renewable energy into these systems. IT represents the fastest growing sectors in energy usage and greenhouse gas pollution. Over the last decade there are dramatic improvements in the energy efficiency of IT systems, but the efficiency improvements do not necessarily lead to reduction in energy consumption because more servers are demanded. Further, little effort has been put in making IT more sustainable, and most of the improvements are from improved "engineering" rather than improved "algorithms". In contrast, my work focuses on developing algorithms with rigorous theoretical analysis that improve the sustainability of IT. In particular, this thesis seeks to exploit the flexibilities of cloud workloads both (i) in time by scheduling delay-tolerant workloads and (ii) in space by routing requests to geographically diverse data centers. These opportunities allow data centers to adaptively respond to renewable availability, varying cooling efficiency, and fluctuating energy prices, while still meeting performance requirements. The design of the enabling algorithms is however very challenging because of limited information, non-smooth objective functions and the need for distributed control. Novel distributed algorithms are developed with theoretically provable guarantees to enable the "follow the renewables" routing. Moving from theory to practice, I helped HP design and implement industry's first Net-zero Energy Data Center.

The second thrust of this thesis is to use IT systems to improve the sustainability and efficiency of our energy infrastructure through data center demand response. The main challenges as we integrate more renewable sources to the existing power grid come from the fluctuation and unpredictability of renewable generation. Although energy storage and reserves can potentially solve the issues, they are very costly. One promising alternative is to make the cloud data centers demand responsive. The potential of such an approach is huge.

To realize this potential, we need adaptive and distributed control of cloud data centers and new electricity market designs for distributed electricity resources. My work is progressing in both directions. In particular, I have designed online algorithms with theoretically guaranteed performance for data center operators to deal with uncertainties under popular demand response programs. Based on local control rules of customers, I have further designed new pricing schemes for demand response to align the interests of customers, utility companies, and the society to improve social welfare.

Effects of starvation on energy density of juvenile chum salmon (Oncorhynchus keta) captured in marine waters of Southeastern Alaska

We conducted laboratory starvation experiments on juvenile chum salmon (Oncorhynchus keta) captured in the neritic marine waters of northern Southeast Alaska in June and July 2003. Temporal changes in fish energy density (whole body energy content [WBEC], cal/g dry weight), percent moisture content, wet weight (g), length (mm), and size-related condition residuals were measured in the laboratory and were then compared to long-term field data. Laboratory water temperatures and salinities averaged 9°C and 32 psu in both months. Trends in response variables were similar for both experimental groups, although sampling intervals were limited in July because fewer fish were available (n= 54) than in June (n=101). Overall, for June (45-d experimental period, 9 intervals), WBEC, wet weight, and condition residuals decreased and percent moisture content increased, whereas fork length did not change. For July (20-d experimental period, 5 intervals), WBEC and condition residuals decreased, percent moisture content and fork length increased, and wet weight did not change. WBEC, percent moisture content, and condition residuals fell outside the norm of longterm data ranges within 10–15 days of starvation, and may be more useful than fork length and wet weight for detecting fish condition responses to suboptimal environments.

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

468 p.

Study of the turbulence in the air-side and water-side boundary layers in experimental laboratory wind induced surface waves

This study detailed the structure of turbulence in the air-side and water-side boundary layers in wind-induced surface waves. Inside the air boundary layer, the kurtosis is always greater than 3 (the value for normal distribution) for both horizontal and vertical velocity fluctuations. The skewness for the horizontal velocity is negative, but the skewness for the vertical velocity is always positive. On the water side, the kurtosis is always greater than 3, and the skewness is slightly negative for the horizontal velocity and slightly positive for the vertical velocity. The statistics of the angle between the instantaneous vertical fluctuation and the instantaneous horizontal velocity in the air is similar to those obtained over solid walls. Measurements in water show a large variance, and the peak is biased towards negative angles. In the quadrant analysis, the contribution of quadrants Q2 and Q4 is dominant on both the air side and the water side. The non-dimensional relative contributions and the concentration match fairly well near the interface. Sweeps in the air side (belonging to quadrant Q4) act directly on the interface and exert pressure fluctuations, which, in addition to the tangential stress and form drag, lead to the growth of the waves. The water drops detached from the crest and accelerated by the wind can play a major role in transferring momentum and in enhancing the turbulence level in the water side.On the air side, the Reynolds stress tensor's principal axes are not collinear with the strain rate tensor, and show an angle α σ≈=-20°to-25°. On the water side, the angle is α σ≈=-40°to-45°. The ratio between the maximum and the minimum principal stresses is σ a/σ b=3to4 on the air side, and σ a/σ b=1.5to3 on the water side. In this respect, the air-side flow behaves like a classical boundary layer on a solid wall, while the water-side flow resembles a wake. The frequency of bursting on the water side increases significantly along the flow, which can be attributed to micro-breaking effects - expected to be more frequent at larger fetches. © 2012 Elsevier B.V.

Tests of an adaptive QM/MM calculation on free energy profiles of chemical reactions in solution.

We present reaction free energy calculations using the adaptive buffered force mixing quantum mechanics/molecular mechanics (bf-QM/MM) method. The bf-QM/MM method combines nonadaptive electrostatic embedding QM/MM calculations with extended and reduced QM regions to calculate accurate forces on all atoms, which can be used in free energy calculation methods that require only the forces and not the energy. We calculate the free energy profiles of two reactions in aqueous solution: the nucleophilic substitution reaction of methyl chloride with a chloride anion and the deprotonation reaction of the tyrosine side chain. We validate the bf-QM/MM method against a full QM simulation, and show that it correctly reproduces both geometrical properties and free energy profiles of the QM model, while the electrostatic embedding QM/MM method using a static QM region comprising only the solute is unable to do so. The bf-QM/MM method is not explicitly dependent on the details of the QM and MM methods, so long as it is possible to compute QM forces in a small region and MM forces in the rest of the system, as in a conventional QM/MM calculation. It is simple, with only a few parameters needed to control the QM calculation sizes, and allows (but does not require) a varying and adapting QM region which is necessary for simulating solutions.

Coupling model uncertainty for coupled rainfall/runoff and surface water quality models in river problems

Coupled hydrology and water quality models are an important tool today, used in the understanding and management of surface water and watershed areas. Such problems are generally subject to substantial uncertainty in parameters, process understanding, and data. Component models, drawing on different data, concepts, and structures, are affected differently by each of these uncertain elements. This paper proposes a framework wherein the response of component models to their respective uncertain elements can be quantified and assessed, using a hydrological model and water quality model as two exemplars. The resulting assessments can be used to identify model coupling strategies that permit more appropriate use and calibration of individual models, and a better overall coupled model response. One key finding was that an approximate balance of water quality and hydrological model responses can be obtained using both the QUAL2E and Mike11 water quality models. The balance point, however, does not support a particularly narrow surface response (or stringent calibration criteria) with respect to the water quality calibration data, at least in the case examined here. Additionally, it is clear from the results presented that the structural source of uncertainty is at least as significant as parameter-based uncertainties in areal models. © 2012 John Wiley & Sons, Ltd.

Design trade-offs for energy regeneration and control in vehicle suspensions

This paper investigates the fundamental trade-offs involved in designing energy-regenerative suspensions, in particular, focusing on efficiency of power extraction and its effect on vehicle dynamics and control. It is shown that typical regenerative devices making use of linear-to-rotational elements can be modelled as a parallel arrangement of an inerter and a dissipative admittance. Taking account of typical adjustable parameters of the generator, it is shown, for a given suspension damping coefficient, that the power efficiency ratio scales with inertance. For a typical passenger vehicle, it is shown that there is a feasible compromise, namely that good efficiency is achievable with an inertance value that is not detrimental to vehicle performance. A prototype is designed and tested with a resistive termination and experimental results show good agreement between ideal and experimental admittances. The possibility to use dynamic (rather than purely resistive) loads to improve vehicle control without limiting the energy recovery is discussed. © 2013 Copyright Taylor and Francis Group, LLC.

A thermalization energy analysis of the threshold voltage shift in amorphous indium gallium zinc oxide thin film transistors under simultaneous negative gate bias and illumination

It has been previously observed that thin film transistors (TFTs) utilizing an amorphous indium gallium zinc oxide (a-IGZO) semiconducting channel suffer from a threshold voltage shift when subjected to a negative gate bias and light illumination simultaneously. In this work, a thermalization energy analysis has been applied to previously published data on negative bias under illumination stress (NBIS) in a-IGZO TFTs. A barrier to defect conversion of 0.65-0.75 eV is extracted, which is consistent with reported energies of oxygen vacancy migration. The attempt-to-escape frequency is extracted to be 10 6-107 s-1, which suggests a weak localization of carriers in band tail states over a 20-40 nm distance. Models for the NBIS mechanism based on charge trapping are reviewed and a defect pool model is proposed in which two distinct distributions of defect states exist in the a-IGZO band gap: these are associated with states that are formed as neutrally charged and 2+ charged oxygen vacancies at the time of film formation. In this model, threshold voltage shift is not due to a defect creation process, but to a change in the energy distribution of states in the band gap upon defect migration as this allows a state formed as a neutrally charged vacancy to be converted into one formed as a 2+ charged vacancy and vice versa. Carrier localization close to the defect migration site is necessary for the conversion process to take place, and such defect migration sites are associated with conduction and valence band tail states. Under negative gate bias stressing, the conduction band tail is depleted of carriers, but the bias is insufficient to accumulate holes in the valence band tail states, and so no threshold voltage shift results. It is only under illumination that the quasi Fermi level for holes is sufficiently lowered to allow occupation of valence band tail states. The resulting charge localization then allows a negative threshold voltage shift, but only under conditions of simultaneous negative gate bias and illumination, as observed experimentally as the NBIS effect. © 2014 AIP Publishing LLC.

Effects of animal-plant protein ratio in extruded and expanded diets on nitrogen and energy budgets of juvenile Chinese soft-shelled turtle (Pelodiscus sinensis Wiegmann)

In this study, we investigated the effects of animal-plant protein ratio in extruded and expanded diets on nutrient digestibility, nitrogen and energy budgets of juvenile soft-shelled turtle (Pelodiscus sinensis). Four extruded and expanded feeds (diets 1-4) were formulated with different animal-plant protein ratios (diet 1, 1.50:1; diet 2, 2.95:1; diet 3, 4.92:1; diet 4, 7.29:1). The apparent digestibility coefficients (ADCs) of dry matter and crude lipid for diet 1 were significantly lower than those for diets 2-4. There was no significant difference in crude protein digestibility among diets 1-4. The ADC of carbohydrate was significantly increased with the increase in animal-plant protein. Although nitrogen intake rate, faecal nitrogen loss rate and excretory nitrogen loss rate of turtles fed diet 1 were significantly higher than those fed diets 2-4, nitrogen retention rate, net protein utilization and biological value of protein in these turtles were significantly lower than those fed diets 2-4. In addition, energy intake rate, excretory energy loss rate and heat production rate of turtles fed diet 1 were also significantly higher than those fed diets 2-4. Faecal energy loss was significantly reduced with the increase in the animal-plant protein ratio. The ADC of energy and assimilation efficiency of energy significantly increased with a higher animal-plant protein ratio. The growth efficiency of energy in the group fed diet 1 was significantly lower than those in the groups fed diets 2-4. Together, our results suggest that the optimum animal-plant protein ratio in extruded and expanded diets is around 3:1.

Structures of energy spectrum and persistent currents in mesoscopic rings with radial potential barrier in magnetic field

The energy spectrum and the persistent currents are calculated for finite-width mesoscopic annular structures with radial potential barrier in the presence of a magnetic field. The introduction of the tunneling barrier leads to the creation of extra edge states around the barrier and the occurrence of oscillatory structures superimposed on the bulk Landau level plateaus in the energy spectrum. We found that the Fermi energy E-F increases with the number of electrons N emerging many kinks. The single eigenstate persistent current exhibits complicated structures with vortex-like texture, ''bifurcation'', and multiple ''furcation'' patterns as N is increased. The total currents versus N display wild fluctuations.

Criterion of aerodynamic performance of large-scale offshore horizontal axis wind turbines

With the background of offshore wind energy projects, this paper studies aerodynamic performance and geometric characteristics of large capacity wind turbine rotors ( 1 to 10 MW), and the main characteristic parameters such as the rated wind speed, blade tip speed, and rotor solidity. We show that the essential criterion of a high-performance wind turbine is a highest possible annual usable energy pattern factor and a smallest possible dimension, capturing the maximum wind energy and producing the maximum annual power. The influence of the above-mentioned three parameters on the pattern factor and rotor geometry of wind turbine operated in China's offshore meteorological environment is investigated. The variation patterns of aerodynamic and geometric parameters are obtained, analyzed, and compared with each other. The present method for aerodynamic analysis and its results can form a basis for evaluating aerodynamic performance of large-scale offshore wind turbine rotors.

Multi-electron processes in C-13(6+) ions with neon collisions in energy region 4.15-11.08 keV/u

The cross-section ratios of double-, triple-, quadruple-, and the total multi-electron processes to the single electron capture process sigma(DE)/sigma(SC), sigma(TE)/sigma(SC), sigma(QE)/sigma(SC) and sigma(ME)/sigma(SC)) as well as the relative ratios among reaction channels in double-electron active, triple-electron active and quadruple- electron active are measured in C-13(6+) -Ne collision in the energy region of 4.15-11.08 keV/u by employing position-sensitive and time-of-flight coincident techniques. It is determined that the cross-section ratios sigma(DE)/sigma(SC), sigma(TE)/sigma(SC), sigma(QE)/sigma(SC) and sigma(ME)/sigma(SC) are approximately the constants of 0.20 +/- 0.03, 0.16 +/- 0.04, 0.06 +/- 0.02 and 0.42 +/- 0.05. These values are obviously smaller than the predictions of the molecular Coulomb over-the-barrier model (MCBM) [J. Phys. B 23 (1990) 4293], the extended classical over-the-barrier model (ECBM) [J. Phys. B 19 (1986) 2925] and the semiempirical scaling laws (SL) [Phys. Rev. A 54 (1996) 4127]. However, the relative ratios among partial processes of DE, TE and QE are found to depend on collision energy, which suggests that the collision dynamics depends on the collision velocity. The limitation of velocity-independent character of ECBM, MCBM and SL is undoubtedly shown.