High-Order Harmonics from Bow Wave Caustics Driven by a High-Intensity Laser

We propose a new mechanism of high-order harmonic generation during an interaction of a high-intensity laser pulse with underdense plasma. A tightly focused laser pulse creates a cavity in plasma pushing electrons aside and exciting the wake wave and the bow wave. At the joint of the cavity wall and the bow wave boundary, an annular spike of electron density is formed. This spike surrounds the cavity and moves together with the laser pulse. Collective motion of electrons in the spike driven by the laser field generates high-order harmonics. A strong localization of the electron spike, its robustness to oscillations imposed by the laser field and, consequently, its ability to produce high-order harmonics is explained by catastrophe theory. The proposed mechanism explains the experimental observations of high-order harmonics with the 9 TW J-KAREN laser (JAEA, Japan) and the 120 TW Astra Gemini laser (CLF RAL, UK) [A. S. Pirozhkov, et al., arXiv:1004.4514 (2010); A. S. Pirozhkov et al, AIP Proceedings, this volume]. The theory is corroborated by high-resolution two- and three-dimensional particle-in-cell simulations.

Unsteady Simulations of a Film Cooling Flow from an Inclined Cylindrical Jet

Film cooling is extensively used to provide protection against the severe thermal environment in gas turbine engines. Most of the computational studies on film cooling flow have been done using steady Reynolds-averaged Navier–Stokes calculation procedures. However, the flowfield associated with a jet in a crossflow is highly unsteady and complex with different types of vortical structures. In this paper, a computational investigation about the unsteady phenomena of a jet in a crossflow is performed using detached eddy simulation. Detailed computation of a single row of 35 deg round holes on a flat plate has been obtained for a 1.0 blowing ratio and a 2.0 density ratio. First, time-step size, grid resolution, and computational domain tests for an unsteady simulation have been conducted. Comparison between the results of unsteady Reynolds-averaged Navier–Stokes calculation, detached eddy simulation, and large eddy simulation is also performed. Comparison of the time-averaged detached eddy simulation prediction with the measured film-cooling effectiveness shows that the detached eddy simulation prediction is reasonable. From present detached eddy simulations, the influential coherent vortical structures of a film cooling flow can be seen. The unsteady physics of jet in a crossflow interactions and a jet liftoff in film cooling flows have been explained.

Improved description of chemical barriers with generalized gradient approximations (GGAs) and meta-GGAs

The performance of exchange and correlation (xc) functionals of the generalized gradient approximation (GGA) type and of the meta-GGA type in the calculation of chemical reactions is related to topological features of the electron density which, in turn, are connected to the orbital structure of chemical bonds within the Kohn-Sham (KS) theory. Seventeen GGA and meta-GGA xc functionals are assessed for 15 hydrogen abstraction reactions and 3 symmetrical S(N)2 reactions. Systems that are problematic for standard GGAs characteristically have enhanced values of the dimensionless gradient argument s(sigma)(2) with local maxima in the bonding region. The origin of this topological feature is the occupation of valence KS orbitals with an antibonding or essentially nonbonding character. The local enhancement of s(sigma)(2) yields too negative exchange-correlation energies with standard GGAs for the transition state of the S(N)2 reaction, which leads to the reduced calculated reaction barriers. The unwarranted localization of the effective xc hole of the standard GGAs, i.e., the nondynamical correlation that is built into them but is spurious in this case, wields its effect by their s(sigma)(2) dependence. Barriers are improved for xc functionals with the exchange functional OPTX as x component, which has a modified dependence on s(sigma)(2). Standard GGAs also underestimate the barriers for the hydrogen abstraction reactions. In this case the barriers are improved by correlation functionals, such as the Laplacian-dependent (LAP3) functional, which has a modified dependence on the Coulomb correlation of the opposite- and like-spin electrons. The best overall performance is established for the combination OLAP3 of OPTX and LAP3.

Modelling the H Lyman lines in evolved late-type stars

