Exchange energy shifts under dense plasma conditions

The variation of the energy interval between the intercombination line ( 1s2p(P-3(1))-> 1s(2)) and the resonance line ( 1s2p(P-1(1))-> 1s(2)) of He-like aluminium with plasma density and temperature is investigated. Since such energy interval is equivalent to the exchange energy of the state 1s2p(P-3(1)), we consider the dependence of this energy shift on the plasma environment. It was found that the shifts of exchange energy increase ( decrease) with the increase of electron density ( electron temperature), and the shifts of exchange energy become more sensitive to the electron density as the electron temperature decreases, i. e. in the strongly coupled plasma regime. An approximately linear relation is found between the shifts of exchange energy and the electron density. The results show that dense plasma effects are very important for the simulation of the spectral fine structure. The relative shifts between the intercombination ( 1s2p(P-3(1))-> 1s(2)) and the resonance line ( 1s2p(P-1(1))-> 1s(2)) are discussed for diagnostic applications.

Black holes in the early universe, in compact binaries, and as energy sources inside solar-type stars

This thesis consists of three separate studies of roles that black holes might play in our universe.

In the first part we formulate a statistical method for inferring the cosmological parameters of our universe from LIGO/VIRGO measurements of the gravitational waves produced by coalescing black-hole/neutron-star binaries. This method is based on the cosmological distance-redshift relation, with "luminosity distances" determined directly, and redshifts indirectly, from the gravitational waveforms. Using the current estimates of binary coalescence rates and projected "advanced" LIGO noise spectra, we conclude that by our method the Hubble constant should be measurable to within an error of a few percent. The errors for the mean density of the universe and the cosmological constant will depend strongly on the size of the universe, varying from about 10% for a "small" universe up to and beyond 100% for a "large" universe. We further study the effects of random gravitational lensing and find that it may strongly impair the determination of the cosmological constant.

In the second part of this thesis we disprove a conjecture that black holes cannot form in an early, inflationary era of our universe, because of a quantum-field-theory induced instability of the black-hole horizon. This instability was supposed to arise from the difference in temperatures of any black-hole horizon and the inflationary cosmological horizon; it was thought that this temperature difference would make every quantum state that is regular at the cosmological horizon be singular at the black-hole horizon. We disprove this conjecture by explicitly constructing a quantum vacuum state that is everywhere regular for a massless scalar field. We further show that this quantum state has all the nice thermal properties that one has come to expect of "good" vacuum states, both at the black-hole horizon and at the cosmological horizon.

In the third part of the thesis we study the evolution and implications of a hypothetical primordial black hole that might have found its way into the center of the Sun or any other solar-type star. As a foundation for our analysis, we generalize the mixing-length theory of convection to an optically thick, spherically symmetric accretion flow (and find in passing that the radial stretching of the inflowing fluid elements leads to a modification of the standard Schwarzschild criterion for convection). When the accretion is that of solar matter onto the primordial hole, the rotation of the Sun causes centrifugal hangup of the inflow near the hole, resulting in an "accretion torus" which produces an enhanced outflow of heat. We find, however, that the turbulent viscosity, which accompanies the convective transport of this heat, extracts angular momentum from the inflowing gas, thereby buffering the torus into a lower luminosity than one might have expected. As a result, the solar surface will not be influenced noticeably by the torus's luminosity until at most three days before the Sun is finally devoured by the black hole. As a simple consequence, accretion onto a black hole inside the Sun cannot be an answer to the solar neutrino puzzle.

Time-resolved investigation of low-density plasma channels produced by a kilohertz femtosecond laser in air

By employing pump-probe back longitudinal diffractometry, the electron density and decay dynamics of a weak plasma channel created by a 1-KHz fs laser in air has been investigated. With ultrashort laser pulses of 50 fs and low energy of 0.6 mJ, we observe weak plasma channels with a length similar to 2 cm in air. An analytical reconstruction method of electron density has been analyzed, which is sensitive to the phase shift and channel size. The electron density in the weak plasma channel is extracted to be about 4x10(16) cm(-3). The diameters of the plasma channel and the filament are about 50 and 150 mu m, respectively, and the measurable electron density can be extended to less than 10(15) cm(-3). Moreover, a different time-frequency technique called linearly chirped longitudinal diffractometry is proposed to time-resolved investigate ultrafast ionization dynamics of laser-irradiated gas, laser interaction with cluster beam, etc.

