Metasomatism and magmatic assimilation at a gabbro-limestone contact, Christmas Mountains, Big Bend region, Texas

A composite stock of alkaline gabbro and syenite is intrusive into limestone of the Del Carmen, Sue Peake and Santa Elena Formations at the northwest end of the Christmas Mountains. There is abundant evidence of solution of wallrock by magma but nowhere are gabbro and limestone in direct contact. The sequence of lithologies developed across the intrusive contact and across xenoliths is gabbro, pyroxenite, calc-silicate skarn, marble. Pyroxenite is made up of euhedral crystals of titanaugite and sphene in a leucocratic matrix of nepheline, Wollastonite and alkali feldspar. The uneven modal distribution of phases in pyroxenite and the occurrence' of nepheline syenite dikes, intrusive into pyroxenite and skarn, suggest that pyroxenite represents an accumulation of clinopyroxene "cemented" together by late-solidifying residual magma of nepheline syenite composition. Assimilation of limestone by gabbroic magma involves reactions between calcite and magma and/or crystals in equilibrium with magma and crystallization of phases in which the magma is saturated, to supply energy for the solution reaction. Gabbroic magma was saturated with plagioclase and clinopyroxene at the time of emplacement. The textural and mineralogic features of pyroxenite can be produced by the reaction 2( 1-X) CALCITE + AN_XAB_l-X = (1-X) NEPHELINE+ 2(1-X) WOLLASTONITE+ X ANORTHITE+ 2(1-X) CO₂. Plagioclase in pyroxenite has corroded margins and is rimmed by nepheline, suggestive of resorption by magma. Anorthite and wollastonite enter solid solution in titanaugite. For each mole of calcite dissolved, approximately one mole of clinopyroxene was crystallized. Thus the amount of limestone that may be assimilated is limited by the concentration of potential clinopyroxene in the magma. Wollastonite appears as a phase when magma has been depleted in iron and magnesium by crystallization of titanaugite. The predominance of mafic and ultramafic compositions among contaminated rocks and their restriction to a narrow zone along the intrusive contact provides little evidence for the generation of a significant volume of desilicated magma as a result of limestone assimilation.

Within 60 m of the intrusive contact with the gabbro, nodular chert in the Santa Elena Limestone reacted with the enveloping marble to form spherical nodules of high-temperature calc-silicate minerals. The phases wollastonite, rankinite, spurrite, tilleyite and calcite, form a series of sharply-bounded, concentric monomineralic and two-phase shells which record a step-wise decrease in silica content from the core of a nodule to its rim. Mineral zones in the nodules vary 'with distance from the gabbro as follows:

0-5 m CALCITE + SPURRITE + RANKINITE + WOLLASTONITE
5-16 m CALCITE + TILLEYITE ± SPURRITE + RANKINITE + WOLLASTONITE
16-31 m CALCITE + TILLEYITE + WOLLASTONITE
31-60 m CALCITE + WOLLASTONITE
60-plus CALCITE + QUARTZ

The mineral of a one-phase zone is compatible with the phases bounding it on either side but these phases are incompatible in the same volume of P-T-X_CO2.

Growth of a monomineralio zone is initiated by reaction between minerals of adjacent one-phase zones which become unstable with rising temperature to form a thin layer of a new single phase that separates the reactants and is compatible with both of them. Because the mineral of the new zone is in equilibrium with the phases at both of its contacts, gradients in the chemical potentials of the exchangeable components are established across it. Although zone boundaries mark discontinuities in the gradients of bulk composition, two-phase equilibria at the contacts demonstrate that the chemical potentials are continuous. Hence, Ca, Si and CO₂ were redistributed in the growing nodule by diffusion. A monomineralic zone grows at the expense of an adjacent zone by reaction between diffusing components and the mineral of the adjacent zone. Equilibria between two phases at zone boundaries buffers the chemical potentials of the diffusing species. Thus, within a monomineralic zone, the chemical potentials of the diffusing components are controlled external to the local assemblage by the two-phase equilibria at the zone boundaries.

