Structures and properties of the potassium-doped carbon clusters KCn/KCn +/KCn −(n = 1−10)

We have reported a theoretical study on the interaction mechanism between dust particles in the presence of asymmetric ion flow and an external magnetic field in complex plasma. The recent experimental and numerical results on the particle-wake interaction ensures the dominance of the wake effect in the subsonic regime of plasma flow using the cold ion approximation. The recent developments in dusty plasma research and its growing interest towards more realistic magnetized dusty plasma scenarios also demand serious attention to study the wake effect both in the sub and supersonic regimes in the presence of a magnetic field. It is a challenging task to develop a correct, quantitative theory of wake potential for different regimes of magnetic field and ion flow velocity. Analytic expressions for the wake potential have been reported in this paper for both subsonic and supersonic regimes in the presence of an external magnetic field along with Debye-Hückel type potentials. The results show that the wake potential plays a dominant role in the subsonic regime and its strength increases with an increase in magnetic field. The behaviour of the wake potential is found to have an interesting effect on the Coulomb crystallization of dust grains and is studied with the help of molecular dynamic (MD) simulation.     

An efficient fluid–solid coupling algorithm for single-phase flows

We present a simple and efficient fluid–solid coupling method in two and three spatial dimensions. In particular, we consider the numerical approximation of the Navier–Stokes equations on irregular domains and propose a novel approach for solving the Hodge projection step on arbitrary shaped domains. This method is straightforward to implement and leads to a symmetric positive definite linear system for both the projection step and for the implicit treatment of the viscosity. We demonstrate the accuracy of our method in the L¹$L^1$ and L^∞$L^{\infty }$ norms and present its removing the errors associated with the conventional rasterization-type discretizations. We apply this method to the simulation of a flow past a cylinder in two spatial dimensions and show that our method can reproduce the known stable and unstable regimes as well as correct lift and drag forces. We also apply this method to the simulation of a flow past a sphere in three spatial dimensions at low and moderate Reynolds number to reproduce the known steady axisymmetric and non-axisymmetric flow regimes. We further apply this algorithm to the coupling of flows with moving rigid bodies.     

Gyrotron generation of broadband chaotic radiation under overlapping of high- and low-frequency resonances

Significant (up to 10%) broadening of the band of noise-like radiation in gyrotrons is possible due to specific tuning of gyrofrequency relative to the critical frequency of the working mode when the high- and low-frequency cyclotron resonances are substantially separated at the given translation and transverse velocities of electrons. The generation bands at the resonances are overlapped under operation above the threshold. The analysis with allowance for finiteness of the electron transit time in the interaction space and, hence, different slopes of the dispersion characteristics of electron beam and wave is critically important for correct investigation of the generation regimes in the framework of the average approach. The results are verified using direct 3D particle-in-cell simulation of the chaotic generation regimes of the gyrotron generation in the 8-mm-wavelength range.     

Exact Renormalization Group Analysis of Turbulent Transport by the Shear Flow

The exact renormalization group (RG) method initiated by Wilson and further developed by Polchinski is used to study the shear flow model proposed by Avellaneda and Majda as a simplified model for the diffusive transport of a passive scalar by a turbulent velocity field. It is shown that this exact RG method is capable of recovering all the scaling regimes as the spectral parameters of velocity statistics vary, found by Avellaneda and Majda in their rigorous study of this model. This gives further confidence that the RG method, if implemented in the right way instead of using drastic truncations as in the Yakhot-Orszag’s approximate RG scheme, does give the correct prediction for the large scale behaviors of solutions of stochastic partial differential equations (PDE). We also derive the analog of the “large eddy simulation” models when a finite amount of small scales are eliminated from the problem.     

Broad angle phase matching in subdiffractive photonic crystals

We show that photonic crystals made of materials with normal dispersion allow broad angular range phase matching in nonlinear wave mixing processes if tuned to the subdiffractive (or equivalently self-collimated) beam propagation regimes for the frequencies of both interacting waves. This allows efficient parametric mixing of narrow beams. We demonstrate this idea by numerical simulation of the second harmonic generation in two-dimensional photonic crystal in particular nonlinear material (AlGaAs) in planar waveguide geometry.     

