Towards realistic simulations of lava dome growth using the level set method

The level set method has been implemented in a computational volcanology context. New techniques are presented to solve the advection equation and the reinitialisation equation. These techniques are based upon an algorithm developed in the finite difference context, but are modified to take advantage of the robustness of the finite element method. The resulting algorithm is tested on a well documented Rayleigh–Taylor instability benchmark [19], and on an axisymmetric problem where the analytical solution is known. Finally, the algorithm is applied to a basic study of lava dome growth.

Using the level set method to model endogenous lava dome growth

Modeling volcanic phenomena is complicated by free-surfaces often supporting large rheological gradients. Analytical solutions and analogue models provide explanations for fundamental characteristics of lava flows. But more sophisticated models are needed, incorporating improved physics and rheology to capture realistic events. To advance our understanding of the flow dynamics of highly viscous lava in Peléean lava dome formation, axi-symmetrical Finite Element Method (FEM) models of generic endogenous dome growth have been developed. We use a novel technique, the level-set method, which tracks a moving interface, leaving the mesh unaltered. The model equations are formulated in an Eulerian framework. In this paper we test the quality of this technique in our numerical scheme by considering existing analytical and experimental models of lava dome growth which assume a constant Newtonian viscosity. We then compare our model against analytical solutions for real lava domes extruded on Soufrière, St. Vincent, W.I. in 1979 and Mount St. Helens, USA in October 1980 using an effective viscosity. The level-set method is found to be computationally light and robust enough to model the free-surface of a growing lava dome. Also, by modeling the extruded lava with a constant pressure head this naturally results in a drop in extrusion rate with increasing dome height, which can explain lava dome growth observables more appropriately than when using a fixed extrusion rate. From the modeling point of view, the level-set method will ultimately provide an opportunity to capture more of the physics while benefiting from the numerical robustness of regular grids.

Quantum trajectory simulations of the fluorescence intensity from a two-level atom driven by a multichromatic field

The quantum trajectories method is illustrated for the resonance fluorescence of a two-level atom driven by a multichromatic field. We discuss the method for the time evolution of the fluorescence intensity in the presence of bichromatic and trichromatic driving fields. We consider the special case wherein one multichromatic field component is strong and resonant with the atomic transition whereas the other components are much weaker and arbitrarily detuned from the atomic resonance. We find that the phase-dependent modulations of the Rabi oscillations, recently observed experimentally [Q. Wu, D. J. Gauthier, and T. W. Mossberg, Phys. Rev. A 49, R1519 (1994)] for the special case when the weaker component of a bichromatic driving field is detuned from the atomic resonance by the strong-field Rabi frequency, appear also for detunings close to the subharmonics of the Rabi frequency. Furthermore, we show that for the atom initially prepared in one of the dressed states of the strong field component the modulations are not sensitive to the phase. We extend the calculations to the case of a trichromatic driving field and find that apart from the modulations of the amplitude there is a modulation of the frequency of the Rabi oscillations. Moreover, the time evolution of the fluorescence intensity depends on the phase regardless of the initial conditions and a phase-dependent suppression of the Rabi oscillations can be observed when the sideband fields are tuned to the subharmonics of the strong-field Rabi frequency. [S1050-2947(98)03501-X].

Index of refraction of a two-level atom in a squeezed vacuum field

We study the index of refraction of a two-level atom replacing the usually applied coherent driving fields by a squeezed vacuum field. This system can produce a large index of refraction accompanied by vanishing absorption when the carrier frequency of the squeezed vacuum is detuned from the atomic resonance. (C) 1998 Elsevier Science B.V.

Squeezing-induced complete transparency in two-level systems

The modification of the statistical properties of vacuum fluctuations, via quadrature squeezing, can dramatically reduce the absorptive and dispersive properties of two-level atoms. We show that for some range of parameter values the system exhibits zero absorption accompanied by zero dispersion of the probe field. This complete transparency is attributed to the coherent population oscillations induced by the squeezed vacuum.

Stationary lineshape of a two-level atom in a narrow-bandwidth squeezed vacuum

