Techniques for Analysis of Certain Nonlinear Systems

This paper suggests the use of simple transformations like ÿ=kx, kx2 for second-order nonlinear differential equations to effect rapid plotting of the phase-plane trajectories. The method is particularly helpful in determining quickly the trajectory slopes along simple curves in any desired region of the phase plane. New planes such as the tÿ-x, tÿ2-x are considered for the study of some groups of nonlinear time-varying systems. Suggestions for solving certain higher-order nonlinear systems are also made.

Arbitrary Lagrangian-Eulerian finite-element method for computation of two-phase flows with soluble surfactants

A finite-element scheme based on a coupled arbitrary Lagrangian-Eulerian and Lagrangian approach is developed for the computation of interface flows with soluble surfactants. The numerical scheme is designed to solve the time-dependent Navier-Stokes equations and an evolution equation for the surfactant concentration in the bulk phase, and simultaneously, an evolution equation for the surfactant concentration on the interface. Second-order isoparametric finite elements on moving meshes and second-order isoparametric surface finite elements are used to solve these equations. The interface-resolved moving meshes allow the accurate incorporation of surface forces, Marangoni forces and jumps in the material parameters. The lower-dimensional finite-element meshes for solving the surface evolution equation are part of the interface-resolved moving meshes. The numerical scheme is validated for problems with known analytical solutions. A number of computations to study the influence of the surfactants in 3D-axisymmetric rising bubbles have been performed. The proposed scheme shows excellent conservation of fluid mass and of the total mass of the surfactant. (C) 2012 Elsevier Inc. All rights reserved.

A Numerical Scheme for Three-Dimensional Front Propagation and Control of Jordan Mode

As an example of a front propagation, we study the propagation of a three-dimensional nonlinear wavefront into a polytropic gas in a uniform state and at rest. The successive positions and geometry of the wavefront are obtained by solving the conservation form of equations of a weakly nonlinear ray theory. The proposed set of equations forms a weakly hyperbolic system of seven conservation laws with an additional vector constraint, each of whose components is a divergence-free condition. This constraint is an involution for the system of conservation laws, and it is termed a geometric solenoidal constraint. The analysis of a Cauchy problem for the linearized system shows that when this constraint is satisfied initially, the solution does not exhibit any Jordan mode. For the numerical simulation of the conservation laws we employ a high resolution central scheme. The second order accuracy of the scheme is achieved by using MUSCL-type reconstructions and Runge-Kutta time discretizations. A constrained transport-type technique is used to enforce the geometric solenoidal constraint. The results of several numerical experiments are presented, which confirm the efficiency and robustness of the proposed numerical method and the control of the Jordan mode.

Two-dimensional flow past circular cylinders using finite volume lattice Boltzmann formulation

In this article, an extension to the total variation diminishing finite volume formulation of the lattice Boltzmann equation method on unstructured meshes was presented. The quadratic least squares procedure is used for the estimation of first-order and second-order spatial gradients of the particle distribution functions. The distribution functions were extrapolated quadratically to the virtual upwind node. The time integration was performed using the fourth-order RungeKutta procedure. A grid convergence study was performed in order to demonstrate the order of accuracy of the present scheme. The formulation was validated for the benchmark two-dimensional, laminar, and unsteady flow past a single circular cylinder. These computations were then investigated for the low Mach number simulations. Further validation was performed for flow past two circular cylinders arranged in tandem and side-by-side. Results of these simulations were extensively compared with the previous numerical data. Copyright (C) 2011 John Wiley & Sons, Ltd.

Theoretical analysis of photoinduced H-atom elimination in thiophenol

