Strength of tantalum at high pressures through Richtmyer-Meshkov laser compression experiments and simulations

Strength at extreme pressures (>1 Mbar or 100 GPa) and high strain rates (106-108 s-1) of materials is not well characterized. The goal of the research outlined in this thesis is to study the strength of tantalum (Ta) at these conditions. The Omega Laser in the Laboratory for Laser Energetics in Rochester, New York is used to create such extreme conditions. Targets are designed with ripples or waves on the surface, and these samples are subjected to high pressures using Omega’s high energy laser beams. In these experiments, the observational parameter is the Richtmyer-Meshkov (RM) instability in the form of ripple growth on single-mode ripples. The experimental platform used for these experiments is the “ride-along” laser compression recovery experiments, which provide a way to recover the specimens having been subjected to high pressures. Six different experiments are performed on the Omega laser using single-mode tantalum targets at different laser energies. The energy indicates the amount of laser energy that impinges the target. For each target, values for growth factor are obtained by comparing the profile of ripples before and after the experiment. With increasing energy, the growth factor increased.

Engineering simulations are used to interpret and correlate the measurements of growth factor to a measure of strength. In order to validate the engineering constitutive model for tantalum, a series of simulations are performed using the code Eureka, based on the Optimal Transportation Meshfree (OTM) method. Two different configurations are studied in the simulations: RM instabilities in single and multimode ripples. Six different simulations are performed for the single ripple configuration of the RM instability experiment, with drives corresponding to laser energies used in the experiments. Each successive simulation is performed at higher drive energy, and it is observed that with increasing energy, the growth factor increases. Overall, there is favorable agreement between the data from the simulations and the experiments. The peak growth factors from the simulations and the experiments are within 10% agreement. For the multimode simulations, the goal is to assist in the design of the laser driven experiments using the Omega laser. A series of three-mode and four-mode patterns are simulated at various energies and the resulting growth of the RM instability is computed. Based on the results of the simulations, a configuration is selected for the multimode experiments. These simulations also serve as validation for the constitutive model and the material parameters for tantalum that are used in the simulations.

By designing samples with initial perturbations in the form of single-mode and multimode ripples and subjecting these samples to high pressures, the Richtmyer-Meshkov instability is investigated in both laser compression experiments and simulations. By correlating the growth of these ripples to measures of strength, a better understanding of the strength of tantalum at high pressures is achieved.

A quantitative investigation of the solar modulation of cosmic-ray protons and helium nuclei

The differential energy spectra of cosmic-ray protons and He nuclei have been measured at energies up to 315 MeV/nucleon using balloon- and satellite-borne instruments. These spectra are presented for solar quiet times for the years 1966 through 1970. The data analysis is verified by extensive accelerator calibrations of the detector systems and by calculations and measurements of the production of secondary protons in the atmosphere.

The spectra of protons and He nuclei in this energy range are dominated by the solar modulation of the local interstellar spectra. The transport equation governing this process includes as parameters the solar-wind velocity, V, and a diffusion coefficient, K(r,R), which is assumed to be a scalar function of heliocentric radius, r, and magnetic rigidity, R. The interstellar spectra, j_D, enter as boundary conditions on the solutions to the transport equation. Solutions to the transport equation have been calculated for a broad range of assumed values for K(r,R) and j_D and have been compared with the measured spectra.

It is found that the solutions may be characterized in terms of a dimensionless parameter, ψ(r,R) = ^∞∫_r V dr'/(K(r',R). The amount of modulation is roughly proportional to ψ. At high energies or far from the Sun, where the modulation is weak, the solution is determined primarily by the value of ψ (and the interstellar spectrum) and is not sensitive to the radial dependence of the diffusion coefficient. At low energies and for small r, where the effects of adiabatic deceleration are found to be large, the spectra are largely determined by the radial dependence of the diffusion coefficient and are not very sensitive to the magnitude of ψ or to the interstellar spectra. This lack of sensitivity to j_D implies that the shape of the spectra at Earth cannot be used to determine the interstellar intensities at low energies.