The stationary lineshape of a two-level atom driven by low-intensity narrow-bandwidth squeezed light is shown to exhibit significant differences in behaviour compared to the lineshape for broadband squeezed light. We find that for narrow-bandwidth squeezed light the lineshape is composed of two Lorentzians whose amplitudes depend on the squeezing correlations. Moreover, one of the Lorentzians has a negative weight which leads to narrowing of the line. These features are absent in the broadband case, where the stationary lineshape is the same as for a thermal field. (C) 1998 Elsevier Science B.V.

Analytical solution for the Mollow and Autler-Townes probe absorption spectra of a three-level atom in a squeezed vacuum

The Mellow and Autler-Townes probe absorption spectra of a three-level atom in a cascade configuration with the lower transition coherently driven and also coupled to a narrow bandwidth squeezed-vacuum field are studied. Analytical studies of the modifications caused by the finite squeezed-vacuum bandwidth to the spectra are made for the case when the Rabi frequency of the driving field is much larger than the natural linewidth. The squeezed vacuum center frequency and the driving laser frequency are assumed equal. We show that the spectral features depend on the bandwidth of a squeezed vacuum field and whether the sources of the squeezing field are degenerate (DPA) or nondegenerate (NDPA) parametric amplifiers. In a broadband or narrow bandwidth squeezed vacuum generated by a NDPA, the central component of the Mellow spectrum can be significantly narrower than that in the normal vacuum. When the source of the squeezed vacuum is a DPA, the central feature is insensitive to squeezing. The Rabi sidebands, however, can be significantly narrowed only in the squeezed vacuum produced by the DPA. The two lines of the Autler-Townes absorption spectrum can be narrowed only in a narrow bandwidth squeezed vacuum, whereas they are independent of the phase and are always broadened in a broadband squeezed vacuum.

Relationship between peptide selectivities of human transporters associated with antigen processing and HLA class I molecules

Efficiency of presentation of a peptide epitope by a MHC class I molecule depends on two parameters: its binding to the MHC molecule and its generation by intracellular Ag processing. In contrast to the former parameter, the mechanisms underlying peptide selection in Ag processing are poorly understood. Peptide translocation by the TAP transporter is required for presentation of most epitopes and may modulate peptide supply to MHC class I molecules. To study the role of human TAP for peptide presentation by individual HLA class I molecules, we generated artificial neural networks capable of predicting the affinity of TAP for random sequence 9-mer peptides. Using neural network-based predictions of TAP affinity, we found that peptides eluted from three different HLA class I molecules had higher TAP affinities than control peptides with equal binding affinities for the same HLA class I molecules, suggesting that human TAP may contribute to epitope selection. In simulated TAP binding experiments with 408 HLA class I binding peptides, HLA class I molecules differed significantly with respect to TAP affinities of their ligands, As a result, some class I molecules, especially HLA-B27, may be particularly efficient in presentation of cytosolic peptides with low concentrations, while most class I molecules may predominantly present abundant cytosolic peptides.

Shift of the subharmonic resonances and suppression of fluorescence in a two-level atom driven by a bichromatic field

We study the behavior of a two-level atom that is driven by a bichromatic field consisting of a strong resonant component and a weaker tunable component. In addition to the splitting of the energy levels (the multiphoton AC Stark effect), we find that the weaker component also shifts the subharmonic resonances, an effect we attribute to a dynamic Stark shift. When the weaker component is tuned to a shifted resonance, no fluorescence occurs at either the frequency of the strong component or the three-photon mixing frequency. Results are obtained with numerical techniques and explained in terms of the dressed-atom model of the system. (C) 1998 Optical Society of America [S0740-3224(98)01508-2] OCIS codes: 270.4180, 270.6620, 270.0270.

Multiphoton ac Stark effect in a bichromatically driven two-level atom

We study the interaction of a two-level atom with two lasers of different frequencies and amplitudes: a strong laser of Rabi frequency 2 Ohm(1) on resonance with the atomic transition, and a weaker laser detuned by subharmonics (2 Ohm(1)/n) of the Rabi frequency of the first. We find that under these conditions the second laser couples the dressed states created by the first in an n-photon process, resulting in doubly dressed states and in a ''multiphoton ac Stark'' effect. We calculate the eigenstates of the doubly dressed atom and their energies, and illustrate the role of this multiphoton ac Stark effect in its fluorescence, absorption, and Autler-Townes spectra. [S1050-2947(98)07607-0].

Two-level atom in a squeezed vacuum with finite bandwidth

