Exact results and open questions in first principle functional RG

Some aspects of the functional RG (FRG) approach to pinned elastic manifolds (of internal dimension d) at finite temperature T > 0 are reviewed and reexamined in this much expanded version of Le Doussal (2006) [67]. The particle limit d = 0 provides a test for the theory: there the FRG is equivalent to the decaying Burgers equation, with viscosity ν ∼ T-both being formally irrelevant. An outstanding question in FRG, i.e. how temperature regularizes the otherwise singular flow of T = 0 FRG, maps to the viscous layer regularization of inertial range Burgers turbulence (i.e. to the construction of the inviscid limit). Analogy between Kolmogorov scaling and FRG cumulant scaling is discussed. First, multi-loop FRG corrections are examined and the direct loop expansion at T > 0 is shown to fail already in d = 0, a hierarchy of ERG equations being then required (introduced in Balents and Le Doussal (2005) [36]). Next we prove that the FRG function R(u) and higher cumulants defined from the field theory can be obtained for any d from moments of a renormalized potential defined in an sliding harmonic well. This allows to measure the fixed point function R(u) in numerics and experiments. In d = 0 the beta function (of the inviscid limit) is obtained from first principles to four loop. For Sinai model (uncorrelated Burgers initial velocities) the ERG hierarchy can be solved and the exact function R(u) is obtained. Connections to exact solutions for the statistics of shocks in Burgers and to ballistic aggregation are detailed. A relation is established between the size distribution of shocks and the one for droplets. A droplet solution to the ERG functional hierarchy is found for any d, and the form of R(u) in the thermal boundary layer is related to droplet probabilities. These being known for the d = 0 Sinai model the function R(u) is obtained there at any T. Consistency of the ?=4-d expansion in one and two loop FRG is studied from first principles, and connected to shock and droplet relations which could be tested in numerics.     

New high negative dispersion photonic crystal fiber

A new high negative dispersion photonic crystal fiber is proposed. It has double-core structure. The inner core has a circle germanium-doped region. The outer core is formed by removing the 3rd ring air-holes around the core. There are two ring air-holes between the two cores, Diameter of the 1st ring air holes is bigger than that of the 2nd ring air-holes, this can make mode coupling between inner mode and outer mode and showed that the high negative PCF is the result of this structure characteristics. There are honeycomb photonic lattice in the PCF's cladding. The influence of the structure parameters deviated from the design those on the chromatic dispersion are evaluated. When the structure parameters Λ=1.50 μm, d_core=2.10 μm, d₁=0.90 μm, d₂=0.44 μm and d₃=1.04 μm, the dispersion coefficient D is −1320 ps/(nm·km) at 1550 nm. This is a new kind of chromatic dispersion compensation PCF.     

Thermal fluctuations in pinned elastic systems: field theory of rare events and droplets

Using the functional renormalization group (FRG) we study the thermal fluctuations of elastic objects (displacement field u, internal dimension d) pinned by a random potential at low temperature T, as prototypes for glasses. A challenge is how the field theory can describe both typical (minimum energy T = 0) configurations, as well as thermal averages which, at any non-zero T as in the phenomenological droplet picture, are dominated by rare degeneracies between low lying minima. We show that this occurs through an essentially non-perturbative thermal boundary layer (TBL) in the (running) effective action Γ [u] at T > 0 for which we find a consistent scaling ansatz to all orders. The TBL describes how temperature smoothes the singularities of the T = 0 theory and contains the physics of rare thermal excitations (droplets). The formal structure of this TBL, which involves all cumulants of the coarse grained disorder, is first explored around d = 4 using a one-loop Wilson RG. Next, a more systematic exact RG (ERG) method is employed, and first tested on d = 0 models where it can be pushed quite far. There we obtain precise relations between TBL quantities and droplet probabilities (those are constrained by exact identities which are then checked against recent exact results). Our analysis is then extended to higher d, where we illustrate how the TBL scaling remains consistent to all orders in the ERG and how droplet picture results can be retrieved. Since correlations are determined deep in the TBL (by derivatives of Γ [u] at u = 0), it remains to be understood (in any d) how they can be retrieved (as u = 0⁺ limits in the non-analytic T = 0 effective action), i.e., how to recover a T = 0 critical theory. This formidable “matching problem” is solved in detail for d = 0, N = 1 by studying the (partial) TBL structure of higher cumulants when points are brought together. We thereby obtain the β-function at T = 0, all ambiguities removed, displayed here up to four loops. A discussion of the d > 4 case and an exact solution at large d are also provided.     

