A Cerenkov-ΔΕ-Cerenkov detector for high energy cosmic ray isotopes and an accelerator study of ^(40)Ar & ^(56)Fe fragmentation

This thesis has two major parts. The first part of the thesis will describe a high energy cosmic ray detector -- the High Energy Isotope Spectrometer Telescope (HEIST). HEIST is a large area (0.25 m²sr) balloon-borne isotope spectrometer designed to make high-resolution measurements of isotopes in the element range from neon to nickel (10 ≤ Z ≤ 28) at energies of about 2 GeV/nucleon. The instrument consists of a stack of 12 NaI(Tl) scintilla tors, two Cerenkov counters, and two plastic scintillators. Each of the 2-cm thick NaI disks is viewed by six 1.5-inch photomultipliers whose combined outputs measure the energy deposition in that layer. In addition, the six outputs from each disk are compared to determine the position at which incident nuclei traverse each layer to an accuracy of ~2 mm. The Cerenkov counters, which measure particle velocity, are each viewed by twelve 5-inch photomultipliers using light integration boxes.

HEIST-2 determines the mass of individual nuclei by measuring both the change in the Lorentz factor (Δγ) that results from traversing the NaI stack, and the energy loss (ΔΕ) in the stack. Since the total energy of an isotope is given by Ε = γM, the mass M can be determined by M = ΔΕ/Δγ. The instrument is designed to achieve a typical mass resolution of 0.2 amu.

The second part of this thesis presents an experimental measurement of the isotopic composition of the fragments from the breakup of high energy ⁴⁰Ar and ⁵⁶Fe nuclei. Cosmic ray composition studies rely heavily on semi-empirical estimates of the cross-sections for the nuclear fragmentation reactions which alter the composition during propagation through the interstellar medium. Experimentally measured yields of isotopes from the fragmentation of ⁴⁰Ar and ⁵⁶Fe are compared with calculated yields based on semi-empirical cross-section formulae. There are two sets of measurements. The first set of measurements, made at the Lawrence Berkeley Laboratory Bevalac using a beam of 287 MeV/nucleon ⁴⁰Ar incident on a CH₂ target, achieves excellent mass resolution (σ_m ≤ 0.2 amu) for isotopes of Mg through K using a Si(Li) detector telescope. The second set of measurements, also made at the Lawrence Berkeley Laboratory Bevalac, using a beam of 583 MeV/nucleon ⁵⁶FeFe incident on a CH₂ target, resolved Cr, Mn, and Fe fragments with a typical mass resolution of ~ 0.25 amu, through the use of the Heavy Isotope Spectrometer Telescope (HIST) which was later carried into space on ISEE-3 in 1978. The general agreement between calculation and experiment is good, but some significant differences are reported here.

Discrete Modeling of Granular Media: A NURBS-based Approach

This dissertation is concerned with the development of a new discrete element method (DEM) based on Non-Uniform Rational Basis Splines (NURBS). With NURBS, the new DEM is able to capture sphericity and angularity, the two particle morphological measures used in characterizing real grain geometries. By taking advantage of the parametric nature of NURBS, the Lipschitzian dividing rectangle (DIRECT) global optimization procedure is employed as a solution procedure to the closest-point projection problem, which enables the contact treatment of non-convex particles. A contact dynamics (CD) approach to the NURBS-based discrete method is also formulated. By combining particle shape flexibility, properties of implicit time-integration, and non-penetrating constraints, we target applications in which the classical DEM either performs poorly or simply fails, i.e., in granular systems composed of rigid or highly stiff angular particles and subjected to quasistatic or dynamic flow conditions. The CD implementation is made simple by adopting a variational framework, which enables the resulting discrete problem to be readily solved using off-the-shelf mathematical programming solvers. The capabilities of the NURBS-based DEM are demonstrated through 2D numerical examples that highlight the effects of particle morphology on the macroscopic response of granular assemblies under quasistatic and dynamic flow conditions, and a 3D characterization of material response in the shear band of a real triaxial specimen.

Source characteristics of recent and historic earthquakes in Central and Southern California: results from forward modeling