The photoinduced hydrogen elimination reaction in thiophenol via the conical intersections of the dissociative (1)pi sigma* excited state with the bound (1)pi pi* excited state and the electronic ground state has been investigated with ab initio electronic-structure calculations and time-dependent quantum wave-packet calculations. A screening of the coupling constants of the symmetry-allowed coupling modes at the (1)pi pi*-(1)pi sigma* and (1)pi sigma*-S-0 conical intersection shows that the SH torsional mode is by far the most important coupling mode at both conical intersections. A model including three intersecting potential-energy surfaces (S-0, (1)pi pi*, (1)pi sigma*) and two nuclear degrees of freedom (SH stretch and SH torsion) has been constructed on the basis of ab initio complete-active-space self-consistent field and multireference second-order perturbation theory calculations. The nonadiabatic quantum wave-packet dynamics initiated by optical excitation of the (1)pi pi* and (1)pi sigma* states has been explored for this three-state two-coordinate model. The photodissociation dynamics is characterized in terms of snapshots of time-dependent wave packets, time-dependent electronic population probabilities, and the branching ratio of the (2)sigma/(2)pi electronic states of the thiophenoxyl radical. The dependence of the timescale of the photodissociation process and the branching ratio on the initial excitation of the SH stretching and SH torsional vibrations has been analyzed. It is shown that the node structure, which is imposed on the nuclear wave packets by the initial vibrational preparation as well as by the transitions through the conical intersections, has a profound effect on the photodissociation dynamics. The effect of additional weak coupling modes of CC twist (nu(16a)) and ring-distortion (nu(16b)) character has been investigated with three-dimensional and four-dimensional time-dependent wave-packet calculations, and has been found to be minor. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4709608]

Electrical Resistance and Seebeck Coefficient in PbTe Nanowires

We address a physics-based simplified analytical formulation of the diffusive electrical resistance ( (Omega)) and Seebeck coefficient () in a PbTe nanowire dominated by acoustic phonon scattering under the presence of a low static longitudinal electric field. The use of a second-order nonparabolic electron energy band structure involving a geometry-dependent band gap has been selected in principle to demonstrate that the electron mean free path (MFP) in such a system can reach as low as about 8 nm at room temperature for a 10-nm-wide PbTe nanowire. This is followed by the formulation of the carrier back-scattering coefficient for determination of (Omega) and as functions of wire dimensions, temperature, and the field, respectively. The present analytical formulation agrees well with the available experimental data and may find extensive use in determination of various electrothermal transport phenomena in PbTe-based one-dimensional electron devices.

A CMOS Gas Sensor Array Platform With Fourier Transform Based Impedance Spectroscopy

A CMOS gas sensor array platform with digital read-out containing 27 sensor pixels and a reference pixel is presented. A signal conditioning circuit at each pixel includes digitally programmable gain stages for sensor signal amplification followed by a second order continuous time delta sigma modulator for digitization. Each sensor pixel can be functionalized with a distinct sensing material that facilitates transduction based on impedance change. Impedance spectrum (up to 10 KHz) of the sensor is obtained off-chip by computing the fast Fourier transform of sensor and reference pixel outputs. The reference pixel also compensates for the phase shift introduced by the signal processing circuits. The chip also contains a temperature sensor with digital readout for ambient temperature measurement. A sensor pixel is functionalized with polycarbazole conducting polymer for sensing volatile organic gases and measurement results are presented. The chip is fabricated in a 0.35 CMOS technology and requires a single step post processing for functionalization. It consumes 57 mW from a 3.3 V supply.

Directional Diffusion Regulator (DDR) for some numerical solvers of hyperbolic conservation laws