Values of ψ determined from electron data were used to calculate the spectra of protons and He nuclei near Earth. Interstellar spectra of the form j_D α (W - 0.25m)^-2.65 for both protons and He nuclei were found to yield the best fits to the measured spectra for these values of ψ, where W is the total energy and m is the rest energy. A simple model for the diffusion coefficient was used in which the radial and rigidity dependence are separable and K is independent of radius inside a modulation region which has a boundary at a distance D. Good agreement was found between the measured and calculated spectra for the years 1965 through 1968, using typical boundary distances of 2.7 and 6.1 A.U. The proton spectra observed in 1969 and 1970 were flatter than in previous years. This flattening could be explained in part by an increase in D, but also seemed to require that a noticeable fraction of the observed protons at energies as high at 50 to 100 MeV be attributed to quiet-time solar emission. The turnup in the spectra at low energies observed in all years was also attributed to solar emission. The diffusion coefficient used to fit the 1965 spectra is in reasonable agreement with that determined from the power spectra of the interplanetary magnetic field (Jokipii and Coleman, 1968). We find a factor of roughly 3 increase in ψ from 1965 to 1970, corresponding to the roughly order of magnitude decrease in the proton intensity at 250 MeV. The change in ψ might be attributed to a decrease in the diffusion coefficient, or, if the diffusion coefficient is essentially unchanged over that period (Mathews et al., 1971), might be attributed to an increase in the boundary distance, D.

First principles based multiparadigm modeling of electronic structures and dynamics

Electronic structures and dynamics are the key to linking the material composition and structure to functionality and performance.

An essential issue in developing semiconductor devices for photovoltaics is to design materials with optimal band gaps and relative positioning of band levels. Approximate DFT methods have been justified to predict band gaps from KS/GKS eigenvalues, but the accuracy is decisively dependent on the choice of XC functionals. We show here for CuInSe₂ and CuGaSe₂, the parent compounds of the promising CIGS solar cells, conventional LDA and GGA obtain gaps of 0.0-0.01 and 0.02-0.24 eV (versus experimental values of 1.04 and 1.67 eV), while the historically first global hybrid functional, B3PW91, is surprisingly the best, with band gaps of 1.07 and 1.58 eV. Furthermore, we show that for 27 related binary and ternary semiconductors, B3PW91 predicts gaps with a MAD of only 0.09 eV, which is substantially better than all modern hybrid functionals, including B3LYP (MAD of 0.19 eV) and screened hybrid functional HSE06 (MAD of 0.18 eV).

The laboratory performance of CIGS solar cells (> 20% efficiency) makes them promising candidate photovoltaic devices. However, there remains little understanding of how defects at the CIGS/CdS interface affect the band offsets and interfacial energies, and hence the performance of manufactured devices. To determine these relationships, we use the B3PW91 hybrid functional of DFT with the AEP method that we validate to provide very accurate descriptions of both band gaps and band offsets. This confirms the weak dependence of band offsets on surface orientation observed experimentally. We predict that the CBO of perfect CuInSe₂/CdS interface is large, 0.79 eV, which would dramatically degrade performance. Moreover we show that band gap widening induced by Ga adjusts only the VBO, and we find that Cd impurities do not significantly affect the CBO. Thus we show that Cu vacancies at the interface play the key role in enabling the tunability of CBO. We predict that Na further improves the CBO through electrostatically elevating the valence levels to decrease the CBO, explaining the observed essential role of Na for high performance. Moreover we find that K leads to a dramatic decrease in the CBO to 0.05 eV, much better than Na. We suggest that the efficiency of CIGS devices might be improved substantially by tuning the ratio of Na to K, with the improved phase stability of Na balancing phase instability from K. All these defects reduce interfacial stability slightly, but not significantly.

A number of exotic structures have been formed through high pressure chemistry, but applications have been hindered by difficulties in recovering the high pressure phase to ambient conditions (i.e., one atmosphere and room temperature). Here we use dispersion-corrected DFT (PBE-ulg flavor) to predict that above 60 GPa the most stable form of N₂O (the laughing gas in its molecular form) is a 1D polymer with an all-nitrogen backbone analogous to cis-polyacetylene in which alternate N are bonded (ionic covalent) to O. The analogous trans-polymer is only 0.03-0.10 eV/molecular unit less stable. Upon relaxation to ambient conditions both polymers relax below 14 GPa to the same stable non-planar trans-polymer, accompanied by possible electronic structure transitions. The predicted phonon spectrum and dissociation kinetics validate the stability of this trans-poly-NNO at ambient conditions, which has potential applications as a new type of conducting polymer with all-nitrogen chains and as a high-energy oxidizer for rocket propulsion. This work illustrates in silico materials discovery particularly in the realm of extreme conditions.