The importance of partial redistribution (PRD) in the modelling of the Lyman a and Lyman ß emission lines of hydrogen in stellar atmospheres is examined using simple atmospheric models of a range of late-type stars. These models represent the subgiant Procyon (F5 IV-V), and the two giants ß Gem (K0 III) and a Tau (K5 III). These stars are selected to span a wide range of surface gravities: 1.25 <log g <4.00. The calculations are performed using the computer code MULTI with the modifications made by Hubeny & Lites. It is found that PRD effects are highly significant, both in the direct prediction of the Lyman line profiles and in the application of hydrostatic equilibrium to calculate the atmospheric electron density in static atmospheric models.

PIC simulation of a thermal anisotropy-driven Weibel instability in a circular rarefaction wave

The expansion of an initially unmagnetized planar rarefaction wave has recently been shown to trigger a thermal anisotropy-driven Weibel instability (TAWI), which can generate magnetic fields from noise levels. It is examined here whether the TAWI can also grow in a curved rarefaction wave. The expansion of an initially unmagnetized circular plasma cloud, which consists of protons and hot electrons, into a vacuum is modelled for this purpose with a two-dimensional particle-in-cell (PIC) simulation. It is shown that the momentum transfer from the electrons to the radially accelerating protons can indeed trigger a TAWI. Radial current channels form and the aperiodic growth of a magnetowave is observed, which has a magnetic field that is oriented orthogonal to the simulation plane. The induced electric field implies that the electron density gradient is no longer parallel to the electric field. Evidence is presented here that this electric field modification triggers a environments, which are needed to explain the electromagnetic emissions by astrophysical jets. It is outlined how this instability could be examined experimentally.second magnetic instability, which results in a rotational low-frequency magnetowave. The relevance of the TAWI is discussed for the growth of small-scale magnetic fields in astrophysical

Investigation of magnesium laser ablated plumes with Thomson scattering

Optical Thomson scattering has been implemented as a diagnostic of laser ablated plumes generated with second harmonic Nd:YAG laser radiation at 532 nm. Thomson scattering data with both spatial and temporal resolution has been collected, giving both electron density, and temperature distributions within the plume as a function of time. Although the spatial profiles do not match very well for simple models assuming either isothermal or isentropic expansion, consideration of the measured ablated mass indicates an isothermal expansion fits better than an isentropic expansion and indeed, at late time, the spatial profile of temperature is almost consistent with an isothermal approximation.

Comparing Optical Thomson Scattering With Emission Spectroscopy In He Plasmas Induced By Laser Breakdown At 1 Atmosphere

We have carried out optical Thomson scattering measurements from a laser induced breakdown in He at 1 atmosphere. The breakdown was created with a Nd:YAG laser with 9ns pulse duration and 400mJ pulse energy focused into a chamber filled with He. A second harmonic Nd: YAG laser with 9ns pulses and up to 80mJ energy was used to obtain temporally and spatially resolved data on the electron density and temperature. In parallel experiments, we measured the emission of the 447.1nm line from He I. Initial results suggest good agreement between densities inferred but full Abel inversion is needed for conclusive results.

Proton ordering energetics in ice phases

Results from first-principles calculations on the subtle energetics of proton ordering in ice phases are shown only to depend on the electrostatic components of the total energy. Proton ordered ice phases can therefore be predicted using electronic structure methods or a tailored potential model. However, analysis of the electron density reveals that high order multipole components, up to hexadecapole, are needed to adequately capture total energy differences between proton ordered and disordered phases. This suggests that current potential models may be unable to reproduce the position of proton ordered ice phases in the phase diagram without extensions to describe high order electrostatics. (c) 2006 Elsevier B.V. All rights reserved.

Dual Enlargement of Gold Nanoparticles: From Mechanism to Scanometric Detection of Pathogenic Bacteria