A computational tool called ``Directional Diffusion Regulator (DDR)'' is proposed to bring forth real multidimensional physics into the upwind discretization in some numerical schemes of hyperbolic conservation laws. The direction based regulator when used with dimension splitting solvers, is set to moderate the excess multidimensional diffusion and hence cause genuine multidimensional upwinding like effect. The basic idea of this regulator driven method is to retain a full upwind scheme across local discontinuities, with the upwind bias decreasing smoothly to a minimum in the farthest direction. The discontinuous solutions are quantified as gradients and the regulator parameter across a typical finite volume interface or a finite difference interpolation point is formulated based on fractional local maximum gradient in any of the weak solution flow variables (say density, pressure, temperature, Mach number or even wave velocity etc.). DDR is applied to both the non-convective as well as whole unsplit dissipative flux terms of some numerical schemes, mainly of Local Lax-Friedrichs, to solve some benchmark problems describing inviscid compressible flow, shallow water dynamics and magneto-hydrodynamics. The first order solutions consistently improved depending on the extent of grid non-alignment to discontinuities, with the major influence due to regulation of non-convective diffusion. The application is also experimented on schemes such as Roe, Jameson-Schmidt-Turkel and some second order accurate methods. The consistent improvement in accuracy either at moderate or marked levels, for a variety of problems and with increasing grid size, reasonably indicate a scope for DDR as a regular tool to impart genuine multidimensional upwinding effect in a simpler framework. (C) 2012 Elsevier Inc. All rights reserved.

Lithium D-Isoascorbate Monohydrate, a New Nonlinear Optical Material

Single crystals of lithium D-isoascorbate monohydrate (LDAM), (C6H7O6Li center dot H2O), are grown by a solution growth method. The crystal structure of LDAM is solved using single crystal X-ray diffraction. The space group is orthorhombic P2(1)2(1)2(1) with four formula units per unit cell and lattice parameters a = 7.7836(3) angstrom, b = 8.7456(3) angstrom, and c = 11.0368(4) angstrom. Solubility of the material in water is determined thermogravimetrically and found to have a positive temperature coefficient of solubility. Large optical quality single crystals are subsequently grown from aqueous solution by a slow cooling method. The crystal has a bulky prismatic habit and among the prominent faces the c face appears as the only principal morphological face. The crystal exhibits a (010) cleavage. Dielectric spectroscopy reveals a nearly Debye type Cole-Cole behavior with anisotropy in relaxation. Optical transmission range is found to be from 300 to 1400 nm. The principal refractive indices of this biaxial crystal, measured using Brewster's angle method, at wavelengths 405, 543, and 632.8 nm, show high dispersion. The crystal is negative biaxial with 2V(z) = 107.8 degrees (405 nm) and belongs to the Hobden class 3. Theoretically generated type 1 and type 2 second order phase matching curves match very well with the experimental results. The second-order nonlinear coefficient d(14) was determined to be 7 x 10(-13) m/V. For the optimum phase matching direction (type 2), the second-order effective nonlinear coefficient and the walk off angle are determined to be 0.84 times d(14) and 3.5 degrees respectively. The crystal possesses high multiple surface damage thresholds of 18 GW/cm(2) and 8 GW/cm(2) at laser wavelengths 1064 and 532 nm, respectively.

A physics-based flexural phonon-dependent thermal conductivity model for single layer graphene

In this paper, we address a physics-based closed-form analytical model of flexural phonon-dependent diffusive thermal conductivity (kappa) of suspended rectangular single layer graphene sheet. A quadratic dependence of the out-of-plane phonon frequency, generally called flexural phonons, on the phonon wave vector has been taken into account to analyze the behavior of kappa at lower temperatures. Such a dependence has further been used for the determination of second-order three-phonon Umklapp and isotopic scatterings. We find that these behaviors in our model are best explained through the upper limit of Debye cut-off frequency in the second-order three-phonon Umklapp scattering of the long phonon waves that actually remove the thermal conductivity singularity by contributing a constant scattering rate at low frequencies and note that the out-of-plane Gruneisen parameter for these modes need not be too high. Using this, we clearly demonstrate that. follows a T-1.5 and T-2 law at lower and higher temperatures in the absence of isotopes, respectively. However in their presence, the behavior of kappa sharply deviates from the T-2 law at higher temperatures. The present geometry-dependent model of kappa is found to possess an excellent match with various experimental data over a wide range of temperatures which can be put forward for efficient electro-thermal analyses of encased/supported graphene.