Modeling non-adiabatic electron dynamics has been a long-standing challenge for computational chemistry and materials science, and the eFF method presents a cost-efficient alternative. However, due to the deficiency of FSG representation, eFF is limited to low-Z elements with electrons of predominant s-character. To overcome this, we introduce a formal set of ECP extensions that enable accurate description of p-block elements. The extensions consist of a model representing the core electrons with the nucleus as a single pseudo particle represented by FSG, interacting with valence electrons through ECPs. We demonstrate and validate the ECP extensions for complex bonding structures, geometries, and energetics of systems with p-block character (C, O, Al, Si) and apply them to study materials under extreme mechanical loading conditions.

Despite its success, the eFF framework has some limitations, originated from both the design of Pauli potentials and the FSG representation. To overcome these, we develop a new framework of two-level hierarchy that is a more rigorous and accurate successor to the eFF method. The fundamental level, GHA-QM, is based on a new set of Pauli potentials that renders exact QM level of accuracy for any FSG represented electron systems. To achieve this, we start with using exactly derived energy expressions for the same spin electron pair, and fitting a simple functional form, inspired by DFT, against open singlet electron pair curves (H₂ systems). Symmetric and asymmetric scaling factors are then introduced at this level to recover the QM total energies of multiple electron pair systems from the sum of local interactions. To complement the imperfect FSG representation, the AMPERE extension is implemented, and aims at embedding the interactions associated with both the cusp condition and explicit nodal structures. The whole GHA-QM+AMPERE framework is tested on H element, and the preliminary results are promising.

Biological Activity of Rhodium Metalloinsertors and the Design of Bifunctional Conjugates

The Barton laboratory has established that octahedral rhodium complexes bearing the sterically expansive 5,6-chrysene diimine ligand can target thermodynamically destabilized sites, such as base pair mismatches, in DNA with high affinity and selectivity. These complexes approach DNA from the minor groove, ejecting the mismatched base pairs from the duplex in a binding mode termed metalloinsertion. In recent years, we have shown that these metalloinsertor complexes also exhibit cytotoxicity preferentially in cancer cells that are deficient in the mismatch repair (MMR) machinery.

Here, we establish that a sensitive structure-activity relationship exists for rhodium metalloinsertors. We studied the relationship between the chemical structures of metalloinsertors and their effect on biological activity for ten complexes with similar DNA binding affinities, but wide variation in their lipophilicity. Drastic differences were observed in the selectivities of the complexes for MMR-deficient cells. Compounds with hydrophilic ligands were highly selective, exhibiting preferential cytotoxicity in MMR-deficient cells at low concentrations and short incubation periods, whereas complexes with lipophilic ligands displayed poor cell-selectivity. It was discovered that all of the complexes localized to the nucleus in concentrations sufficient for mismatch binding; however, highly lipophilic complexes also exhibited high mitochondrial uptake. Significantly, these results support the notion that mitochondrial DNA is not the desired target for our metalloinsertor complexes; instead, selectivity stems from targeting mismatches in genomic DNA.

We have also explored the potential for metalloinsertors to be developed into more complex structures with multiple functionalities that could either enhance their overall potency or impart mismatch selectivity onto other therapeutic cargo. We have constructed a family of bifunctional metalloinsertor conjugates incorporating cis-platinum, each unique in its chemical structure, DNA binding interactions, and biological activity. The study of these complexes in MMR-deficient cells has established that the cell-selective biological activity of rhodium metalloinsertors proceeds through a critical cellular pathway leading to necrosis.