Mineralogically zoned calc-silicate skarn occurs as a narrow band that separates pyroxenite and marble along the intrusive contact and forms a rim on marble xenoliths in gabbro. Skarn consists of melilite or idocrase pseudomorphs of melili te, one or two . stoichiometric calcsilicate phases and accessory Ti-Zr garnet, perovskite and magnetite. The sequence of mineral zones from pyroxenite to marble, defined by a characteristic calc-silicate, is wollastonite, rankinite, spurrite, calcite. Mineral assemblages of adjacent skarn zones are compatible and the set of zones in a skarn band defines a facies type, indicating that the different mineral assemblages represent different bulk compositions recrystallized under identical conditions. The number of phases in each zone is less than the number that might be expected to result from metamorphism of a general bulk composition under conditions of equilibrium, trivariant in P, T and u_CO2. The "special" bulk composition of each zone is controlled by reaction between phases of the zones bounding it on either side. The continuity of the gradients of composition of melilite and garnet solid solutions across the skarn is consistent with the local equilibrium hypothesis and verifies that diffusion was the mechanism of mass transport. The formula proportions of Ti and Zr in garnet from skarn vary antithetically with that of Si Which systematically decreases from pyroxenite to marble. The chemical potential of Si in each skarn zone was controlled by the coexisting stoichiometric calc-silicate phases in the assemblage. Thus the formula proportion of Si in garnet is a direct measure of the chemical potential of Si from point to point in skarn. Reaction between gabbroic magma saturated with plagioclase and clinopyroxene produced nepheline pyroxenite and melilite-wollastonite skarn. The calcsilicate zones result from reaction between calcite and wollastonite to form spurrite and rankinite.

Broad-bandwidth semi-noncollinear optical parametric amplification in periodically poled LiNbO3 based on tilted quasi-phase-matched gratings

On the basis of noncollinear optical parametric amplification in periodically poled lithium niobate (PPLN) which is realized by quasi-phase matching (QPM) technology, we consider the possibility of semi-noncollinear phase matching between collinear and noncollinear geometries by tilting a PPLN-crystal's parallel grating at a sure angle. Numerical simulation with proper parameters shows that we can achieve a broader optical parametric amplification (OPA) bandwidth than that of noncollinear geometry. About 121 nm at a signal wavelength of 800 and 70 nm at a signal wavelength of 1064 nm under optimal conditions are obtained when the crystal length is 9 mm.

Selective metallization on insulator surfaces with femtosecond laser pulses

We report selective metallization on surfaces of insulators ( glass slides and lithium niobate crystal) based on femtosecond laser modification combined with electroless plating. The process is mainly composed of four steps: (1) formation of silver nitrate thin films on the surfaces of glass or crystal substrates; (2) generation of silver particles in the irradiated area by femtosecond laser direct writing; (3) removal of unirradiated silver nitrate films; and (4) selective electroless plating in the modified area. We discuss the mechanism of selective metallization on the insulators. Moreover, we investigate the electrical and adhesive properties of the copper microstructures patterned on the insulator surfaces, showing great potential of integrating electrical functions into lab-on-a-chip devices. (C) 2007 Optical Society of America.

周期性极化铌酸锂晶体中半非共线型宽带光学参变放大理论研究

基于周期性极化铌酸锂晶体(PPLN)的准相位匹配光参变放大过程,通过倾斜周期极化铌酸锂晶体中极化域(极化光栅)一定角度,实现了介于共线匹配方式和非共线匹配方式之间的一种半非共线型准相位匹配方式,并以该匹配方式下的各光矢量几何关系得出相位匹配曲线,找到在特定抽运光和信号光波长下能获得宽带增益放大的周期极化长度。并研究其极化倾斜角度与温度特性。模拟计算表明,在合适的角度与温度条件下,该方式可以532 nm抽运光抽运的信号光在800 nm和1064 nm处均获得宽带光参变放大。

Surface effects on spinwave resonance in thin magnetic films

Over the past few decades, ferromagnetic spinwave resonance in magnetic thin films has been used as a tool for studying the properties of magnetic materials. A full understanding of the boundary conditions at the surface of the magnetic material is extremely important. Such an understanding has been the general objective of this thesis. The approach has been to investigate various hypotheses of the surface condition and to compare the results of these models with experimental data. The conclusion is that the boundary conditions are largely due to thin surface regions with magnetic properties different from the bulk. In the calculations these regions were usually approximated by uniform surface layers; the spins were otherwise unconstrained except by the same mechanisms that exist in the bulk (i.e., no special "pinning" at the surface atomic layer is assumed). The variation of the ferromagnetic spinwave resonance spectra in YIG films with frequency, temperature, annealing, and orientation of applied field provided an excellent experimental basis for the study.