The long- and short-period body waves of a number of moderate earthquakes occurring in central and southern California recorded at regional (200-1400 km) and teleseismic (> 30°) distances are modeled to obtain the source parameters-focal mechanism, depth, seismic moment, and source time history. The modeling is done in the time domain using a forward modeling technique based on ray summation. A simple layer over a half space velocity model is used with additional layers being added if necessary-for example, in a basin with a low velocity lid.

The earthquakes studied fall into two geographic regions: 1) the western Transverse Ranges, and 2) the western Imperial Valley. Earthquakes in the western Transverse Ranges include the 1987 Whittier Narrows earthquake, several offshore earthquakes that occurred between 1969 and 1981, and aftershocks to the 1983 Coalinga earthquake (these actually occurred north of the Transverse Ranges but share many characteristics with those that occurred there). These earthquakes are predominantly thrust faulting events with the average strike being east-west, but with many variations. Of the six earthquakes which had sufficient short-period data to accurately determine the source time history, five were complex events. That is, they could not be modeled as a simple point source, but consisted of two or more subevents. The subevents of the Whittier Narrows earthquake had different focal mechanisms. In the other cases, the subevents appear to be the same, but small variations could not be ruled out.

The recent Imperial Valley earthquakes modeled include the two 1987 Superstition Hills earthquakes and the 1969 Coyote Mountain earthquake. All are strike-slip events, and the second 1987 earthquake is a complex event With non-identical subevents.

In all the earthquakes studied, and particularly the thrust events, constraining the source parameters required modeling several phases and distance ranges. Teleseismic P waves could provide only approximate solutions. P_(nl) waves were probably the most useful phase in determining the focal mechanism, with additional constraints supplied by the SH waves when available. Contamination of the SH waves by shear-coupled PL waves was a frequent problem. Short-period data were needed to obtain the source time function.

In addition to the earthquakes mentioned above, several historic earthquakes were also studied. Earthquakes that occurred before the existence of dense local and worldwide networks are difficult to model due to the sparse data set. It has been noticed that earthquakes that occur near each other often produce similar waveforms implying similar source parameters. By comparing recent well studied earthquakes to historic earthquakes in the same region, better constraints can be placed on the source parameters of the historic events.

The Lompoc earthquake (M=7) of 1927 is the largest offshore earthquake to occur in California this century. By direct comparison of waveforms and amplitudes with the Coalinga and Santa Lucia Banks earthquakes, the focal mechanism (thrust faulting on a northwest striking fault) and long-period seismic moment (10^(26) dyne cm) can be obtained. The S-P travel times are consistent with an offshore location, rather than one in the Hosgri fault zone.

Historic earthquakes in the western Imperial Valley were also studied. These events include the 1942 and 1954 earthquakes. The earthquakes were relocated by comparing S-P and R-S times to recent earthquakes. It was found that only minor changes in the epicenters were required but that the Coyote Mountain earthquake may have been more severely mislocated. The waveforms as expected indicated that all the events were strike-slip. Moment estimates were obtained by comparing the amplitudes of recent and historic events at stations which recorded both. The 1942 event was smaller than the 1968 Borrego Mountain earthquake although some previous studies suggested the reverse. The 1954 and 1937 earthquakes had moments close to the expected value. An aftershock of the 1942 earthquake appears to be larger than previously thought.

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.

Characterization and modeling of premixed turbulent n-heptane flames in the thin reaction zone regime

n-heptane/air premixed turbulent flames in the high-Karlovitz portion of the thin reaction zone regime are characterized and modeled in this thesis using Direct Numerical Simulations (DNS) with detailed chemistry. In order to perform these simulations, a time-integration scheme that can efficiently handle the stiffness of the equations solved is developed first. A first simulation with unity Lewis number is considered in order to assess the effect of turbulence on the flame in the absence of differential diffusion. A second simulation with non-unity Lewis numbers is considered to study how turbulence affects differential diffusion. In the absence of differential diffusion, minimal departure from the 1D unstretched flame structure (species vs. temperature profiles) is observed. In the non-unity Lewis number case, the flame structure lies between that of 1D unstretched flames with "laminar" non-unity Lewis numbers and unity Lewis number. This is attributed to effective Lewis numbers resulting from intense turbulent mixing and a first model is proposed. The reaction zone is shown to be thin for both flames, yet large chemical source term fluctuations are observed. The fuel consumption rate is found to be only weakly correlated with stretch, although local extinctions in the non-unity Lewis number case are well correlated with high curvature. These results explain the apparent turbulent flame speeds. Other variables that better correlate with this fuel burning rate are identified through a coordinate transformation. It is shown that the unity Lewis number turbulent flames can be accurately described by a set of 1D (in progress variable space) flamelet equations parameterized by the dissipation rate of the progress variable. In the non-unity Lewis number flames, the flamelet equations suggest a dependence on a second parameter, the diffusion of the progress variable. A new tabulation approach is proposed for the simulation of such flames with these dimensionally-reduced manifolds.