Adsorption in cylindrical pores: Mixed lattice-site/off-site Monte Carlo simulations in pores with heterogeneous wall structure

We present results of grand canonical Monte Carlo simulations of adsorption in cylindrical pores with rough surface modeled by lattice-site approach. Each site is characterized by two parameters: structural and energetic, which locally modify the structure and energy properties of the surface. There are three types of sites, randomly distributed over the wall: attractive, neutral and repulsive with respect to the smooth pore model. The results presented here show how this model affects the mechanism of adsorption and how it changes the forms of adsorption isotherm. We compare our numerical results with the experimental data of adsorption of a simple fluid (CH₄, T = 77 K) in cylindrical silica pore of diameter d = 4 nm (MCM-41 material).     

A two dimensional theoretical model for describing gain coefficient in N2-lasers and the model validity for CVL and excimer lasers

Based on our previously reported measurements on the gain-value in a N₂- laser and numerical calculations, we introduce a method to obtain an analytical expression for the small signal gain, g₀, where the dependency of g₀ on the laser geometrical configuration, including electrodes length and gap separation, is demonstrated. For this study one- and two-dimensional approaches for the photon density have been applied independently to determine gain-parameter, where for explaining the observed dependency of the gain-parameter on the laser electrodes separation, d_AMP, which was found experimentally and explained by an empirical expression of the type g₀ = r + q/d_AMP, with r and q some constants, realization of introducing an extra dimension along the gap separation was unavoidable. For the electrodes length, l_AMP, we have already shown that an empirical equation of the type g₀ = m + n/l_AMP, with m and n some constants, is consistent with the measurements corresponding to N₂-lasers. With this realization, it is proved that the gain-parameter in N₂-lasers can be written as g₀^{above threshold} = m″ + g_0z(γ_L^z, z) + g_0y(γ_L^y, y), where it consists of two independent gain-values along the electrodes length (z) and gap separation (y) with the corresponding power losses given by γ_L^z and γ_L^y. m″ is a very small quantity showing that laser is operating slightly above the threshold. The results of this calculation are consistent with our recent measurements and also other reported N₂-laser gain values measured under moderate current density conditions. To check the validity of the model for other types of lasers, the reported gain-values for copper vapor lasers of different laser tube radii, R_AMP, and tube lengths, l_AMP, have been examined using the one-dimensional model of either g₀(R_AMP) or g₀(l_AMP) and the consistency with the observed measurements was found to be quite satisfactory. The model was also found to be valid for the excimer lasers of different types, different gas mixtures and pressures at a constant input operational voltage. Due to limited numbers of the reported experimental measurements, for the graphs preparation of g₀(l_AMP) in excimer lasers we used the observed data at V₀ = 30 kV and also some variations of the input voltages in the range of ΔV ≅ 20 kV have been adopted. The results for both cases were found to be consistent with the proposed one-dimensional model.     

Effect of the surface integral in the torque equation of the electromagnetic angular momentum on the Faraday rotation

A circularly polarized plane wave of infinite transverse extent (δ = ∞) has no spin angular momentum, while a realistic light does carry it. This paradox originates from the presence (δ = ∞) and absence (δ ≈ 0) of the surface integral in the total angular momentum J. The same holds for the torque equation of dJ/dt, so that δ is also connected with the relative Faraday rotation angle Θ_F/θ_F when a radius (a) of a cylindrical medium with optical activity is only a little larger than that (b) of light beam, where Θ_F is the Faraday rotation angle and θ_F is the intrinsic Faraday rotation angle of a medium. It is shown here that it is possible to estimate δ for a realistic light from the drastic variation in Θ_F/θ_F near b/a = 1.     

Bose-Einstein solitons in time-dependent linear potential

In this paper, Bose-Einstein soliton solutions of the nonlinear Schrödinger equation with time-dependent linear potential are considered. Based on the F-expansion method, we present a number of Jacobian elliptic function solutions. Particular cases of these solutions, where the elliptic function modulus equals 1 and 0, are various localized solutions and trigonometric functions, respectively. Specially, for V_ext = ZF(T) = Z[mg + Hcos (ω₁T)], we discussed the Bose-Einstein condensate trapped in the coupling external field with considering the effect of gravity; for F(T) = constant, it describes the wave (Langmuir or electromagnetic) in a linearly inhomogeneous plasma with cubic nonlinearly.     