A mechanism of dual enlargement of gold nanoparticles (AuNPs) comprising two steps is described. In the first step, the AuNPs are enlarged by depositing Au atoms on their crystalline faces. In this process, the particles are not only enlarged but they are also observed to multiply: new Au nuclei are formed by the budding and division of the enlarged particles. In the second step, a silver enhancement is subsequently performed by the deposition of silver atoms on the enlarged and newly formed AuNPs to generate bimetallic Au@Ag core-shell structures. The dual nanocatalysis greatly enhances the electron density of the nanostructures, leading to a stronger intensity for colorimetric discrimination as well as better sensitivity for quantitative measurement. Based on this, a simple scanometric assay for the on-slide detection of the food-born pathogen Campylobacter jejuni is developed. After capturing the target bacteria, gold-tagged immunoprobes are added to create a signal on a solid substrate. The signal is then amplified by the dual enlargement process, resulting in a strong color intensity that can easily be recognized by the unaided eye, or measured by an inexpensive flatbed scanner. In this paper, dual nanocatalysis is reported for the first time. It provides a valuable mechanistic insight into the development of a simple and cost-effective detection format.

Conjugate heat transfer and film cooling of a multi-stage cooling scheme

Cooling techniques play a key role in improving efficiency and power output of modern gas turbines. The conjugate technique of film and impingement cooling schemes is considered in this study. The Multi-Stage Cooling Scheme (MSCS) involves coolant passing from inside to outside turbine blade through two stages. The first stage; the coolant passes through first hole to internal gap where the impinging jet cools the external layer of the blade. Finally, the coolant passes through the internal gap to the second hole which has specific designed geometry for external film cooling. The effect of design parameters, such as, offset distance between two-stage holes, gap height, and inclination angle of the first hole, on upstream conjugate heat transfer rate and downstream film cooling effectiveness performance are investigated computationally. An Inconel 617 alloy with variable properties is selected for the solid material. The conjugate heat transfer and film cooling characteristics of MSCS are analyzed across blowing ratios of Br = 1 and 2 for density ratio, 2. This study presents upstream wall temperature distributions due to conjugate heat transfer for different gap design parameters. The maximum film cooling effectiveness with upstream conjugate heat transfer is less than adiabatic film cooling effectiveness by 24–34%. However, the full coverage of cooling effectiveness in spanwise direction can be obtained using internal cooling with conjugate heat transfer, whereas adiabatic film cooling effectiveness has narrow distribution.

Unsteady heat transfer analysis of a film cooling flow

Unsteady heat transfer in a turbine blade film cooling flow is studied using detached eddy simulation (DES). Detailed computation of a single row of 35 degree round holes on a flat plate has been obtained for a blowing ratio of 1.0 and a density ratio of 2.0. The instantaneous flow fields and heat transfer distributions are found to be highly unsteady and oscillatory in nature. The fluctuation of the adiabatic effectiveness and heat transfer coefficient, for example, can be as high as 15 and 50 percent of the time-averaged value, respectively. The correlation between the coherent vortical structures and the unsteady heat transfer is carefully examined. It is shown that the fluctuations in the adiabatic effectiveness and heat transfer coefficient are mainly caused by the spanwise fluctuation of the coolant jet and the thermal turbulent boundary layer accompanying the unsteady flow structures.

Interacting kinetics of neutral and ionic species in an atmospheric-pressure helium-oxygen plasma with humid air impurities