Influence of target thickness on the generation of high-density ion bunches by ultrashort circularly polarized laser pulses

Ion acceleration by ultrashort circularly polarized laser pulse in a solid-density target is investigated using two-dimensional particle-in-cell simulation. The ions are accelerated and compressed by the continuously extending space-charge field created by the evacuation and compression of the target electrons by the laser light pressure. For a sufficiently thin target, the accelerated and compressed ions can reach and exit from the rear surface as a high-density high-energy ion bunch. The peak ion energy depends on the target thickness and reaches maximum when the compressed ion layer can just reach the rear target surface. The compressed ion layer exhibits lateral striation which can be suppressed by using a sharp-rising laser pulse. (c) 2008 American Institute of Physics.

Direct acceleration of solid-density plasma bunch by ultraintense laser

The interaction of a petawatt laser with a small solid-density plasma bunch is studied by particle-in-cell simulation. It is shown that when irradiated by a laser of intensity >10(21) W/cm(2), a dense plasma bunch of micrometer size can be efficiently accelerated. The kinetic energy of the ions in the high-density region of the plasma bunch can exceed ten MeV at a density in the 10(23)-cm(-3) level. Having a flux density orders of magnitude higher than that of the traditional charged-particle pulses, the laser-accelerated plasma bunch can have a wide range of applications. In particular, such a dense energetic plasma bunch impinging on the compressed fuel in inertial fusion can significantly enhance the nuclear-reaction cross section and is thus a promising alternative for fast ignition.

Optimum hot electron production with low-density foams for laser fusion by fast ignition

We propose a foam cone-in-shell target design aiming at optimum hot electron production for the fast ignition. A thin low-density foam is proposed to cover the inner tip of a gold cone inserted in a fuel shell. An intense laser is then focused on the foam to generate hot electrons for the fast ignition. Element experiments demonstrate increased laser energy coupling efficiency into hot electrons without increasing the electron temperature and beam divergence with foam coated targets in comparison with solid targets. This may enhance the laser energy deposition in the compressed fuel plasma.

Nanostructured tungsten trioxide photoanodes for solar energy conversion

Nanostructured tungsten trioxide (WO₃) photoelectrodes are potential candidates for the anodic portion of an integrated solar water-splitting device that generates hydrogen fuel and oxygen from water. These nanostructured materials can potentially offer improved performance in photooxidation reactions compared to unstructured materials because of enhancements in light scattering, increases in surface area, and their decoupling of the directions of light absorption and carrier collection. To evaluate the presence of these effects and their contributions toward energy conversion efficiency, a variety of nanostructured WO₃ photoanodes were synthesized by electrodeposition within nanoporous templates and by anodization of tungsten foils. A robust fabrication process was developed for the creation of oriented WO₃ nanorod arrays, which allows for control nanorod diameter and length. Films of nanostructured WO₃ platelets were grown via anodization, the morphology of the films was controlled by the anodization conditions, and the current-voltage performance and spectral response properties of these films were studied. The observed photocurrents were consistent with the apparent morphologies of the nanostructured arrays. Measurements of electrochemically active surface area and other physical characteristics were correlated with observed differences in absorbance, external quantum yield, and photocurrent density for the anodized arrays. The capability to quantify these characteristics and relate them to photoanode performance metrics can allow for selection of appropriate structural parameters when designing photoanodes for solar energy conversion.

Short-pulse laser absorption in very steep plasma density gradients

An analytical fluid model for resonance absorption during the oblique incidence by femtosecond laser pulses on a small-scale-length density plasma [k(0)L is an element of(0.1,10)] is proposed. The physics of resonance absorption is analyzed more clearly as we separate the electric field into an electromagnetic part and an electrostatic part. It is found that the characteristics of the physical quantities (fractional absorption, optimum angle, etc.) in a small-scale-length plasma are quite different from the predictions of classical theory. Absorption processes are generally dependent on the density scale length. For shorter scale length or higher laser intensity, vacuum heating tends to be dominant. It is shown that the electrons being pulled out and then returned to the plasma at the interface layer by the wave field can lead to a phenomenon like wave breaking. This can lead to heating of the plasma at the expanse of the wave energy. It is found that the optimum angle is independent of the laser intensity while the absorption rate increases with the laser intensity, and the absorption rate can reach as high as 25%. (c) 2006 American Institute of Physics.