High-order incompressible large-eddy simulation of fully inhomogeneous turbulent flows

The subgrid-scale (SGS) eddy-viscosity model developed by Vreman [Phys. Fluids 16 (2004) 3670] and its dynamic version [Phys. Fluids 19 (2007) 065110] are tested in large-eddy simulations (LES) of the turbulent flow in an Re = 12,000 lid-driven cubical cavity by comparison to the direct numerical simulation (DNS) data of Leriche and Gavrilakis [Phys. Fluids 12 (2000) 1363]. This appears to be the first test of this class of model to flows without any homogeneous flow directions, which is typical of flows in complex geometries. Additional LES predictions at Re = 18,000 and Re = 22,000 are compared to the DNS data of Leriche [J. Sci. Comp. 27 (2006)]. The new LES framework yielded excellent agreement for both the mean velocity and Reynolds stress profiles and matches DNS data better than the more traditional Smagorinsky-based SGS models.     

Evidence for diffusion-limited kinetics during electromigration-induced step bunching on Si(1 1 1)

We have measured how the initial terrace width l₀ on vicinal Si(1 1 1) surfaces influences the rate of step bunching and the minimum terrace width within a bunch when direct-current heated at 940-1290 °C. A comparison of this data with analytic solutions and numerical simulations of the conventional “sharp-step” model give strong evidence that the kinetic length d is relatively small (d < ∼20 nm) in both temperature regime I (∼850-950 °C) and regime III (∼1200-1300 °C), in which step-down current is required for step bunching. This indicates that surface mass transport is diffusion-limited in both regimes I and III when l₀ > 20 nm, and hence that the adatom attachment- and terrace diffusion-hopping rates are of comparable magnitude. We also observe similar scaling with initial terrace width in temperature regime II (∼1040-1190 °C), in which step-up current is required for bunching, suggesting a similar step bunching mechanism in all three temperature regimes.     

One-dimensional models of internal rotation in CX3NO molecules (X = H, D, F)

Various non-empirical methods for estimating the parameters of one-dimensional internal rotation potentials and energies of torsional transitions were compared for the CX₃NO molecules (X = H, D, F) in the ground (S₀) and lowest excited singlet (S₁) electronic states. The potential energy surfaces were studied by the ab initio MR-AQCC/cc-pVTZ, MR-AQCC/cc-pVTZ(-f), MP2/6-311G(2d), and MP2/6-311G(d,p) methods. The one-dimensional internal rotation problem was solved using the following models: (1) geometry optimization at a given internal rotation coordinate; (2) intrinsic reaction path; (3) gradient extremal; and (4) use of only the data on potential energy surface stationary points. Special attention was paid to the problem of calculation of kinematic coefficient. In all cases, the calculated torsional energies for CX₃NO molecules (X = H, D, F) are in agreement with experiment. The results from different methods for constructing torsional cross-sections of the potential energy surface are virtually equivalent and differ insignificantly from the results of calculations within the framework of the simplest model, hence, estimates of the barrier to internal rotation are of most importance. It was found that a change in the zero-point energy could give a correction to the internal rotation potential as large as 15% of the potential barrier. However, in the case under consideration the calculations in the harmonic approximation taking into account this correction do not improve the agreement between the calculated torsional transitions and the experimental data.     

The influence of a demagnetizing field on hysteresis losses in a dense assembly of superparamagnetic nanoparticles

The electrodynamic method is used to measure the hysteresis losses of a dense assembly of magnetite nanoparticles with an average diameter D=25 nm in the frequency range f=10–150 kHz and for magnetic field amplitudes H₀=100–300 Oe. It is found that the specific loss power is determined by a demagnetizing factor of a whole sample. It diminishes approximately 4.5 times when the sample aspect ratio decreases from L/d=11.4 to L/d≈1, where L and d are the sample length and diameter, respectively. For H₀≤300 Oe the maximal specific loss power 120 W/g is obtained for the sample with L/d=11.4 at f=120 kHz. For comparison, the assembly specific absorption rate has been determined also by means of direct measurement of the temperature difference between the inner and outer surfaces of a flat cuvette containing magnetic nanoparticles. For both methods of measurement close values for the specific absorption rate are obtained for samples with similar demagnetizing factors.     

Quasi-particle excitations around a pair of half-quantum vortices in p- wave superconductors