The resonance fluorescence of a two-level atom driven by a coherent laser field and damped by a finite bandwidth squeezed vacuum is analysed. We extend the Yeoman and Barnett technique to a non-zero detuning of the driving field from the atomic resonance and discuss the role of squeezing bandwidth and the detuning in the level shifts, widths and intensities of the spectral lines. The approach is valid for arbitrary values of the Rabi frequency and detuning but for the squeezing bandwidths larger than the natural linewidth in order to satisfy the Markoff approximation. The narrowing of the spectral lines is interpreted in terms of the quadrature-noise spectrum. We find that, depending on the Rabi frequency, detuning and the squeezing phase, different factors contribute to the line narrowing. For a strong resonant driving field there is no squeezing in the emitted field and the fluorescence spectrum exactly reveals the noise spectrum. In this case the narrowing of the spectral lines arises from the noise reduction in the input squeezed vacuum. For a weak or detuned driving field the fluorescence exhibits a large squeezing and, as a consequence, the spectral lines have narrowed linewidths. Moreover, the fluorescence spectrum can be asymmetric about the central frequency despite the symmetrical distribution of the noise. The asymmetry arises from the absorption of photons by the squeezed vacuum which reduces the spontaneous emission. For an appropriate choice of the detuning some of the spectral lines can vanish despite that there is no population trapping. Again this process can be interpreted as arising from the absorption of photons by the squeezed vacuum. When the absorption is large it may compensate the spontaneous emission resulting in the vanishing of the fluorescence lines.

Split tolerance to a viral antigen expressed in thymic epithelium and keratinocytes

When expressed as a transgene from the keratin 14 (K14) promoter in an MHC class II-deficient mouse, I-Ab expressed in thymic cortical epithelium promotes positive but not negative selection of I-Ab-restricted CD4(+) T cells (Laufer, T. M. et al., Nature 1996. 383:81-85). Transgenic mice expressing the E7 protein of human papilloma virus 16 from the K14 promoter were studied to determine the consequence of expression of a cytoplasmic/nuclear protein from the K14 promoter. K14E7-transgenic mice express E7 in the thymus and skin without evidence for autoimmunity to E7. Repeated immunization of FVB(H-2(q)) or F1(C57BV6JxFVB) mice with E7 elicited similar antibody responses to the defined B cell epitopes of E7 in K14E7-transgenic and non-transgenic animals. In contrast, for each genetic background, a single immunization with E7 elicited demonstrable T cell proliferative responses to the major promiscuous T helper epitope of E7 in the transgenic but not the non-transgenic animals. Further,E7-immunized non-transgenic F1 (FVBxC57BL/6J) animals developed strong E7-specific cytotoxic T lymphocyte (CTL) responses and were protected against challenge with E7(+) tumors, whereas similarly immunized K14E7-transgenic animals had a markedly reduced CTL response to E7 and no E7-specific tumor protection was observed, although the antibody and CTL response to ovalbumin was normal. Expression of E7 protein as a transgene from the K14 promoter in the skin and thymus thus induces E7-specific tolerance in the cytotoxic T effector repertoire, together with expansion of the E7-specific T helper repertoire. These findings demonstrate that limited tissue distribution of an autoantigen may result in split tolerance to that autoantigen.

Potential strategies utilised by papillomavirus to evade host immunity

The co-evolution of papillomaviruses (PV) and their mammalian hosts has produced mechanisms by which PV might avoid specific and non-specific host immune responses. Low level expression of PV proteins in infected basal epithelial cells, together with an absence of inflammation and of virus-induced cell lysis, restricts the opportunity for effective PV protein presentation to immunocytes by dendritic cells. Additionally, PV early proteins, by a range of mechanisms, may restrict the efficacy of antigen presentation by these cells. Should an immune response be induced to PV antigens, resting keratinocytes (KC) appear resistant to interferon-gamma-enhanced mechanisms of cytotoxic T-lymphocyte (CTL)-mediated lysis, and expression of PV antigens by resting KC can tolerise PV-specific CTL. Thus, KC, in the absence of inflammation, may represent an immunologically privileged site for PV infection. Together, these mechanisms play a parr in allowing persistence of PV-induced proliferative skin lesions for months to years, even in immunocompetent hosts.