This thesis can be divided into two parts. The first part is ferromagnetic resonance theory; the second part is the comparison of calculated with experimental data in YIG films. Both are essential in understanding the conclusion that surface regions with properties different from the bulk are responsible for the resonance phenomena associated with boundary conditions.

The theoretical calculations have been made by finding the wave vectors characteristic of the magnetic fields inside the magnetic medium, and then combining the fields associated with these wave vectors in superposition to match the specified boundary conditions. In addition to magnetic boundary conditions required for the surface layer model, two phenomenological magnetic boundary conditions are discussed in detail. The wave vectors are easily found by combining the Landau-Lifshitz equations with Maxwell's equations. Mode positions are most easily predicted from the magnetic wave vectors obtained by neglecting damping, conductivity, and the displacement current. For an insulator where the driving field is nearly uniform throughout the sample, these approximations permit a simple yet accurate calculation of the mode intensities. For metal films this calculation may be inaccurate but the mode positions are still accurately described. The techniques necessary for calculating the power absorbed by the film under a specific excitation including the effects of conductivity, displacement current and damping are also presented.

In the second part of the thesis the properties of magnetic garnet materials are summarized and the properties believed associated with the two surface regions of a YIG film are presented. Finally, the experimental data and calculated data for the surface layer model and other proposed models are compared. The conclusion of this study is that the remarkable variety of spinwave spectra that arises from various preparation techniques and subsequent treatments can be explained by surface regions with magnetic properties different from the bulk.

Ferromagnetic Resonance Force Microscopy

A new approach to magnetic resonance was introduced in 1992 based upon detection of spin-induced forces by J. Sidles [1]. This technique, now called magnetic resonance force microscopy (MRFM), was first demonstrated that same year via electron paramagnetic resonance (EPR) by D. Rugar et al. [2]. This new method combines principles of magnetic resonance with those of scanned probe technology to detect spin resonance through mechanical, rather than inductive, means. In this thesis the development and use of ferromagnetic resonance force microscopy (FMRFM) is described. This variant of MRFM, which allows investigation of ferromagnetic samples, was first demonstrated in 1996 by Z. Zhang et al. [3]. FMRFM enables characterization of (a) the dynamic magnetic properties of microscale magnetic devices, and (b) the spatial dependence of ferromagnetic resonance within a sample. Both are impossible with conventional ferromagnetic resonance techniques.

Ferromagnetically coupled systems, however, pose unique challenges for force detection. In this thesis the attainable spatial resolution - and the underlying physical mechanisms that determine it - are established. We analyze the dependence of the magnetostatic modes upon sample dimensions using a series of microscale yttrium iron garnet (YIG) samples. Mapping of mode amplitudes within these sample is attained with an unprecedented spatial resolution of 15μm. The modes, never before analyzed on this scale, fit simple models developed in this thesis for samples of micron dimensions. The application of stronger gradient fields induces localized perturbation of the ferromagnetic resonance modes. The first demonstrations of this effect are presented in this study, and a simple theoretical model is developed to explain our observations. The results indicate that the characteristics of the locally-detected ferromagnetic modes are still largely determined by the external fields and dimensions of the entire sample, rather than by the localized interaction volume (i.e., the locale most strongly affected by the local gradient field). Establishing this is a crucial first step toward understanding FMRFM in the high gradient field limit where the dispersion relations become locally determined. In this high gradient field regime, FMRFM imaging becomes analogous with that of EPR MRFM.

FMRFM has also been employed to characterize magnetic multilayers, similar to those utilized in giant magnetoresistance (GMR) devices, on a lateral scale 40 x 40μm. This is orders of magnitude smaller than possible via conventional methods. Anisotropy energies, thickness, and interface qualities of individual layers have been resolved.

This initial work clearly demonstrates the immense and unique potential that FMRFM offers for characterizing advanced magnetic nanostructures and magnetic devices.

Coherent scattering of light by Bose-Einstein condensation in a tight confinement

We have observed strong scattering of a probe light by dilute Bose-Einstein condensate (BEC) Rb-87 gas in a tight magnetic trap. The scattering light forms fringes at the image plane. It is found that we can infer the real size of the condensation and the number of the atoms by modelling the imaging system. We present a quantitative calculation of light scattering by the condensed atoms. The calculation shows that the experimental results agree well with the prediction of the generalized diffraction theory, and thus we can directly observe the phase transition of BEC in a tight trap.