Triplet superconductors such as Sr₂RuO₄ and Na_xCoO₂⋅yH₂O are now found to be p-wave (p_x ± ip_y) or f-wave ((p_x ± ip_y) cos cp_z) superconductors. It was phenomenologically suggested that in these p-wave or f-wave superconductors, a pair of half-quantum vortices (HQVs) becomes stable. Using the Bogoliubov-de Gennes equation, previously we have analyzed quasi-particle excitations around an HQV at one end of a d-soliton for simplicity. In next study, we will investigate the stability of the pair of HQV’s, which are connected by the d-soliton. For this purpose, we have developed a new numerical method to solve the Bogoliubov-de Gennes equation for two vortices state, using Mathieu functions.     

A lattice Boltzmann model with an amending function forsimulating nonlinear partial differential equations

This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations (NPDEs) in the form $u_t+\alpha uu_{xx}+\beta u^n u_x+\gamma u_{xxx}+\xi u_{xxxx}=0$. This model is different from existing models because it lets the time step be equivalent to the square of the space step and derives higher accuracy and nonlinear terms in NPDEs. With the Chapman--Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The numerical results agree well with the analytical solutions.     

Optical sensing characteristics of the photonic crystal fiber filled with magnetic fluid

A novel photonic crystal fiber sensing theory filled with magnetic fluid is proposed based on the change of the MF refractive index under varied magnetic field. The magnetically induced tuning of the magnetic fluid filled PCF propagation properties were investigated by the full-vector finite element method with a perfectly matched layer. Theoretical calculations show that both the effective refractive index and the effective mode area increase vs. the increased magnetic field, and the PCF filled MF with larger d/Λ is more sensitive to magnetic field. When the wavelength λ = 1550 nm, the duty ratio d/Λ = 0.9, d/Λ = 0.6, the effective refractive indexes increase respectively from 1.598279 to 1.617572, from 1.61948 to 1.632484, and the effective mode areas increase respectively from 3.561115 μm² to 7.052360 μm², from 6.167494 μm² to 37.221998 μm² as the magnetic field changes from 25 Oe to 175 Oe. This scheme provides theoretical foundation to use magnetic field to control light in photonic crystal fiber and also offers a potential method for magnetic field sensing based on the TIR-PCF.     

An efficient numerical method for computing dynamics of spin F = 2 Bose–Einstein condensates

In this paper, we extend the efficient time-splitting Fourier pseudospectral method to solve the generalized Gross–Pitaevskii (GP) equations, which model the dynamics of spin F = 2 Bose–Einstein condensates at extremely low temperature. Using the time-splitting technique, we split the generalized GP equations into one linear part and two nonlinear parts: the linear part is solved with the Fourier pseudospectral method; one of nonlinear parts is solved analytically while the other one is reformulated into a matrix formulation and solved by diagonalization. We show that the method keeps well the conservation laws related to generalized GP equations in 1D and 2D. We also show that the method is of second-order in time and spectrally accurate in space through a one-dimensional numerical test. We apply the method to investigate the dynamics of spin F = 2 Bose–Einstein condensates confined in a uniform/nonuniform magnetic field.     

Quantum-mechanical vs. semi-classical spectral-line widths and shifts from the line core in the non-impact region for the Ar-perturbed/ K-radiator system