We further explored the underlying mechanisms surrounding the biological response to mismatch recognition by metalloinsertors in the genome. Immunofluorescence assays of MMR-deficient and MMR-proficient cells revealed that a critical biomarker for DNA damage, phosphorylation of histone H2AX (γH2AX) rapidly accumulates in response to metalloinsertor treatment, signifying the induction of double strand breaks in the genome. Significantly, we have discovered that our metalloinsertor complexes selectively inhibit transcription in MMR-deficient cells, which may be a crucial checkpoint in the eventual breakdown of the cell via necrosis. Additionally, preliminary in vivo studies have revealed the capability of these compounds to traverse the complex environments of multicellular organisms and accumulate in MMR-deficient tumors. Our ever-increasing understanding of metalloinsertors, as well as the development of new generations of complexes both monofunctional and bifunctional, enables their continued progress into the clinic as promising new chemotherapeutic agents.

A study of T=1, T=3/2, and T=2 states in some light nuclei using (He^3,n) reactions

The (He³, n) reactions on B¹¹, N¹⁵, O¹⁶, and O¹⁸ targets have been studied using a pulsed-beam time-of-flight spectrometer. Special emphasis was placed upon the determination of the excitation energies and properties of states with T = 1 (in Ne¹⁸), T = 3/2 (in N¹³ and F¹⁷) and T = 2 (in Ne²⁰). The identification of the T = 3/2 and T = 2 levels is based on the structure of these states as revealed by intensities and shapes of angular distributions. The reactions are interpreted in terms of double stripping theory. Angular distributions have been compared with plane and distorted wave stripping theories. Results for the four reactions are summarized below:

1) O¹⁶ (He³, n). The reaction has been studied at incident energies up to 13.5 MeV and two previously unreported levels in Ne¹⁸ were observed at E_x = 4.55 ± .015 MeV (Γ = 70 ± 30 keV) and E_x = 5.14 ± .018 MeV (Γ = 100 ± 40 keV).

2) B¹¹ (He³, n). The reaction has been studied at incident energies up to 13.5 MeV. Three T = 3/2 levels in N¹³ have been identified at E_x = 15.068 ± .008 MeV (Γ ˂ 15 keV), E_x = 18.44 ± .04, and E_x 18.98 ± .02 MeV (Γ = 40 ± 20 keV).

3) N¹⁵ (He³, n). The reaction has been studied at incident energies up to 11.88 MeV. T = 3/2 levels in F¹⁷ have been identified at E_x = 11.195 ± .007 MeV (Γ ˂ 20 keV), E_x = 12.540 ± .010 MeV (Γ ˂ 25 keV), and E_x = 13.095 ± .009 MeV (Γ ˂ 25 keV).

4) O¹⁸ (He³, n). The reaction has been studied at incident energies up to 9.0 MeV. The excitation energy of the lowest T = 2 level in Ne²⁰ has been found to be 16.730 ± .006 MeV (Γ ˂ 20 keV).

Angular distributions of the transitions leading to the above higher isospin states are well described by double stripping theory. Analog correspondences are established by comparing the present results with recent studies (t, p) and (He³, p) reactions on the same targets.

Beta and gamma vibrational bands in deformed nuclei

Energies and relative intensities of gamma transitions in ¹⁵²Sm, ¹⁵²Gd, ¹⁵⁴Gd, ¹⁶⁶Er, and ²³²U following radioactive decay have been measured with a Ge(Li) spectrometer. A peak fitting program has been developed to determine gamma ray energies and relative intensities with precision sufficient to give a meaningful test of nuclear models. Several previously unobserved gamma rays were placed in the nuclear level schemes. Particular attention has been paid to transitions from the beta and gamma vibrational bands, since the gamma ray branching ratios are sensitive tests of configuration mixing in the nuclear levels. As the reduced branching ratios depend on the multipolarity of the gamma transitions, experiments were performed to measure multipole mixing ratios for transitions from the gamma vibrational band. In ¹⁵⁴Gd, angular correlation experiments showed that transitions from the gamma band to the ground state band were predominantly electric quadrupole, in agreement with the rotational model. In ²³²U, the internal conversion spectrum has been studied with a Si(Li) spectrometer constructed for electron spectroscopy. The strength of electric monopole transitions and the multipolarity of some gamma transitions have been determined from the measured relative electron intensities.