Hydrodynamic binary coalescence of droplets under air flow in a hydrophobic microchannel

Based on the volume of fluid(VOF) method, we conduct a numerical simulation to study the hydrodynamic binary coalescence of droplets under air flow in a hydrophobic rectangular microchannel. Two distinct regimes, coalescence followed by sliding motion and that followed by detaching motion, are identified and discussed. Additionally, the detailed hydrodynamic information behind the binary coalescence is provided, based on which a dynamic mechanical analysis is conducted to reveal the hydrodynamic mechanisms underlying these two regimes. The simulation results indicate that the sliding motion of droplets is driven by the drag force and restrained by the adhesion force induced by the interfacial tension along the main flow direction. The detachment(i.e., upward motion) of the droplet is driven by the lift force associated with an aerodynamic lifting pressure difference imposed on the coalescent droplet, and also restrained by the adhesion force perpendicular to the main flow direction. Especially, the lift force is mainly induced by an aerodynamic lifting pressure difference imposed on the coalescent droplet. Two typical regimes can be quantitatively recognized by a regime diagram depending on Re and We. The higher Re and We respectively lead to relatively larger lift forces and smaller adhesion forces acting on the droplet, both of which are helpful to detachment of the coalesced droplet.     

Simulation of femtosecond pulse propagation in sub-micron diameter tapered fibers

Ultrashort pulse propagation and supercontinuum generation in tapered and microstructured optical fibers is usually simulated using the corrected nonlinear Schrödinger equation. One of the underlying approximations is the use of a wavelength-independent effective area or, equivalently, of a constant nonlinear coefficient . In very thin waveguide structures with strong light confinement, including silica wires and sub-micron tapered fibers and some microstructured fibers, the validity of such an approximation comes into question. In this paper we present an improved model in which all modal properties are fully taken into account as functions of the wavelength. We use comparative numerical simulation to identify certain regimes in which an improved model is needed for quantitatively correct results. PACS 02.70.Hm; 02.60.Cb; 42.65.Wi; 42.65.Re; 42.81.Pp     

p-Wave cold collisions in an optical lattice clock

We study ultracold collisions in fermionic ytterbium by precisely measuring the energy shifts they impart on the atoms' internal clock states. Exploiting Fermi statistics, we uncover p-wave collisions, in both weakly and strongly interacting regimes. With the higher density afforded by two-dimensional lattice confinement, we demonstrate that strong interactions can lead to a novel suppression of this collision shift. In addition to reducing the systematic errors of lattice clocks, this work has application to quantum information and quantum simulation with alkaline-earth atoms.     

Nonlinear dynamics of modulation instability in distributed resonators under external harmonic driving

We study the regimes of complex field dynamics upon modulation instability in distributed nonlinear resonators under external harmonic driving. Two regimes are considered: the regime of a nonlinear ring cavity, described by nonlinear Schrödinger equation (NLS) with a delayed boundary condition, and the regime of a one-dimensional Fabri-Perot cavity, described by a system of coupled NLS for the forward and backward waves. Theoretical stability analysis of stationary forced oscillations is carried out. The results of numerical simulation of transition to chaos with increasing input intensity are presented.     

Self-Modulation and Chaotic Regimes of Generation in a Relativistic Backward-Wave Oscillator with End Reflections

We study self-modulation and chaotic regimes of generation in a relativistic backward-wave oscillator in the presence of reflections of radiation from the boundaries of a slow-wave structure. The onset of self-modulation is examined in detail in the cases of weak and strong reflections. A numerical simulation of transition-to-chaos scenarios is performed over a wide range of parameters. The relation of bifurcation transitions between different regimes to the formation of spatio-temporal structures in the electron beam is discussed.     

Quantitative modeling of laser speckle imaging

We have analyzed the image formation and dynamic properties in laser speckle imaging (LSI) both experimentally and with Monte Carlo simulation. We show for the case of a liquid inclusion that the spatial resolution and the signal itself are both significantly affected by scattering from the turbid environment. Multiple scattering leads to blurring of the dynamic inhomogeneity as detected by LSI. The presence of a nonfluctuating component of scattered light results in the significant increase in the measured image contrast and complicates the estimation of the relaxation time. We present a refined processing scheme that allows a correct estimation of the relaxation time from LSI data.     