Oxygen isotope ratios in coexisting minerals of regionally metamorphosed rocks

The O¹⁸/O¹⁶ ratios of coexisting minerals from a number of regionally metamorphosed rocks have been measured, using a bromine pentafluoride extraction-technique. Listed in order of their increasing tendency to concentrate O¹⁸, the minerals analyzed are magnetite, ilmenite, chlorite, biotite, garnet, hornblende, kyanite, muscovite, feldspar, and quartz. The only anomalous sequence detected occurs in a xenolith of schist, in which quartz, muscovite, biotite, and ilmenite, but not garnet, have undergone isotopic exchange with surrounding trondjemite.

With few exceptions, quartz-magnetite and quartz-ilmenite fractionations decrease with increasing metamorphic grade determined by mineral paragenesis and spatial distribution. This consistency does not apply to quartz-magnetite and quartz-ilmenite fractionations obtained from rocks in which petrographic evidence of retrogradation is present.

Whereas measured isotopic fractionations among quartz, garnet, ilmenite, and magnetite are approximately related to metamorphic grade, fractionations between these minerals and biotite or muscovite show poor correlation with grade. Variations in muscovite-biotite fractionations are relatively small. These observations are interpreted to mean that muscovite and biotite are affected by retrograde re-equilibration to a greater extent than the anhydrous minerals analyzed.

Measured quartz-ilmenite fractionations range from 12 permil in the biotite zone of central Vermont to 6.5 permil in the sillimanite-orthoclase zone of southeastern Connecticut. Analyses of natural assemblages from the kyanite and sillimanite zones suggest that equilibrium quartz-ilmenite fractionations are approximately 8 percent smaller than corresponding quartz-magnetite fractionations. Employing the quartz-magnetite geothermometer calibrated by O'Neil and Clayton (1964), a temperature of 560°C was obtained for kyanite-bearing schists from Addison County, Vermont. Extending the calibration to quartz-ilmenite fractionations, a temperature of 600°C was obtained for kyanite-schists from Shoshone County, Idaho. At these temperatures kyanite is stable only at pressures exceeding 11 kbars (Bell, 1963), corresponding to lithostatic loads of over 40 km.

The Equation of State and Petrogenesis of Komatiite

(1) Equation of State of Komatiite

The equation of state (EOS) of a molten komatiite (27 wt% MgO) was detennined in the 5 to 36 GPa pressure range via shock wave compression from 1550°C and 0 bar. Shock wave velocity, U_S, and particle velocity, U_P, in km/s follow the linear relationship U_S = 3.13(±0.03) + 1.47(±0.03) U_P. Based on a calculated density at 1550°C, 0 bar of 2.745±0.005 glee, this U_S-U_P relationship gives the isentropic bulk modulus K_S = 27.0 ± 0.6 GPa, and its first and second isentropic pressure derivatives, K'_S = 4.9 ± 0.1 and K"_S = -0.109 ± 0.003 GPa^-1.

The calculated liquidus compression curve agrees within error with the static compression results of Agee and Walker [1988a] to 6 GPa. We detennine that olivine (FO₉₄) will be neutrally buoyant in komatiitic melt of the composition we studied near 8.2 GPa. Clinopyroxene would also be neutrally buoyant near this pressure. Liquidus garnet-majorite may be less dense than this komatiitic liquid in the 20-24 GPa interval, however pyropic-garnet and perovskite phases are denser than this komatiitic liquid in their respective liquidus pressure intervals to 36 GPa. Liquidus perovskite may be neutrally buoyant near 70 GPa.

At 40 GPa, the density of shock-compressed molten komatiite would be approximately equal to the calculated density of an equivalent mixture of dense solid oxide components. This observation supports the model of Rigden et al. [1989] for compressibilities of liquid oxide components. Using their theoretical EOS for liquid forsterite and fayalite, we calculate the densities of a spectrum of melts from basaltic through peridotitic that are related to the experimentally studied komatiitic liquid by addition or subtraction of olivine. At low pressure, olivine fractionation lowers the density of basic magmas, but above 14 GPa this trend is reversed. All of these basic to ultrabasic liquids are predicted to have similar densities at 14 GPa, and this density is approximately equal to the bulk (PREM) mantle. This suggests that melts derived from a peridotitic mantle may be inhibited from ascending from depths greater than 400 km.