Ground-state energy shifts of H-like Ti under dense and hot plasma conditions

In this study, by adopting the ion sphere model, the self-consistent. field method is used with the Poisson-Boltzmann equation and the Dirac equation to calculate the ground-state energies of H-like Ti at a plasma electron density from 10(22) cm(-3) to 10(24) cm(-3) and the electron temperature from 100 eV to 3600 eV. The ground-state energy shifts of H-like Ti show different trends with the electron density and the electron temperature. It is shown that the energy shifts increase with the increase in the electron density and decrease with the increase in the electron temperature. The energy shifts are sensitive to the electron density, but only sensitive to the low electron temperature. In addition, an accurately fitting formula is obtained to fast estimate the ground-state energies of H-like Ti. Such fitted formula can also be used to estimate the critical electron density of pressure ionization for the ground state of H-like Ti.

Density functional theory embedding for correlated wavefunctions

Methods that exploit the intrinsic locality of molecular interactions show significant promise in making tractable the electronic structure calculation of large-scale systems. In particular, embedded density functional theory (e-DFT) offers a formally exact approach to electronic structure calculations in which the interactions between subsystems are evaluated in terms of their electronic density. In the following dissertation, methodological advances of embedded density functional theory are described, numerically tested, and applied to real chemical systems.

First, we describe an e-DFT protocol in which the non-additive kinetic energy component of the embedding potential is treated exactly. Then, we present a general implementation of the exact calculation of the non-additive kinetic potential (NAKP) and apply it to molecular systems. We demonstrate that the implementation using the exact NAKP is in excellent agreement with reference Kohn-Sham calculations, whereas the approximate functionals lead to qualitative failures in the calculated energies and equilibrium structures.

Next, we introduce density-embedding techniques to enable the accurate and stable calculation of correlated wavefunction (CW) in complex environments. Embedding potentials calculated using e-DFT introduce the effect of the environment on a subsystem for CW calculations (WFT-in-DFT). We demonstrate that WFT-in-DFT calculations are in good agreement with CW calculations performed on the full complex.

We significantly improve the numerics of the algorithm by enforcing orthogonality between subsystems by introduction of a projection operator. Utilizing the projection-based embedding scheme, we rigorously analyze the sources of error in quantum embedding calculations in which an active subsystem is treated using CWs, and the remainder using density functional theory. We show that the embedding potential felt by the electrons in the active subsystem makes only a small contribution to the error of the method, whereas the error in the nonadditive exchange-correlation energy dominates. We develop an algorithm which corrects this term and demonstrate the accuracy of this corrected embedding scheme.

Advanced density functional theory methods for materials science

In this work we chiefly deal with two broad classes of problems in computational materials science, determining the doping mechanism in a semiconductor and developing an extreme condition equation of state. While solving certain aspects of these questions is well-trodden ground, both require extending the reach of existing methods to fully answer them. Here we choose to build upon the framework of density functional theory (DFT) which provides an efficient means to investigate a system from a quantum mechanics description.

Zinc Phosphide (Zn₃P₂) could be the basis for cheap and highly efficient solar cells. Its use in this regard is limited by the difficulty in n-type doping the material. In an effort to understand the mechanism behind this, the energetics and electronic structure of intrinsic point defects in zinc phosphide are studied using generalized Kohn-Sham theory and utilizing the Heyd, Scuseria, and Ernzerhof (HSE) hybrid functional for exchange and correlation. Novel 'perturbation extrapolation' is utilized to extend the use of the computationally expensive HSE functional to this large-scale defect system. According to calculations, the formation energy of charged phosphorus interstitial defects are very low in n-type Zn₃P₂ and act as 'electron sinks', nullifying the desired doping and lowering the fermi-level back towards the p-type regime. Going forward, this insight provides clues to fabricating useful zinc phosphide based devices. In addition, the methodology developed for this work can be applied to further doping studies in other systems.