Observation of the nuclear magnetic octupole moment of Yb-173 from precise measurements of the hyperfine structure in the P-3(2) state

We measure hyperfine structure in the metastable P-3(2) state of Yb-173 and extract the nuclear magnetic octupole moment. We populate the state using dipole-allowed transitions through the P-3(1) and S-3(1) states. We measure frequencies of hyperfine transitions of the P-3(2) -> S-3(1) line at 770 nm using a Rb-stabilized ring cavity resonator with a precision of 200 kHz. Second-order corrections due to perturbations from the nearby P-3(1) and P-1(1) states are below 30 kHz. We obtain the hyperfine coefficients as A = -742.11(2) MHz and B = 1339.2(2) MHz, which represent a two orders-of-magnitude improvement in precision, and C = 0.54(2) MHz. From atomic structure calculations, we obtain the nuclear moments quadrupole Q = 2.46(12) b and octupole Omega = -34.4(21) b x mu(N). DOI: 10.1103/PhysRevA.87.012512

On the treewidth of MDS and Reed-Muller codes

The treewidth of a linear code is the least constraint complexity of any of its cycle-free graphical realizations. This notion provides a useful parametrization of the maximum-likelihood decoding complexity for linear codes. In this paper, we compute exact expressions for the treewidth of maximum distance separable codes, and first- and second-order Reed-Muller codes. These results constitute the only known explicit expressions for the treewidth of algebraic codes.

A finite-rate-of-innovation signal sampling perspective of the source-filter model of speech production

We consider the speech production mechanism and the asso- ciated linear source-filter model. For voiced speech sounds in particular, the source/glottal excitation is modeled as a stream of impulses and the filter as a cascade of second-order resonators. We show that the process of sampling speech signals can be modeled as filtering a stream of Dirac impulses (a model for the excitation) with a kernel function (the vocal tract response),and then sampling uniformly. We show that the problem of esti- mating the excitation is equivalent to the problem of recovering a stream of Dirac impulses from samples of a filtered version. We present associated algorithms based on the annihilating filter and also make a comparison with the classical linear prediction technique, which is well known in speech analysis. Results on synthesized as well as natural speech data are presented.

Monaural sound source separation using covariance profile of partials

This paper addresses the problem of separation of pitched sounds in monaural recordings. We present a novel feature for the estimation of parameters of overlapping harmonics which considers the covariance of partials of pitched sounds. Sound templates are formed from the monophonic parts of the mixture recording. A match for every note is found among these templates on the basis of covariance profile of their harmonics. The matching template for the note provides the second order characteristics for the overlapped harmonics of the note. The algorithm is tested on the RWC music database instrument sounds. The results clearly show that the covariance characteristics can be used to reconstruct overlapping harmonics effectively.

Narrowing of resonances in electromagnetically induced transparency and absorption using a Laguerre-Gaussian control beam

We study the phenomenon of electromagnetically induced transparency and absorption (EITA) using a control laser with a Laguerre-Gaussian (LG) profile instead of the usual Gaussian profile, and observe significant narrowing of the resonance widths. Aligning the probe beam to the central hole in the doughnut-shaped LG control beam allows simultaneously a strong control intensity required for high signal-to-noise ratio and a low intensity in the probe region required to get narrow resonances. Experiments with an expanded Gaussian control and a second-order LG control show that transit time and orbital angular momentum do not play a significant role. This explanation is borne out by a density-matrix analysis with a radially varying control Rabi frequency. We observe these resonances using degenerate two-level transitions in the D-2 line of Rb-87 in a room temperature vapor cell, and an EIA resonance with width up to 20 times below the natural linewidth for the F = 2 -> F' = 3 transition. Thus the use of LG beams should prove advantageous in all applications of EITA and other kinds of pump-probe spectroscopy as well.