The EOS of ultrabasic magmas was used to model adiabatic melting in a peridotitic mantle. If komatiites are formed by >15% partial melting of a peridotitic mantle, then komatiites generated by adiabatic melting come from source regions in the lower transition zone (≈500-670 km) or the lower mantle (>670 km). The great depth of incipient melting implied by this model, and the melt density constraint mentioned above, suggest that komatiitic volcanism may be gravitationally hindered. Although komatiitic magmas are thought to separate from their coexisting crystals at a temperature =200°C greater than that for modern MORBs, their ultimate sources are predicted to be diapirs that, if adiabatically decompressed from initially solid mantle, were more than 700°C hotter than the sources of MORBs and derived from great depth.

We considered the evolution of an initially molten mantle, i.e., a magma ocean. Our model considers the thermal structure of the magma ocean, density constraints on crystal segregation, and approximate phase relationships for a nominally chondritic mantle. Crystallization will begin at the core-mantle boundary. Perovskite buoyancy at > 70 GPa may lead to a compositionally stratified lower mantle with iron-enriched mangesiowiistite content increasing with depth. The upper mantle may be depleted in perovskite components. Olivine neutral buoyancy may lead to the formation of a dunite septum in the upper mantle, partitioning the ocean into upper and lower reservoirs, but this septum must be permeable.

(2) Viscosity Measurement with Shock Waves

We have examined in detail the analytical method for measuring shear viscosity from the decay of perturbations on a corrugated shock front The relevance of initial conditions, finite shock amplitude, bulk viscosity, and the sensitivity of the measurements to the shock boundary conditions are discussed. The validity of the viscous perturbation approach is examined by numerically solving the second-order Navier-Stokes equations. These numerical experiments indicate that shock instabilities may occur even when the Kontorovich-D'yakov stability criteria are satisfied. The experimental results for water at 15 GPa are discussed, and it is suggested that the large effective viscosity determined by this method may reflect the existence of ice VII on the Rayleigh path of the Hugoniot This interpretation reconciles the experimental results with estimates and measurements obtained by other means, and is consistent with the relationship of the Hugoniot with the phase diagram for water. Sound waves are generated at 4.8 MHz at in the water experiments at 15 GPa. The existence of anelastic absorption modes near this frequency would also lead to large effective viscosity estimates.

(3) Equation of State of Molybdenum at 1400°C

Shock compression data to 96 GPa for pure molybdenum, initially heated to 1400°C, are presented. Finite strain analysis of the data gives a bulk modulus at 1400°C, K'_S. of 244±2 GPa and its pressure derivative, K'_OS of 4. A fit of shock velocity to particle velocity gives the coefficients of U_S = C_O+S U_P to be C_O = 4.77±0.06 km/s and S = 1.43±0.05. From the zero pressure sound speed, C_O, a bulk modulus of 232±6 GPa is calculated that is consistent with extrapolation of ultrasonic elasticity measurements. The temperature derivative of the bulk modulus at zero pressure, θK_OSθT|_P, is approximately -0.012 GPa/K. A thermodynamic model is used to show that the thermodynamic Grüneisen parameter is proportional to the density and independent of temperature. The Mie-Grüneisen equation of state adequately describes the high temperature behavior of molybdenum under the present range of shock loading conditions.

Development and Application of Embedding Methods for the Simulation of Large Chemical Systems

The high computational cost of correlated wavefunction theory (WFT) calculations has motivated the development of numerous methods to partition the description of large chemical systems into smaller subsystem calculations. For example, WFT-in-DFT embedding methods facilitate the partitioning of a system into two subsystems: a subsystem A that is treated using an accurate WFT method, and a subsystem B that is treated using a more efficient Kohn-Sham density functional theory (KS-DFT) method. Representation of the interactions between subsystems is non-trivial, and often requires the use of approximate kinetic energy functionals or computationally challenging optimized effective potential calculations; however, it has recently been shown that these challenges can be eliminated through the use of a projection operator. This dissertation describes the development and application of embedding methods that enable accurate and efficient calculation of the properties of large chemical systems.