Accurate determination of high pressure and temperature equations of state is fundamental in a variety of fields. However, it is often very difficult to cover a wide range of temperatures and pressures in an laboratory setting. Here we develop methods to determine a multi-phase equation of state for Ta through computation. The typical means of investigating thermodynamic properties is via ’classical’ molecular dynamics where the atomic motion is calculated from Newtonian mechanics with the electronic effects abstracted away into an interatomic potential function. For our purposes, a ’first principles’ approach such as DFT is useful as a classical potential is typically valid for only a portion of the phase diagram (i.e. whatever part it has been fit to). Furthermore, for extremes of temperature and pressure quantum effects become critical to accurately capture an equation of state and are very hard to capture in even complex model potentials. This requires extending the inherently zero temperature DFT to predict the finite temperature response of the system. Statistical modelling and thermodynamic integration is used to extend our results over all phases, as well as phase-coexistence regions which are at the limits of typical DFT validity. We deliver the most comprehensive and accurate equation of state that has been done for Ta. This work also lends insights that can be applied to further equation of state work in many other materials.

A model for energy and morphology of crystalline grain boundaries with arbitrary geometric character

It has been well-established that interfaces in crystalline materials are key players in the mechanics of a variety of mesoscopic processes such as solidification, recrystallization, grain boundary migration, and severe plastic deformation. In particular, interfaces with complex morphologies have been observed to play a crucial role in many micromechanical phenomena such as grain boundary migration, stability, and twinning. Interfaces are a unique type of material defect in that they demonstrate a breadth of behavior and characteristics eluding simplified descriptions. Indeed, modeling the complex and diverse behavior of interfaces is still an active area of research, and to the author's knowledge there are as yet no predictive models for the energy and morphology of interfaces with arbitrary character. The aim of this thesis is to develop a novel model for interface energy and morphology that i) provides accurate results (especially regarding "energy cusp" locations) for interfaces with arbitrary character, ii) depends on a small set of material parameters, and iii) is fast enough to incorporate into large scale simulations.

In the first half of the work, a model for planar, immiscible grain boundary is formulated. By building on the assumption that anisotropic grain boundary energetics are dominated by geometry and crystallography, a construction on lattice density functions (referred to as "covariance") is introduced that provides a geometric measure of the order of an interface. Covariance forms the basis for a fully general model of the energy of a planar interface, and it is demonstrated by comparison with a wide selection of molecular dynamics energy data for FCC and BCC tilt and twist boundaries that the model accurately reproduces the energy landscape using only three material parameters. It is observed that the planar constraint on the model is, in some cases, over-restrictive; this motivates an extension of the model.

In the second half of the work, the theory of faceting in interfaces is developed and applied to the planar interface model for grain boundaries. Building on previous work in mathematics and materials science, an algorithm is formulated that returns the minimal possible energy attainable by relaxation and the corresponding relaxed morphology for a given planar energy model. It is shown that the relaxation significantly improves the energy results of the planar covariance model for FCC and BCC tilt and twist boundaries. The ability of the model to accurately predict faceting patterns is demonstrated by comparison to molecular dynamics energy data and experimental morphological observation for asymmetric tilt grain boundaries. It is also demonstrated that by varying the temperature in the planar covariance model, it is possible to reproduce a priori the experimentally observed effects of temperature on facet formation.

Finally, the range and scope of the covariance and relaxation models, having been demonstrated by means of extensive MD and experimental comparison, future applications and implementations of the model are explored.

Nuclear reactions leading to energy levels of Ar37 and K37

In view of recent interest in the Cl³⁷ (ʋ solar’^e-)Ar³⁷ reaction cross section, information on some aspects of mass 37 nuclei has been obtained using the K³⁹ (d, ∝) Ar³⁷ and Cl³⁵ (He³, p) Ar³⁷ reactions. Ar³⁷ levels have been found at 0, 1.41, 1.62, 2.22, 2.50, 2.80, 3.17, 3.27, 3.53, 3.61, 3.71, (3.75), (3.90), 3.94, 4.02, (4.21), 4.28, 4.32, 4.40, 4.45, 4.58, 4.63, 4.74, 4.89, 4.98, 5.05, 5.10, 5.13, 5.21, 5.35, 5.41, 5.44, 5.54, 5.58, 5.67, 5.77, and 5.85 MeV (the underlined values correspond to previously tabulated levels). The nuclear temperature calculated from the Ar³⁷ level density is 1.4 MeV. Angular distributions of the lowest six levels with the K³⁹ (d, ∝) Ar³⁷ reaction at E_d = 10 MeV indicate a dominant direct interaction mechanism and the inapplicability of the 2I + 1 rule of the statistical model. Comparison of the spectra obtained with the K³⁹ (d, ∝) Ar³⁷ and Cl³⁵ (He³, p) Ar³⁷ reactions leads to the suggestion that the 5.13-MeV level is the T = 3/2 Cl³⁷ ground state analog. The ground state Q-value of the Ca⁴⁰ (p, ∝) K³⁷ reaction has been measured: -5179 ± 9 keV. This value implies a K³⁷ mass excess of -24804 ± 10 keV. Description of a NMR magnetometer and a sixteen-detector array used in conjunction with a 61-cm double-focusing magnetic spectrometer are included in appendices.