LiNbO3：Fe晶体光栅擦除特性

在光折变过程中，光波对掺杂中心的激发系数S是一个重要的物理量。为了求出光激发系数S对波长A的依赖，研究了在不同的擦除波长条件下，记录在光折变晶体LiNbO3：Fe中的光栅衰减特性。实验结果表明lnP（P为光擦除灵敏度）与擦除波长λ呈线性关系，即波长越短，光栅擦除越快，光擦除灵敏度越高。在单中心电子带传输模型的基础上，理论计算也证实了实验结果。利用带传输模型，给出了理论上光激发系数S对波长的指数依赖关系。同时微观量φμγ（φ为量子效率）也可由光栅擦除实验求出。

组合外加电场提高Fe：LiNbO3中90°记录体全息的衍射效率

针对90°全息记录结构的特殊性，在记录和读出阶段分别施加不同极性的外加电场，实现了Fe：LiNbO3晶体中高衍射效率的全息记录和读出。在库赫塔列夫（Kukhtarev）方程基础上引人两维耦合波理论，对组合外加电场提高掺铁铌酸锂晶体光折变特性的机理进行了探讨。实验和理论计算结果都表明这种分别在记录和读出过程施加不同极性组合的外加电场是在Fe：LiNbO3中实现高衍射效率90°记录结构体全息的有效技术方案。

LiNbO3：Fe：Ni晶体非挥发全息存储研究

采用三种不同的双光记录方案进行了LiNbO3：Fe：Ni晶体全息存储实验．详细研究了饱和衍射效率、固定衍射效率、动态范围和记录灵敏度，以及退火条件对记录的影响。结果表明，氧化LiNbO3：Fe：Ni晶体的饱和衍射效率、固定衍射效率和记录灵敏度比其他报道的双掺杂LiNbO3晶体高。结合掺杂能级图，理论分析了LiNbO3双掺杂晶体深陷阱中心能级的相对位置及其微观光学参量对全息记录性能的影响。LiNbO3：Fe：Ni晶体有望成为一种新的高效率非挥发全息存储材料。

双掺杂LiNbO3：Fe：Ru晶体的全息读出特性

为了解决读出过程中全息数据的擦除问题，研究了新型双掺杂LiNbOa：Fe：Ru晶体的全息读出特性。分析了双中心和单中心记录方案光栅的读出特性，并联立双中心物质方程和耦合波方程进行了模拟计算。结果表明，双中心记录所得到的饱和全息的读出时间常数远低于LiNbO2：Fe：Mn晶体的读出时间常数；单色光记录可以实现有效的全息，且其读出时间常数远大于记录时间常数，表现为准态非挥发读出。分析表明，这可能由于Ru的能级比Mn更靠近Fe，更易被红光激发，从而使得双中心记录所得饱和光栅的存贮持久性降低；单色光记录中红光能够

LiNbO3：Ru：Fe晶体畴反转过程中的电色效应

观察到了双掺杂LiNbO3：Ru：Fe晶体中电色效应伴随畴反转而发生，且与畴反转一样也具有可逆性，两者相辅相成，畴反转导致了晶体变色，电色效应促进了畴反转，系统的实验结果证明了两者的相辅相成性．基于铌酸锂铁电微结构模型，简要解释了其机理．而且发现在极化过程中电色效应促使了畴核的形成，使之不同于非掺杂同成分比铌酸锂晶体的矫顽场大于击穿电场，用恒定直流电场代替脉冲电场也能实现畴反转，这将为周期性极化铌酸锂的制备提供一种新的技术改进．

LiNbO3∶Fe∶Rh晶体近红外非挥发全息记录