Experimental characterization of superradiance in a single-pass high-gain laser-seeded free-electron laser amplifier

In this Letter we report the first experimental characterization of superradiance in a single-pass high-gain free-electron laser (FEL) seeded by a 150 femtosecond (FWHM) Ti:sapphire laser. The nonlinear energy gain after an exponential gain regime was observed. We also measured the evolution of the longitudinal phase space in both the exponential and superradiant regimes. The output FEL pulse duration was measured to be as short as 81 fs, a roughly 50% reduction compared to the input seed laser. The temporal distribution of the FEL radiation as predicted by a numerical simulation was experimentally verified for the first time.     

Saturable absorber mode-locking and laser pulse compression

The theoretical predictions of two alternative models of saturable absorber mode-locking in dye lasers are compared. The predicted mode-locking and pulse compression regimes are considered for both flashlamp-pumped and cw dye lasers, and some questions raised thereby are tested by numerical simulation of pulse propagation.     

Numerical simulation of pool boiling of a Lennard-Jones liquid

We performed a numerical simulation of pool boiling by a molecular dynamics model. In the simulation, a liquid composed of Lennard-Jones particles in a uniform gravitational field is heated by a heat source at the bottom of the system. The model successfully reproduces the change in regimes of boiling from nucleate boiling to film boiling with the increase of the heat source temperature. We present the pool boiling curve by the model, whose general behavior is consistent with those observed in experiments of pool boiling.     

Time scaling regimes in aggregation of magnetic dipolar particles: scattering dichroism results

We report experimental results on the aggregation kinetics in magnetorheological fluids subject to a constant uniaxial magnetic field using the technique of scattering dichroism. We show that the number of aggregated particles displays a long-time power-law dependence with exponents that correspond to two different aggregation regimes. These regimes coincide with 3D and 1D-like aggregation. We also derive the values of both time exponents for the number of aggregated particles.     

CPMG relaxation by diffusion with constant magnetic field gradient in a restricted geometry: numerical simulation and application

Carr-Purcell-Meiboom-Gill (CPMG) measurements are the primary nuclear magnetic resonance (NMR) technique used for evaluating formation properties and reservoir fluid properties in the well logging industry and laboratory sample analysis. The estimation of bulk volume irreducible (BVI), permeability, and fluid type relies on the accurate interpretation of the spin-spin relaxation time (T(2)) distribution. The interpretation is complicated when spin's self-diffusion in an inhomogeneous field and restricted geometry becomes dominant. The combined effects of field gradient, diffusion, and a restricted geometry are not easily evaluated analytically. We used a numerical method to evaluate the dependence of the free and restricted diffusion on the system parameters in the absence of surface relaxation, which usually can be neglected for the non-wetting fluids (e.g., oil or gas). The parameter space that defines the relaxation process is reduced to two dimensionless groups: D* and tau*. Three relaxation regimes: free diffusion, localization, and motionally averaging regimes are identified in the (log(10)D*, log(10)tau*) domain. The hypothesis that the normalized magnetization, M*, relaxes as a single exponential with a constant dimensionless relaxation time T*(2) is justified for most regions of the parameter space. The numerical simulation results are compared with the analytical solutions from the contour plots of T*(2). The locations of the boundaries between different relaxation regimes, derived from equalizing length scales, are challenged by observed discrepancies between numerical and analytical solutions. After adjustment of boundaries by equalizing T*(2), numerical simulation result and analytical solution match each other for every relaxation regime. The parameters, fluid diffusivity and pore length, can be estimated from analytical solutions in the free diffusion and motionally averaging regimes, respectively. Estimation of the parameters near the boundaries of the regimes may require numerical simulation.     

Quantum simulation of quantum field theories in trapped ions

We propose the quantum simulation of fermion and antifermion field modes interacting via a bosonic field mode, and present a possible implementation with two trapped ions. This quantum platform allows for the scalable add up of bosonic and fermionic modes, and represents an avenue towards quantum simulations of quantum field theories in perturbative and nonperturbative regimes.