New quantum-mechanical (QM) and semi-classical (SC) shifts (d's) and widths (HWHM's, w's) were measured from the line core of computed full spectral-line shapes for the Ar-perturbed/K-radiator system (K/Ar). The initial state of our model was based on a 4p²P_3/2,1/2 pseudo-potential for the K/Ar system, and the final state on a zero potential. The Fourier transform of the line shape formed the basis for the computations. Excellent agreement was found between the QM and SC values of d and of w in a high-pressure (P) non-impact region, which was characterized by a √P dependence of w and a P dependence of d. These agreements were shown to be another example of a correspondence between classical (SC) quantities and QM quantities in the limit of large quantum numbers. Typically at P=1×10⁶ Torr and T=400 K, w_QM=448 cm⁻¹ and w_SC=479 cm⁻¹, where the deviation from the mean is ±3.3%. Also, d_QM=−3815 cm⁻¹ and d_SC=−3716 cm⁻¹, where the deviation from the mean is ±1.3%. A new general method was formulated which yielded a definite pressure P₀, which was defined as an upper limit to the low-pressure impact approximation and a lower limit to the non-impact region.     

Two-boson truncation of Pauli-Villars-regulated Yukawa theory

We apply light-front quantization, Pauli-Villars regularization, and numerical techniques to the nonperturbative solution of the dressed-fermion problem in Yukawa theory in 3 + 1 dimensions. The solution is developed as a Fock-state expansion truncated to include at most one fermion and two bosons. The basis includes a negative-metric heavy boson and a negative-metric heavy fermion to provide the necessary cancellations of ultraviolet divergences. The integral equations for the Fock-state wave functions are solved by reducing them to effective one-boson-one-fermion equations for eigenstates with J_z = 1/2. The equations are converted to a matrix equation with a specially tuned quadrature scheme, and the lowest mass state is obtained by diagonalization. Various properties of the dressed-fermion state are then computed from the nonperturbative light-front wave functions. This work is a major step in our development of Pauli-Villars regularization for the nonperturbative solution of four-dimensional field theories and represents a significant advance in the numerical accuracy of such solutions.     

On the well-posedness of the stochastic Allen–Cahn equation in two dimensions

White noise-driven nonlinear stochastic partial differential equations (SPDEs) of parabolic type are frequently used to model physical systems in space dimensions d = 1, 2, 3. Whereas existence and uniqueness of weak solutions to these equations are well established in one dimension, the situation is different for d ? 2. Despite their popularity in the applied sciences, higher dimensional versions of these SPDE models are generally assumed to be ill-posed by the mathematics community. We study this discrepancy on the specific example of the two dimensional Allen–Cahn equation driven by additive white noise. Since it is unclear how to define the notion of a weak solution to this equation, we regularize the noise and introduce a family of approximations. Based on heuristic arguments and numerical experiments, we conjecture that these approximations exhibit divergent behavior in the continuum limit. The results strongly suggest that shrinking the mesh size in simulations of the two-dimensional white noise-driven Allen–Cahn equation does not lead to the recovery of a physically meaningful limit.     

Microwave spectrum, molecular structure, and theoretical calculation of cyclohexanone oxime (chair form)

The microwave spectra of cyclohexanone oxime and d₁ (=NOD) and d₄(2,2,6,6-d₄) derivatives were observed in the frequency range from 8 to 40 GHz in the ground and excited vibrational states. The rotational constants were determined to be A = 3799.844(48), B = 1513.7912(23), and C = 1189.6118(29) MHz for normal species, A = 3791.835(88), B = 1461.0324(47), and C = 1157.5653(53) MHz for d₁ species, and A = 3364.141(49), B = 1487.9551(34), and C = 1154.0965(44) MHz for d₄ species in the ground vibrational state. The planar moments, P_bb (P_bb = (I_c + I_a − I_b)/2) of normal, d₁, and d₄ species were determined to be 111.9885(26), 111.9817(46), and 124.2394(49) uŲ, respectively. The almost same values of P_bb of normal and d₁ species suggest that the hydroxyl hydrogen atom is very close to the a-c plane. From the r_s coordinates of the hydroxyl hydrogen atom, the OH bond was found to be at the trans position with respect to the CN double bond. The conformation of cyclohexanone oxime was determined to be chair form by comparing the observed and calculated rotational constants, ΔI, and planar moments, and taking account of the calculated the relative energy difference, ΔE. The structural parameters, the three bond lengths, three bond angles, and three dihedral angles, were adjusted to the nine rotational constants observed. The bond angle of ∠C₂C₁N is much wider than that of ∠C₆C₁N by about 10°. The dihedral angles of ∠C₁C₂C₃C₄, ∠C₂C₃C₄C₅, and ∠C₃C₄C₅C₆ were determined to be 53.3(5), −57.2(5), and 57.2(5)°. Two vibrational modes were assigned to the ring-bending and ring-twisting ones, which are almost harmonic up to v = 3.     

Phase-matched third-harmonic UV generation using low-order modes in a glass micro-fiber

We observe phase-matched third-harmonic generation at 355 nm in a low-order mode of a sub-micron diameter glass fiber. The third-harmonic signal exhibits a sharp resonance for a fiber diameter d = 0.49 ± 0.02 μm, in excellent agreement with the value d = 0.506 μm predicted by theory. The third-harmonic conversion efficiency is 2 × 10⁻⁶, and is limited by the pump power (1 kW) and effective fiber length (100 μm).