采用双中心记录方案在双掺杂LiNbO3∶Fe∶Rh晶体中实现了近红外非挥发全息记录,研究了LiNbO3∶Fe∶Rh晶体在633 nm,752 nm,799 nm波长下的全息记录性能。结果表明,在使用近红外记录光时,其记录灵敏度随敏化光强的变化趋势与双中心短波长记录时的不同。通过和LiNbO3∶Fe∶Mn等传统双掺杂铌酸锂晶体的近红外波段记录效果对比,发现同时掺杂Fe和Rh可增强晶体对近红外光的吸收,获得更高的浅中心Fe光生伏特系数,从而能够在LiNbO3∶Fe∶Rh晶体中实现近红外波段的光折变全息记录。

Quantificação de metais (Al, Fe, Mn, Pb, Zn e Cu) nos sedimentos do sistema estuarino das Ilhas de Tinharé-Boipeba (Bahia - Brasil)

O sistema estuarino das Ilhas de Tinharé-Boipeba está inserido na região do Baixo Sul baiano (Bahia Brasil), costa nordeste brasileira, e conectam-se com o oceano através de três saídas principais, Morro de São Paulo, Boipeba e Barra dos Carvalhos. Informações relacionadas à presença de metais no sedimento são quase inexistentes para a região. Nestes estuários foram coletadas 40 amostras de sedimento, onde analisou-se a concentração de metais (Al, Cu, Fe, Mn, Pb e Zn) através da extração com água-régia, segundo o protocolo do material de referência BCR-701 (RAURET et al., 2001). Baixas concentrações de metais foram registradas nas proximidades das saídas para o mar e no Canal de Garapuá e Rio Grande. Altas concentrações, porém praticamente dentro dos valores de referência na CONAMA 344/04, foram registradas no Rio Una, na maior parte do Rio Cairu, Rio das Almas, em quase todas as estações, no braço leste do Rio Cairu e na porção intermediária do Rio dos Patos. As concentrações dos metais (mg.kg1) apresentaram valores entre os seguintes intervalos, Al (3,57 x 104 a 3,19 x 102), Cu (1,02 x 102 a 0,35), Fe (4,33 x 104 a 2,05 x 102), Mn (1,44 x 103 a 2,73), Pb (6,67 x 101 a 2,66) e Zn (5,08 x 102 a 3,18). Através da análise estatística ACP (Análise dos Componentes Principais) e dos gráficos de correlação entre a granulometria do sedimento (areia, silte e argila) e os metais, observou-se maiores concentrações de metais com o aumento do percentual de silte e diminuição do percentual de areia. Também se identificou uma forte correlação entre a ocorrência do fósforo e a presença de metais. Acredita-se que a principal espécie química em questão seja o fosfato (PO43 ou P2O5) que é uma base dura, onde a ligação com os metais se dá pelo oxigênio, sendo o caráter iônico relevante. Não foi identificada uma correlação entre a presença de metais e a argila, fato atribuído ao baixo teor dessa granulometria para todas as amostras do estudo e a composição da argila desses estuários

A study of the proposed Santa Fe extension from the present Redondo Branch at El Segundo to Los Angeles Harbor

The Los Angeles Harbor at San Pedro with its natural advantages, and the big development of these now underway, will very soon be the key to the traffic routes of Southern California. The Atchison, Topeka, and Santa Fe railway company realizing this and, not wishing to be caught asleep, has planned to build a line from El Segundo to the harbor. The developments of the harbor are not the only developments taking place in these localities and the proposed new line is intended to include these also.

Nonvolatile holograms in LiNbO3 : Fe : Cu by use of the bleaching effect

We report our observation of a bleaching effect under an ultraviolet exposure in LiNbO3:Fe:Cu crystals. Two three-step recording-transferring-fixing schemes are proposed to record nonvolatile photorefractive holograms in such crystals. In the schemes two red laser beams and an ultraviolet illumination are used selectively to write the charge grating in the shallow-level Fe centers, to develop the charge grating in the deep-level Cu centers by transferring the charge grating in the Fe centers, and to fix only the charge grating in the Cu centers for unerasable read-out. Experimental results, verifications, and an optimal recording scheme are given. A comparison of the lithium niobate crystals of the same double-doping system of Fe:Mn, Ce:Mn, Ce:Cu, and Fe:Cu is outlined. (C) 2002 Optical Society of America.