Energy loss of light ions channeling in silicon

Experimental measurements of rate of energy loss were made for protons of energy .5 to 1.6 MeV channeling through 1 μm thick silicon targets along the <110>, <111>, and <211> axial directions, and the {100}, {110}, {111}, and {211} planar directions. A .05% resolution automatically controlled magnetic spectrometer was used. The data are presented graphically along with an extensive summary of data in the literature. The data taken cover a wider range of channels than has previously been examined, and are in agreement with the data of F. Eisen, et al., Radd. Eff. 13, 93 (1972).

The theory in the literature for channeling energy loss due to interaction with local electrons, core electrons, and distant valence electrons of the crystal atoms is summarized. Straggling is analyzed, and a computer program which calculates energy loss and straggling using this theory and the Moliere approximation to the Thomas Fermi potential, V_TF, and the detailed silicon crystal structure is described. Values for the local electron density Z_loc in each of the channels listed above are extracted from the data by graphical matching of the experimental and computer results.

Zeroth and second order contributions to Z_loc as a function of distance from the center of the channel were computed from ∇²V_TF = 4πρ for various channels in silicon. For data taken in this work and data of F. Eisen, et al., Rad. Eff. 13, 93 (1972), the calculated zeroth order contribution to Z_loc lies between the experimentally extracted Z_loc values obtained by using the peak and the leading edge of the transmission spectra, suggesting that the observed straggling is due both to statistical fluctuations and to path variation.

I. The structure and properties of closed circular duplex DNA. II. Chemically interacting systems at equilibrium in a buoyant density gradient

I. The binding of the intercalating dye ethidium bromide to closed circular SV 40 DNA causes an unwinding of the duplex structure and a simultaneous and quantitatively equivalent unwinding of the superhelices. The buoyant densities and sedimentation velocities of both intact (I) and singly nicked (II) SV 40 DNAs were measured as a function of free dye concentration. The buoyant density data were used to determine the binding isotherms over a dye concentration range extending from 0 to 600 µg/m1 in 5.8 M CsCl. At high dye concentrations all of the binding sites in II, but not in I, are saturated. At free dye concentrations less than 5.4 µg/ml, I has a greater affinity for dye than II. At a critical amount of dye bound I and II have equal affinities, and at higher dye concentration I has a lower affinity than II. The number of superhelical turns, τ, present in I is calculated at each dye concentration using Fuller and Waring's (1964) estimate of the angle of duplex unwinding per intercalation. The results reveal that SV 40 DNA I contains about -13 superhelical turns in concentrated salt solutions.

The free energy of superhelix formation is calculated as a function of τ from a consideration of the effect of the superhelical turns upon the binding isotherm of ethidium bromide to SV 40 DNA I. The value of the free energy is about 100 kcal/mole DNA in the native molecule. The free energy estimates are used to calculate the pitch and radius of the superhelix as a function of the number of superhelical turns. The pitch and radius of the native I superhelix are 430 Å and 135 Å, respectively.

A buoyant density method for the isolation and detection of closed circular DNA is described. The method is based upon the reduced binding of the intercalating dye, ethidium bromide, by closed circular DNA. In an application of this method it is found that HeLa cells contain in addition to closed circular mitochondrial DNA of mean length 4.81 microns, a heterogeneous group of smaller DNA molecules which vary in size from 0.2 to 3.5 microns and a paucidisperse group of multiples of the mitochondrial length.

II. The general theory is presented for the sedimentation equilibrium of a macromolecule in a concentrated binary solvent in the presence of an additional reacting small molecule. Equations are derived for the calculation of the buoyant density of the complex and for the determination of the binding isotherm of the reagent to the macrospecies. The standard buoyant density, a thermodynamic function, is defined and the density gradients which characterize the four component system are derived. The theory is applied to the specific cases of the binding of ethidium bromide to SV 40 DNA and of the binding of mercury and silver to DNA.