The results of the experiments have been compared with the rotational model and several microscopic models. Relative B(E2) strengths for transitions from the gamma band in ²³²U and ¹⁶⁶Er are in good agreement with a single parameter band mixing model, with values of z₂= 0.025(10) and 0.046(2), respectively. Neither the beta nor the gamma band transition strengths in ¹⁵²Sm and ¹⁵⁴Gd can be accounted for by a single parameter theory, nor can agreement be found by considering the large mixing found between the beta and gamma bands. The relative B(E2) strength for transitions from the gamma band to the beta band in ²³²U is found to be five times greater than the strength to the ground state band, indicating collective transitions with strength approximately 15 single particle units.

1. Ultrasonic studies of binary liquid structure in the critical region. Theory and experiment for the 2,6-lutidine/water system. 2. Hartree-Fock calculations of electric polarizabilities of some simple atoms and molecules, and their practicality. 3. Calculation of vibrational transition probabilities in collinear atom-diatom and diatom-diatom collisions with Lennard-Jones interaction

Part 1. Many interesting visual and mechanical phenomena occur in the critical region of fluids, both for the gas-liquid and liquid-liquid transitions. The precise thermodynamic and transport behavior here has some broad consequences for the molecular theory of liquids. Previous studies in this laboratory on a liquid-liquid critical mixture via ultrasonics supported a basically classical analysis of fluid behavior by M. Fixman (e. g., the free energy is assumed analytic in intensive variables in the thermodynamics)--at least when the fluid is not too close to critical. A breakdown in classical concepts is evidenced close to critical, in some well-defined ways. We have studied herein a liquid-liquid critical system of complementary nature (possessing a lower critical mixing or consolute temperature) to all previous mixtures, to look for new qualitative critical behavior. We did not find such new behavior in the ultrasonic absorption ascribable to the critical fluctuations, but we did find extra absorption due to chemical processes (yet these are related to the mixing behavior generating the lower consolute point). We rederived, corrected, and extended Fixman's analysis to interpret our experimental results in these more complex circumstances. The entire account of theory and experiment is prefaced by an extensive introduction recounting the general status of liquid state theory. The introduction provides a context for our present work, and also points out problems deserving attention. Interest in these problems was stimulated by this work but also by work in Part 3.

Part 2. Among variational theories of electronic structure, the Hartree-Fock theory has proved particularly valuable for a practical understanding of such properties as chemical binding, electric multipole moments, and X-ray scattering intensity. It also provides the most tractable method of calculating first-order properties under external or internal one-electron perturbations, either developed explicitly in orders of perturbation theory or in the fully self-consistent method. The accuracy and consistency of first-order properties are poorer than those of zero-order properties, but this is most often due to the use of explicit approximations in solving the perturbed equations, or to inadequacy of the variational basis in size or composition. We have calculated the electric polarizabilities of H₂, He, Li, Be, LiH, and N₂ by Hartree-Fock theory, using exact perturbation theory or the fully self-consistent method, as dictated by convenience. By careful studies on total basis set composition, we obtained good approximations to limiting Hartree-Fock values of polarizabilities with bases of reasonable size. The values for all species, and for each direction in the molecular cases, are within 8% of experiment, or of best theoretical values in the absence of the former. Our results support the use of unadorned Hartree-Pock theory for static polarizabilities needed in interpreting electron-molecule scattering data, collision-induced light scattering experiments, and other phenomena involving experimentally inaccessible polarizabilities.

Part 3. Numerical integration of the close-coupled scattering equations has been carried out to obtain vibrational transition probabilities for some models of the electronically adiabatic H₂-H₂ collision. All the models use a Lennard-Jones interaction potential between nearest atoms in the collision partners. We have analyzed the results for some insight into the vibrational excitation process in its dependence on the energy of collision, the nature of the vibrational binding potential, and other factors. We conclude also that replacement of earlier, simpler models of the interaction potential by the Lennard-Jones form adds very little realism for all the complication it introduces. A brief introduction precedes the presentation of our work and places it in the context of attempts to understand the collisional activation process in chemical reactions as well as some other chemical dynamics.