We unravel the complex chemistry in both the neutral and ionic systems of a radio-frequency-driven atmospheric-pressure plasma in a helium-oxygen mixture (He-0.5% O) with air impurity levels from 0 to 500 ppm of relative humidity from 0% to 100% using a zero-dimensional, time-dependent global model. Effects of humid air impurity on absolute densities and the dominant production and destruction pathways of biologically relevant reactive neutral species are clarified. A few hundred ppm of air impurity crucially changes the plasma from a simple oxygen-dependent plasma to a complex oxygen-nitrogen-hydrogen plasma. The density of reactive oxygen species decreases from 10 to 10 cm, which in turn results in a decrease in the overall chemical reactivity. Reactive nitrogen species (10 cm ), atomic hydrogen and hydroxyl radicals (10-10 cm) are generated in the plasma. With 500 ppm of humid air impurity, the densities of positively charged ions and negatively charged ions slightly increase and the electron density slightly decreases (to the order of 10 cm). The electronegativity increases up to 2.3 compared with 1.5 without air admixture. Atomic hydrogen, hydroxyl radicals and oxygen ions significantly contribute to the production and destruction of reactive oxygen and reactive nitrogen species. © 2013 IOP Publishing Ltd.

Universal tight binding model for chemical reactions in solution and at surfaces. I. Organic molecules

As is now well established, a first order expansion of the Hohenberg-Kohn total energy density functional about a trial input density, namely, the Harris-Foulkes functional, can be used to rationalize a non self consistent tight binding model. If the expansion is taken to second order then the energy and electron density matrix need to be calculated self consistently and from this functional one can derive a charge self consistent tight binding theory. In this paper we have used this to describe a polarizable ion tight binding model which has the benefit of treating charge transfer in point multipoles. This admits a ready description of ionic polarizability and crystal field splitting. It is necessary in constructing such a model to find a number of parameters that mimic their more exact counterparts in the density functional theory. We describe in detail how this is done using a combination of intuition, exact analytical fitting, and a genetic optimization algorithm. Having obtained model parameters we show that this constitutes a transferable scheme that can be applied rather universally to small and medium sized organic molecules. We have shown that the model gives a good account of static structural and dynamic vibrational properties of a library of molecules, and finally we demonstrate the model's capability by showing a real time simulation of an enolization reaction in aqueous solution. In two subsequent papers, we show that the model is a great deal more general in that it will describe solvents and solid substrates and that therefore we have created a self consistent quantum mechanical scheme that may be applied to simulations in heterogeneous catalysis.

Investigating the dynamics of laser induced sparks in atmospheric helium using Rayleigh and Thomson scattering

We have used optical Rayleigh and Thomson scattering to investigate the expansion dynamics of laser induced plasma in atmospheric helium and to map its electron parameters both in time and space. The plasma is created using 9 ns duration, 140 mJ pulses from a Nd:YAG laser operating at 1064 nm, focused with a 10 cm focal length lens, and probed with 7 ns, 80 mJ, and 532 nm Nd:YAG laser pulses. Between 0.4 μs and 22.5 μs after breakdown, the electron density decreases from 3.3 × 1017 cm−3 to 9 × 1013 cm−3, while the temperature drops from 3.2 eV to 0.1 eV. Spatially resolved Thomson scattering data recorded up to 17.5 μs reveal that during this time the laser induced plasma expands at a rate given by R ∼ t0.4 consistent with a non-radiative spherical blast wave. This data also indicate the development of a toroidal structure in the lateral profile of both electron temperature and density. Rayleigh scattering data show that the gas density decreases in the center of the expanding plasma with a central scattering peak reemerging after about 12 μs. We have utilized a zero dimensional kinetic global model to identify the dominant particle species versus delay time and this indicates that metastable helium and the He2 + molecular ion play an important role.

Numerical modelling of low temperature laser produced magnesium plasma

The two-dimensional laser-plasma-interaction hydrodynamic code POLLUX has been used to simulate the ablation of a magnesium target by a 30-ns, 248-nm KrF excimer laser at low laser fluences of ≤10 J cm2. This code, originally written for much higher laser intensities, has been recently extended to include a detailed description of the equation of state in order to treat changes of phase within the target material, and also includes a Thomas Fermi description of the electrons. The simulated temporal and spatial evolution of the plasma plume in the early phase of the expansion (≤100 ns) is compared with experimental interferometric measurements of electron density. The expansion dynamics are in good agreement, although the simulated electron number density is about 2.5 times higher than the experimental values.