Chapter 1 introduces a method for efficiently performing projection-based WFT-in-DFT embedding calculations on large systems. This is accomplished by using a truncated basis set representation of the subsystem A wavefunction. We show that naive truncation of the basis set associated with subsystem A can lead to large numerical artifacts, and present an approach for systematically controlling these artifacts.

Chapter 2 describes the application of the projection-based embedding method to investigate the oxidative stability of lithium-ion batteries. We study the oxidation potentials of mixtures of ethylene carbonate (EC) and dimethyl carbonate (DMC) by using the projection-based embedding method to calculate the vertical ionization energy (IE) of individual molecules at the CCSD(T) level of theory, while explicitly accounting for the solvent using DFT. Interestingly, we reveal that large contributions to the solvation properties of DMC originate from quadrupolar interactions, resulting in a much larger solvent reorganization energy than that predicted using simple dielectric continuum models. Demonstration that the solvation properties of EC and DMC are governed by fundamentally different intermolecular interactions provides insight into key aspects of lithium-ion batteries, with relevance to electrolyte decomposition processes, solid-electrolyte interphase formation, and the local solvation environment of lithium cations.

铌酸锂晶体紫外激光诱导畴反转实验研究

报道了在铌酸锂晶体中实现紫外激光诱导畴反转的实验。在一定外加均匀电场下，铌酸锂晶体中通过波长365nm的紫外激光，由于紫外光的照射降低了矫顽电场只在通光区实现畴反转。研究表明，该方案可用于周期性极化铌酸锂的制备，并有望成为制作精细周期性畴结构的有效技术方案。

掺镁近化学计量比LiTaO3晶体的生长及光学性能

在同成分LiTaO3熔体中掺入一定剂量的K2O,采用顶部籽晶提拉法生长掺镁近化学计量比LiTaO3晶体。对晶体分别进行光谱分析,畴结构和抗光损伤阈值的测定。结果表明：与同成分掺镁LiTaO3晶体相比较,其紫外吸收边出现明显蓝移,红外吸收峰变弱。腐蚀晶片的晶相显微镜观察结果表明：掺镁近化学剂量比晶体的畴结构是较为规则的六边形;晶体的抗光致散射能力明显提高。

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.

Electrochromic effect in domain-inversion process LiNbO3 : Ru : Fe crystals

A reversible electrochromic effect accompanying domain-inversion during the electrical poling process in LiNbO3: Ru: Fe crystals at room temperature has been observed. In electrode area, both electrochromism and domain-inversion occur alternately, and electrochromism is also reversible during back-switch poling, which is experimentally verified and whose mechanism is briefly explained using a microstructure ferroelectric model. In addition, because of the enhancing elcctrochromic effect, different from the undoped LiNbO3 crystals, the coercive riled (21.0 kV/mm or so) measured in LiNbO3: Ru: Fe is lower than its breakdown field, thus providing a possible new technique for realizing the domain-inversion by constant electric field rather than a pulsed one.

Absorption characteristic and nonvolatile holographic recording in LiNbO3 : Cr : Cu crystals

The absorption characteristic of lithium niobate crystals doped with chromium and copper (Cr and Cu) is investigated. We find that there are two apparent absorption bands for LiNbO3:Cr:Cu crystal doped with 0.14 wt.% Cr2O3 and 0.011 wt.% CuO; one is around 480 nm, and the other is around 660 nm. With a decrease in the doping composition of Cr and an increase in the doping composition of Cu, no apparent absorption band in the shorter wavelength range exists. The higher the doping level of Cr, the larger the absorbance around 660 nm. Although a 633 nm red light is located in the absorption band around 660 nm, the absorption at 633 nm does not help the photorefractive process; i.e., unlike other doubly doped crystals, for example, LiNbO3:Fe:Mn crystal, a nonvolatile holographic recording can be realized by a 633 nm red light as the recording light and a 390 nm UV light as the sensitizing light. For LiNbO3:Cr:Cu crystals, by changing the recording light from a 633 nm red light to a 514 nm green light, sensitizing with a 390 nm UV light and a 488 nm blue light, respectively, a nonvolatile holographic recording can be realized. Doping the appropriate Cr (for example, N-Cr = 2.795 X 10(25)m(-3) and N-Cr/N-Cu = 1) benefits the improvement of holographic recording properties. (c) 2005 Optical Society of America.