Diffusion of interacting particles in discrete geometries: Equilibrium and dynamical properties

NASA Astrophysics Data System (ADS)

Becker, T.; Nelissen, K.; Cleuren, B.; Partoens, B.; Van den Broeck, C.

2014-11-01

We expand on a recent study of a lattice model of interacting particles [Phys. Rev. Lett. 111, 110601 (2013), 10.1103/PhysRevLett.111.110601]. The adsorption isotherm and equilibrium fluctuations in particle number are discussed as a function of the interaction. Their behavior is similar to that of interacting particles in porous materials. Different expressions for the particle jump rates are derived from transition-state theory. Which expression should be used depends on the strength of the interparticle interactions. Analytical expressions for the self- and transport diffusion are derived when correlations, caused by memory effects in the environment, are neglected. The diffusive behavior is studied numerically with kinetic Monte Carlo (kMC) simulations, which reproduces the diffusion including correlations. The effect of correlations is studied by comparing the analytical expressions with the kMC simulations. It is found that the Maxwell-Stefan diffusion can exceed the self-diffusion. To our knowledge, this is the first time this is observed. The diffusive behavior in one-dimensional and higher-dimensional systems is qualitatively the same, with the effect of correlations decreasing for increasing dimension. The length dependence of both the self- and transport diffusion is studied for one-dimensional systems. For long lengths the self-diffusion shows a 1 /L dependence. Finally, we discuss when agreement with experiments and simulations can be expected. The assumption that particles in different cavities do not interact is expected to hold quantitatively at low and medium particle concentrations if the particles are not strongly interacting.

Analytical model for advective-dispersive transport involving flexible boundary inputs, initial distributions and zero-order productions

NASA Astrophysics Data System (ADS)

Chen, Jui-Sheng; Li, Loretta Y.; Lai, Keng-Hsin; Liang, Ching-Ping

2017-11-01

A novel solution method is presented which leads to an analytical model for the advective-dispersive transport in a semi-infinite domain involving a wide spectrum of boundary inputs, initial distributions, and zero-order productions. The novel solution method applies the Laplace transform in combination with the generalized integral transform technique (GITT) to obtain the generalized analytical solution. Based on this generalized analytical expression, we derive a comprehensive set of special-case solutions for some time-dependent boundary distributions and zero-order productions, described by the Dirac delta, constant, Heaviside, exponentially-decaying, or periodically sinusoidal functions as well as some position-dependent initial conditions and zero-order productions specified by the Dirac delta, constant, Heaviside, or exponentially-decaying functions. The developed solutions are tested against an analytical solution from the literature. The excellent agreement between the analytical solutions confirms that the new model can serve as an effective tool for investigating transport behaviors under different scenarios. Several examples of applications, are given to explore transport behaviors which are rarely noted in the literature. The results show that the concentration waves resulting from the periodically sinusoidal input are sensitive to dispersion coefficient. The implication of this new finding is that a tracer test with a periodic input may provide additional information when for identifying the dispersion coefficients. Moreover, the solution strategy presented in this study can be extended to derive analytical models for handling more complicated problems of solute transport in multi-dimensional media subjected to sequential decay chain reactions, for which analytical solutions are not currently available.

How cosmic microwave background correlations at large angles relate to mass autocorrelations in space

NASA Technical Reports Server (NTRS)

Blumenthal, George R.; Johnston, Kathryn V.

1994-01-01

The Sachs-Wolfe effect is known to produce large angular scale fluctuations in the cosmic microwave background radiation (CMBR) due to gravitational potential fluctuations. We show how the angular correlation function of the CMBR can be expressed explicitly in terms of the mass autocorrelation function xi(r) in the universe. We derive analytic expressions for the angular correlation function and its multipole moments in terms of integrals over xi(r) or its second moment, J(sub 3)(r), which does not need to satisfy the sort of integral constraint that xi(r) must. We derive similar expressions for bulk flow velocity in terms of xi and J(sub 3). One interesting result that emerges directly from this analysis is that, for all angles theta, there is a substantial contribution to the correlation function from a wide range of distance r and that radial shape of this contribution does not vary greatly with angle.

ODF Maxima Extraction in Spherical Harmonic Representation via Analytical Search Space Reduction

PubMed Central

Aganj, Iman; Lenglet, Christophe; Sapiro, Guillermo

2015-01-01

By revealing complex fiber structure through the orientation distribution function (ODF), q-ball imaging has recently become a popular reconstruction technique in diffusion-weighted MRI. In this paper, we propose an analytical dimension reduction approach to ODF maxima extraction. We show that by expressing the ODF, or any antipodally symmetric spherical function, in the common fourth order real and symmetric spherical harmonic basis, the maxima of the two-dimensional ODF lie on an analytically derived one-dimensional space, from which we can detect the ODF maxima. This method reduces the computational complexity of the maxima detection, without compromising the accuracy. We demonstrate the performance of our technique on both artificial and human brain data. PMID:20879302

Coupled-cluster Green's function: Analysis of properties originating in the exponential parametrization of the ground-state wave function

DOE Office of Scientific and Technical Information (OSTI.GOV)

Peng, Bo; Kowalski, Karol

In this paper we derive basic properties of the Green’s function matrix elements stemming from the exponential coupled cluster (CC) parametrization of the ground-state wave function. We demon- strate that all intermediates used to express retarded (or equivalently, ionized) part of the Green’s function in the ω-representation can be expressed through connected diagrams only. Similar proper- ties are also shared by the first order ω-derivatives of the retarded part of the CC Green’s function. This property can be extended to any order ω-derivatives of the Green’s function. Through the Dyson equation of CC Green’s function, the derivatives of corresponding CCmore » self-energy can be evaluated analytically. In analogy to the CC Green’s function, the corresponding CC self-energy is expressed in terms of connected diagrams only. Moreover, the ionized part of the CC Green’s func- tion satisfies the non-homogeneous linear system of ordinary differential equations, whose solution may be represented in the exponential form. Our analysis can be easily generalized to the advanced part of the CC Green’s function.« less

Some analytical models to estimate maternal age at birth using age-specific fertility rates.

PubMed

Pandey, A; Suchindran, C M

1995-01-01

"A class of analytical models to study the distribution of maternal age at different births from the data on age-specific fertility rates has been presented. Deriving the distributions and means of maternal age at birth of any specific order, final parity and at next-to-last birth, we have extended the approach to estimate parity progression ratios and the ultimate parity distribution of women in the population.... We illustrate computations of various components of the model expressions with the current fertility experiences of the United States for 1970." excerpt

The Effect of Planetary Albedo on Solar Orientation of Spacecraft

NASA Technical Reports Server (NTRS)

Fontana, Anthony

1967-01-01

The analytical expression for the solar orientation error caused by planetary albedo is derived. A typical solar sensor output characteristic is assumed and a computer solution to the analytical is obtained. The computer results are presented for a spacecraft in the vicinity of Earth, Venus, Mars, and the Moon. Each planetary body is assumed to be a spherical diffuse reflector with cylindrical shadows and a constant albedo. The data generated herein permit the selection of an appropriate coarse-sensor to fine-sensor switching angle for solar orientation control systems and facilitate the the interpretation of solar-referenced scientific experiment data.

Radiation force of an arbitrary acoustic beam on an elastic sphere in a fluid

PubMed Central

Sapozhnikov, Oleg A.; Bailey, Michael R.

2013-01-01

A theoretical approach is developed to calculate the radiation force of an arbitrary acoustic beam on an elastic sphere in a liquid or gas medium. First, the incident beam is described as a sum of plane waves by employing conventional angular spectrum decomposition. Then, the classical solution for the scattering of a plane wave from an elastic sphere is applied for each plane-wave component of the incident field. The net scattered field is expressed as a superposition of the scattered fields from all angular spectrum components of the incident beam. With this formulation, the incident and scattered waves are superposed in the far field to derive expressions for components of the radiation stress tensor. These expressions are then integrated over a spherical surface to analytically describe the radiation force on an elastic sphere. Limiting cases for particular types of incident beams are presented and are shown to agree with known results. Finally, the analytical expressions are used to calculate radiation forces associated with two specific focusing transducers. PMID:23363086

Exact results in the large system size limit for the dynamics of the chemical master equation, a one dimensional chain of equations.

PubMed

Martirosyan, A; Saakian, David B

2011-08-01

We apply the Hamilton-Jacobi equation (HJE) formalism to solve the dynamics of the chemical master equation (CME). We found exact analytical expressions (in large system-size limit) for the probability distribution, including explicit expression for the dynamics of variance of distribution. We also give the solution for some simple cases of the model with time-dependent rates. We derived the results of the Van Kampen method from the HJE approach using a special ansatz. Using the Van Kampen method, we give a system of ordinary differential equations (ODEs) to define the variance in a two-dimensional case. We performed numerics for the CME with stationary noise. We give analytical criteria for the disappearance of bistability in the case of stationary noise in one-dimensional CMEs.

An Analysis of the Effect of Surface Heat Exchange on the Thermal Behavior of an Idealized Aquifer Thermal Energy Storage System

NASA Astrophysics Data System (ADS)

Güven, O.; Melville, J. G.; Molz, F. J.

1983-06-01

Analytical expressions are derived for the temperature distribution and the mean temperature of an idealized aquifer thermal energy storage (ATES) system, taking into account the heat exchange at the ground surface and the finite thickness of the overlying layer above the storage aquifer. The analytical expressions for the mean temperature may be used to obtain rough estimates of first-cycle recovery factors for preliminary evaluations of shallow confined or unconfined ATES systems. The results, which are presented in nondimensional plots, indicate that surface heat exchange may have a significant influence on the thermal behavior of shallow ATES systems. Thus it is suggested that the effects of surface heat exchange should be considered carefully and included in the detailed analyses of such ATES systems.

Analytical pair correlations in ideal quantum gases: temperature-dependent bunching and antibunching.

PubMed

Bosse, J; Pathak, K N; Singh, G S

2011-10-01

The fluctuation-dissipation theorem together with the exact density response spectrum for ideal quantum gases has been utilized to yield a new expression for the static structure factor, which we use to derive exact analytical expressions for the temperature-dependent pair distribution function g(r) of the ideal gases. The plots of bosonic and fermionic g(r) display "Bose pile" and "Fermi hole" typically akin to bunching and antibunching as observed experimentally for ultracold atomic gases. The behavior of spin-scaled pair correlation for fermions is almost featureless, but bosons show a rich structure including long-range correlations near T(c). The coherent state at T=0 shows no correlation at all, just like single-mode lasers. The depicted decreasing trend in correlation with decrease in temperature for T

An alternative model for a partially coherent elliptical dark hollow beam

NASA Astrophysics Data System (ADS)

Li, Xu; Wang, Fei; Cai, Yangjian

2011-04-01

An alternative theoretical model named partially coherent hollow elliptical Gaussian beam (HEGB) is proposed to describe a partially coherent beam with an elliptical dark hollow profile. Explicit expression for the propagation factors of a partially coherent HEGB is derived. Based on the generalized Collins formula, analytical formulae for the cross-spectral density and mean-squared beam width of a partially coherent HEGB, propagating through a paraxial ABCD optical system, are derived. Propagation properties of a partially coherent HEGB in free space are studied as a numerical example.

A new method to compute lunisolar perturbations in satellite motions

NASA Technical Reports Server (NTRS)

Kozai, Y.

1973-01-01

A new method to compute lunisolar perturbations in satellite motion is proposed. The disturbing function is expressed by the orbital elements of the satellite and the geocentric polar coordinates of the moon and the sun. The secular and long periodic perturbations are derived by numerical integrations, and the short periodic perturbations are derived analytically. The perturbations due to the tides can be included in the same way. In the Appendix, the motion of the orbital plane for a synchronous satellite is discussed; it is concluded that the inclination cannot stay below 7 deg.

Analytical Model and Optimized Design of Power Transmitting Coil for Inductively Coupled Endoscope Robot.

PubMed

Ke, Quan; Luo, Weijie; Yan, Guozheng; Yang, Kai

2016-04-01

A wireless power transfer system based on the weakly inductive coupling makes it possible to provide the endoscope microrobot (EMR) with infinite power. To facilitate the patients' inspection with the EMR system, the diameter of the transmitting coil is enlarged to 69 cm. Due to the large transmitting range, a high quality factor of the Litz-wire transmitting coil is a necessity to ensure the intensity of magnetic field generated efficiently. Thus, this paper builds an analytical model of the transmitting coil, and then, optimizes the parameters of the coil by enlarging the quality factor. The lumped model of the transmitting coil includes three parameters: ac resistance, self-inductance, and stray capacitance. Based on the exact two-dimension solution, the accurate analytical expression of ac resistance is derived. Several transmitting coils of different specifications are utilized to verify this analytical expression, being in good agreements with the measured results except the coils with a large number of strands. Then, the quality factor of transmitting coils can be well predicted with the available analytical expressions of self- inductance and stray capacitance. Owing to the exact estimation of quality factor, the appropriate coil turns of the transmitting coil is set to 18-40 within the restrictions of transmitting circuit and human tissue issues. To supply enough energy for the next generation of the EMR equipped with a Ø9.5×10.1 mm receiving coil, the coil turns of the transmitting coil is optimally set to 28, which can transfer a maximum power of 750 mW with the remarkable delivering efficiency of 3.55%.

Thermostatistical description of gas mixtures from space partitions

NASA Astrophysics Data System (ADS)

Rohrmann, R. D.; Zorec, J.

2006-10-01

The new mathematical framework based on the free energy of pure classical fluids presented by Rohrmann [Physica A 347, 221 (2005)] is extended to multicomponent systems to determine thermodynamic and structural properties of chemically complex fluids. Presently, the theory focuses on D -dimensional mixtures in the low-density limit (packing factor η<0.01 ). The formalism combines the free-energy minimization technique with space partitions that assign an available volume v to each particle. v is related to the closeness of the nearest neighbor and provides a useful tool to evaluate the perturbations experimented by particles in a fluid. The theory shows a close relationship between statistical geometry and statistical mechanics. New, unconventional thermodynamic variables and mathematical identities are derived as a result of the space division. Thermodynamic potentials μil , conjugate variable of the populations Nil of particles class i with the nearest neighbors of class l are defined and their relationships with the usual chemical potentials μi are established. Systems of hard spheres are treated as illustrative examples and their thermodynamics functions are derived analytically. The low-density expressions obtained agree nicely with those of scaled-particle theory and Percus-Yevick approximation. Several pair distribution functions are introduced and evaluated. Analytical expressions are also presented for hard spheres with attractive forces due to Kac-tails and square-well potentials. Finally, we derive general chemical equilibrium conditions.

Computing sensitivity and selectivity in parallel factor analysis and related multiway techniques: the need for further developments in net analyte signal theory.

PubMed

Olivieri, Alejandro C

2005-08-01

Sensitivity and selectivity are important figures of merit in multiway analysis, regularly employed for comparison of the analytical performance of methods and for experimental design and planning. They are especially interesting in the second-order advantage scenario, where the latter property allows for the analysis of samples with a complex background, permitting analyte determination even in the presence of unsuspected interferences. Since no general theory exists for estimating the multiway sensitivity, Monte Carlo numerical calculations have been developed for estimating variance inflation factors, as a convenient way of assessing both sensitivity and selectivity parameters for the popular parallel factor (PARAFAC) analysis and also for related multiway techniques. When the second-order advantage is achieved, the existing expressions derived from net analyte signal theory are only able to adequately cover cases where a single analyte is calibrated using second-order instrumental data. However, they fail for certain multianalyte cases, or when third-order data are employed, calling for an extension of net analyte theory. The results have strong implications in the planning of multiway analytical experiments.

Modification of the parallel scattering mean free path of cosmic rays in the presence of adiabatic focusing

DOE Office of Scientific and Technical Information (OSTI.GOV)

He, H.-Q.; Schlickeiser, R., E-mail: hqhe@mail.iggcas.ac.cn, E-mail: rsch@tp4.rub.de

The cosmic ray mean free path in a large-scale nonuniform guide magnetic field with superposed magnetostatic turbulence is calculated to clarify some conflicting results in the literature. A new, exact integro-differential equation for the cosmic-ray anisotropy is derived from the Fokker-Planck transport equation. A perturbation analysis of this integro-differential equation leads to an analytical expression for the cosmic ray anisotropy and the focused transport equation for the isotropic part of the cosmic ray distribution function. The derived parallel spatial diffusion coefficient and the associated cosmic ray mean free path include the effect of adiabatic focusing and reduce to the standardmore » forms in the limit of a uniform guide magnetic field. For the illustrative case of isotropic pitch angle scattering, the derived mean free path agrees with the earlier expressions of Beeck and Wibberenz, Bieber and Burger, Kota, and Litvinenko, but disagrees with the result of Shalchi. The disagreement with the expression of Shalchi is particularly strong in the limit of strong adiabatic focusing.« less

Simplified expressions that incorporate finite pulse effects into coherent two-dimensional optical spectra.

PubMed

Do, Thanh Nhut; Gelin, Maxim F; Tan, Howe-Siang

2017-10-14

We derive general expressions that incorporate finite pulse envelope effects into a coherent two-dimensional optical spectroscopy (2DOS) technique. These expressions are simpler and less computationally intensive than the conventional triple integral calculations needed to simulate 2DOS spectra. The simplified expressions involving multiplications of arbitrary pulse spectra with 2D spectral response function are shown to be exactly equal to the conventional triple integral calculations of 2DOS spectra if the 2D spectral response functions do not vary with population time. With minor modifications, they are also accurate for 2D spectral response functions with quantum beats and exponential decay during population time. These conditions cover a broad range of experimental 2DOS spectra. For certain analytically defined pulse spectra, we also derived expressions of 2D spectra for arbitrary population time dependent 2DOS spectral response functions. Having simpler and more efficient methods to calculate experimentally relevant 2DOS spectra with finite pulse effect considered will be important in the simulation and understanding of the complex systems routinely being studied by using 2DOS.

Velocity space scattering coefficients with applications in antihydrogen recombination studies

NASA Astrophysics Data System (ADS)

Chang, Yongbin; Ordonez, C. A.

2000-12-01

An approach for calculating velocity space friction and diffusion coefficients with Maxwellian field particles is developed based on a kernel function derived in a previous paper [Y. Chang and C. A. Ordonez, Phys. Plasmas 6, 2947 (1999)]. The original fivefold integral expressions for the coefficients are reduced to onefold integrals, which can be used for any value of the Coulomb logarithm. The onefold integrals can be further reduced to standard analytical expressions by using a weak coupling approximation. The integral expression for the friction coefficient is used to predict a time scale that describes the rate at which a reflecting antiproton beam slows down within a positron plasma, while both species are simultaneously confined by a nested Penning trap. The time scale is used to consider the possibility of achieving antihydrogen recombination within the trap. The friction and diffusion coefficients are then used to derive an expression for calculating the energy transfer rate between antiprotons and positrons. The expression is employed to illustrate achieving antihydrogen recombination while taking into account positron heating by the antiprotons. The effect of the presence of an electric field on recombination is discussed.

Closed-form recursive formula for an optimal tracker with terminal constraints

NASA Technical Reports Server (NTRS)

Juang, J.-N.; Turner, J. D.; Chun, H. M.

1984-01-01

Feedback control laws are derived for a class of optimal finite time tracking problems with terminal constraints. Analytical solutions are obtained for the feedback gain and the closed-loop response trajectory. Such formulations are expressed in recursive forms so that a real-time computer implementation becomes feasible. Two examples are given to illustrate the validity and usefulness of the formulations.

Simulation of inverse Compton scattering and its implications on the scattered linewidth

NASA Astrophysics Data System (ADS)

Ranjan, N.; Terzić, B.; Krafft, G. A.; Petrillo, V.; Drebot, I.; Serafini, L.

2018-03-01

Rising interest in inverse Compton sources has increased the need for efficient models that properly quantify the behavior of scattered radiation given a set of interaction parameters. The current state-of-the-art simulations rely on Monte Carlo-based methods, which, while properly expressing scattering behavior in high-probability regions of the produced spectra, may not correctly simulate such behavior in low-probability regions (e.g. tails of spectra). Moreover, sampling may take an inordinate amount of time for the desired accuracy to be achieved. In this paper, we present an analytic derivation of the expression describing the scattered radiation linewidth and propose a model to describe the effects of horizontal and vertical emittance on the properties of the scattered radiation. We also present an improved version of the code initially reported in Krafft et al. [Phys. Rev. Accel. Beams 19, 121302 (2016), 10.1103/PhysRevAccelBeams.19.121302], that can perform the same simulations as those present in cain and give accurate results in low-probability regions by integrating over the emissions of the electrons. Finally, we use these codes to carry out simulations that closely verify the behavior predicted by the analytically derived scaling law.

Simulation of inverse Compton scattering and its implications on the scattered linewidth

DOE PAGES

Ranjan, N.; Terzić, B.; Krafft, G. A.; ...

2018-03-06

Rising interest in inverse Compton sources has increased the need for efficient models that properly quantify the behavior of scattered radiation given a set of interaction parameters. The current state-of-the-art simulations rely on Monte Carlo-based methods, which, while properly expressing scattering behavior in high-probability regions of the produced spectra, may not correctly simulate such behavior in low-probability regions (e.g. tails of spectra). Moreover, sampling may take an inordinate amount of time for the desired accuracy to be achieved. Here in this article, we present an analytic derivation of the expression describing the scattered radiation linewidth and propose a model tomore » describe the effects of horizontal and vertical emittance on the properties of the scattered radiation. We also present an improved version of the code initially reported in Krafft et al. [Phys. Rev. Accel. Beams 19, 121302 (2016)], that can perform the same simulations as those present in cain and give accurate results in low-probability regions by integrating over the emissions of the electrons. Finally, we use these codes to carry out simulations that closely verify the behavior predicted by the analytically derived scaling law.« less

Solving Kinetic Equations for the Laser Flash Photolysis Experiment on Nitric Oxide Synthases: Effect of Conformational Dynamics on the Interdomain Electron Transfer.

PubMed

Astashkin, Andrei V; Feng, Changjian

2015-11-12

The production of nitric oxide by the nitric oxide synthase (NOS) enzyme depends on the interdomain electron transfer (IET) between the flavin mononucleotide (FMN) and heme domains. Although the rate of this IET has been measured by laser flash photolysis (LFP) for various NOS proteins, no rigorous analysis of the relevant kinetic equations was performed so far. In this work, we provide an analytical solution of the kinetic equations underlying the LFP approach. The derived expressions reveal that the bulk IET rate is significantly affected by the conformational dynamics that determines the formation and dissociation rates of the docking complex between the FMN and heme domains. We show that in order to informatively study the electron transfer across the NOS enzyme, LFP should be used in combination with other spectroscopic methods that could directly probe the docking equilibrium and the conformational change rate constants. The implications of the obtained analytical expressions for the interpretation of the LFP results from various native and modified NOS proteins are discussed. The mathematical formulas derived in this work should also be applicable for interpreting the IET kinetics in other modular redox enzymes.

Effective refractive index and first-order-mode cutoff conditions in InGaAsP/InP DH laser structures /lambda = 1.2-1.6 microns/

NASA Technical Reports Server (NTRS)

Botez, D.

1982-01-01

A highly accurate analytical expression for the effective refractive index in In GaAsP/InP DH lasers emitting in the 1.2-1.6 micron range is presented. This closed-form expression is used to derive simple wavelength-independent expressions for the first-order mode cutoff conditions of various lateral waveguides. The effective refractive index is a function of emission wavelength and active layer thickness, and the mode cutoff conditions are compared to experimental data from mode-stabilized 1.3 and 1.55 micron DH lasers.

Gene Profiling in Experimental Models of Eye Growth: Clues to Myopia Pathogenesis

PubMed Central

Stone, Richard A.; Khurana, Tejvir S.

2010-01-01

To understand the complex regulatory pathways that underlie the development of refractive errors, expression profiling has evaluated gene expression in ocular tissues of well-characterized experimental models that alter postnatal eye growth and induce refractive errors. Derived from a variety of platforms (e.g. differential display, spotted microarrays or Affymetrix GeneChips), gene expression patterns are now being identified in species that include chicken, mouse and primate. Reconciling available results is hindered by varied experimental designs and analytical/statistical features. Continued application of these methods offers promise to provide the much-needed mechanistic framework to develop therapies to normalize refractive development in children. PMID:20363242

Optimization of cascade blade mistuning. I - Equations of motion and basic inherent properties

NASA Technical Reports Server (NTRS)

Nissim, E.

1985-01-01

Attention is given to the derivation of the equations of motion of mistuned compressor blades, interpolating aerodynamic coefficients by means of quadratic expressions in the reduced frequency. If the coefficients of the quadratic expressions are permitted to assume complex values, excellent accuracy is obtained and Pade rational expressions are obviated. On the basis of the resulting equations, it is shown analytically that the sum of all the real parts of the eigenvalues is independent of the mistuning introduced into the system. Blade mistuning is further treated through the aerodynamic energy approach, and the limiting vibration modes associated with alternative mistunings are identified.

Loading-unloading response of circular GLARE fiber-metal laminates under lateral indentation

NASA Astrophysics Data System (ADS)

Tsamasphyros, George J.; Bikakis, George S.

2015-01-01

GLARE is a Fiber-Metal laminated material used in aerospace structures which are frequently subjected to various impact damages. Hence, the response of GLARE plates subjected to lateral indentation is very important. In this paper, analytical expressions are derived and a non-linear finite element modeling procedure is proposed in order to predict the static load-indentation curves of circular GLARE plates during loading and unloading by a hemispherical indentor. We have recently published analytical formulas and a finite element procedure for the static indentation of circular GLARE plates which are now used during the loading stage. Here, considering that aluminum layers are in a state of membrane yield and employing energy balance during unloading, the unloading path is determined. Using this unloading path, an algebraic equation is derived for calculating the permanent dent depth of the GLARE plate after the indentor's withdrawal. Furthermore, our finite element procedure is modified in order to simulate the unloading stage as well. The derived formulas and the proposed finite element modeling procedure are applied successfully to GLARE 2-2/1-0.3 and to GLARE 3-3/2-0.4 circular plates. The analytical results are compared with corresponding FEM results and a good agreement is found. The analytically calculated permanent dent depth is within 6 % for the GLARE 2 plate, and within 7 % for the GLARE 3 plate, of the corresponding numerically calculated result. No other solution of this problem is known to the authors.

Dynamic NMR under nonstationary conditions: Theoretical model, numerical calculation, and potential of application

DOE Office of Scientific and Technical Information (OSTI.GOV)

Babailov, S. P., E-mail: babajlov@niic.nsc.ru; National Research Tomsk Polytechnic University, Lenin Avenue 30, Tomsk 634050; Purtov, P. A.

An expression has been derived for the time dependence of the NMR line shape for systems with multi-site chemical exchange in the absence of spin-spin coupling, in a zero saturation limit. The dynamics of variation of the NMR line shape with time is considered in detail for the case of two-site chemical exchange. Mathematical programs have been designed for numerical simulation of the NMR spectra of chemical exchange systems. The analytical expressions obtained are useful for NMR line shape simulations for systems with photoinduced chemical exchange.

Transition probability functions for applications of inelastic electron scattering

PubMed Central

Löffler, Stefan; Schattschneider, Peter

2012-01-01

In this work, the transition matrix elements for inelastic electron scattering are investigated which are the central quantity for interpreting experiments. The angular part is given by spherical harmonics. For the weighted radial wave function overlap, analytic expressions are derived in the Slater-type and the hydrogen-like orbital models. These expressions are shown to be composed of a finite sum of polynomials and elementary trigonometric functions. Hence, they are easy to use, require little computation time, and are significantly more accurate than commonly used approximations. PMID:22560709

Energy density and energy flow of surface waves in a strongly magnetized graphene

NASA Astrophysics Data System (ADS)

Moradi, Afshin

2018-01-01

General expressions for the energy density and energy flow of plasmonic waves in a two-dimensional massless electron gas (as a simple model of graphene) are obtained by means of the linearized magneto-hydrodynamic model and classical electromagnetic theory when a strong external magnetic field perpendicular to the system is present. Also, analytical expressions for the energy velocity, wave polarization, wave impedance, transverse and longitudinal field strength functions, and attenuation length of surface magneto-plasmon-polariton waves are derived, and numerical results are prepared.

The impact of inter-fraction dose variations on biological equivalent dose (BED): the concept of equivalent constant dose.

PubMed

Zavgorodni, S

2004-12-07

Inter-fraction dose fluctuations, which appear as a result of setup errors, organ motion and treatment machine output variations, may influence the radiobiological effect of the treatment even when the total delivered physical dose remains constant. The effect of these inter-fraction dose fluctuations on the biological effective dose (BED) has been investigated. Analytical expressions for the BED accounting for the dose fluctuations have been derived. The concept of biological effective constant dose (BECD) has been introduced. The equivalent constant dose (ECD), representing the constant physical dose that provides the same cell survival fraction as the fluctuating dose, has also been introduced. The dose fluctuations with Gaussian as well as exponential probability density functions were investigated. The values of BECD and ECD calculated analytically were compared with those derived from Monte Carlo modelling. The agreement between Monte Carlo modelled and analytical values was excellent (within 1%) for a range of dose standard deviations (0-100% of the dose) and the number of fractions (2 to 37) used in the comparison. The ECDs have also been calculated for conventional radiotherapy fields. The analytical expression for the BECD shows that BECD increases linearly with the variance of the dose. The effect is relatively small, and in the flat regions of the field it results in less than 1% increase of ECD. In the penumbra region of the 6 MV single radiotherapy beam the ECD exceeded the physical dose by up to 35%, when the standard deviation of combined patient setup/organ motion uncertainty was 5 mm. Equivalently, the ECD field was approximately 2 mm wider than the physical dose field. The difference between ECD and the physical dose is greater for normal tissues than for tumours.

A comparison between observed and analytical velocity dispersion profiles of 20 nearby galaxy clusters

NASA Astrophysics Data System (ADS)

Khan, Mohammad S.; Abdullah, Mohamed H.; Ali, Gamal B.

2014-05-01

We derive analytical expression for the velocity dispersion of galaxy clusters, using the statistical mechanical approach. We compare the observed velocity dispersion profiles for 20 nearby ( z≤0.1) galaxy clusters with the analytical ones. It is interesting to find that the analytical results closely match with the observed velocity dispersion profiles only if the presence of the diffuse matter in clusters is taken into consideration. This takes us to introduce a new approach to detect the ratio of diffuse mass, M diff , within a galaxy cluster. For the present sample, the ratio f= M diff / M, where M the cluster's total mass is found to has an average value of 45±12 %. This leads us to the result that nearly 45 % of the cluster mass is impeded outside the galaxies, while around 55 % of the cluster mass is settled in the galaxies.

Corrected Four-Sphere Head Model for EEG Signals.

PubMed

Næss, Solveig; Chintaluri, Chaitanya; Ness, Torbjørn V; Dale, Anders M; Einevoll, Gaute T; Wójcik, Daniel K

2017-01-01

The EEG signal is generated by electrical brain cell activity, often described in terms of current dipoles. By applying EEG forward models we can compute the contribution from such dipoles to the electrical potential recorded by EEG electrodes. Forward models are key both for generating understanding and intuition about the neural origin of EEG signals as well as inverse modeling, i.e., the estimation of the underlying dipole sources from recorded EEG signals. Different models of varying complexity and biological detail are used in the field. One such analytical model is the four-sphere model which assumes a four-layered spherical head where the layers represent brain tissue, cerebrospinal fluid (CSF), skull, and scalp, respectively. While conceptually clear, the mathematical expression for the electric potentials in the four-sphere model is cumbersome, and we observed that the formulas presented in the literature contain errors. Here, we derive and present the correct analytical formulas with a detailed derivation. A useful application of the analytical four-sphere model is that it can serve as ground truth to test the accuracy of numerical schemes such as the Finite Element Method (FEM). We performed FEM simulations of the four-sphere head model and showed that they were consistent with the corrected analytical formulas. For future reference we provide scripts for computing EEG potentials with the four-sphere model, both by means of the correct analytical formulas and numerical FEM simulations.

Universal analytical scattering form factor for shell-, core-shell, or homogeneous particles with continuously variable density profile shape.

PubMed

Foster, Tobias

2011-09-01

A novel analytical and continuous density distribution function with a widely variable shape is reported and used to derive an analytical scattering form factor that allows us to universally describe the scattering from particles with the radial density profile of homogeneous spheres, shells, or core-shell particles. Composed by the sum of two Fermi-Dirac distribution functions, the shape of the density profile can be altered continuously from step-like via Gaussian-like or parabolic to asymptotically hyperbolic by varying a single "shape parameter", d. Using this density profile, the scattering form factor can be calculated numerically. An analytical form factor can be derived using an approximate expression for the original Fermi-Dirac distribution function. This approximation is accurate for sufficiently small rescaled shape parameters, d/R (R being the particle radius), up to values of d/R ≈ 0.1, and thus captures step-like, Gaussian-like, and parabolic as well as asymptotically hyperbolic profile shapes. It is expected that this form factor is particularly useful in a model-dependent analysis of small-angle scattering data since the applied continuous and analytical function for the particle density profile can be compared directly with the density profile extracted from the data by model-free approaches like the generalized inverse Fourier transform method. © 2011 American Chemical Society

Corrected Four-Sphere Head Model for EEG Signals

PubMed Central

Næss, Solveig; Chintaluri, Chaitanya; Ness, Torbjørn V.; Dale, Anders M.; Einevoll, Gaute T.; Wójcik, Daniel K.

2017-01-01

The EEG signal is generated by electrical brain cell activity, often described in terms of current dipoles. By applying EEG forward models we can compute the contribution from such dipoles to the electrical potential recorded by EEG electrodes. Forward models are key both for generating understanding and intuition about the neural origin of EEG signals as well as inverse modeling, i.e., the estimation of the underlying dipole sources from recorded EEG signals. Different models of varying complexity and biological detail are used in the field. One such analytical model is the four-sphere model which assumes a four-layered spherical head where the layers represent brain tissue, cerebrospinal fluid (CSF), skull, and scalp, respectively. While conceptually clear, the mathematical expression for the electric potentials in the four-sphere model is cumbersome, and we observed that the formulas presented in the literature contain errors. Here, we derive and present the correct analytical formulas with a detailed derivation. A useful application of the analytical four-sphere model is that it can serve as ground truth to test the accuracy of numerical schemes such as the Finite Element Method (FEM). We performed FEM simulations of the four-sphere head model and showed that they were consistent with the corrected analytical formulas. For future reference we provide scripts for computing EEG potentials with the four-sphere model, both by means of the correct analytical formulas and numerical FEM simulations. PMID:29093671

Efficient scheme for parametric fitting of data in arbitrary dimensions.

PubMed

Pang, Ning-Ning; Tzeng, Wen-Jer; Kao, Hisen-Ching

2008-07-01

We propose an efficient scheme for parametric fitting expressed in terms of the Legendre polynomials. For continuous systems, our scheme is exact and the derived explicit expression is very helpful for further analytical studies. For discrete systems, our scheme is almost as accurate as the method of singular value decomposition. Through a few numerical examples, we show that our algorithm costs much less CPU time and memory space than the method of singular value decomposition. Thus, our algorithm is very suitable for a large amount of data fitting. In addition, the proposed scheme can also be used to extract the global structure of fluctuating systems. We then derive the exact relation between the correlation function and the detrended variance function of fluctuating systems in arbitrary dimensions and give a general scaling analysis.

Criticality and Connectivity in Macromolecular Charge Complexation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Qin, Jian; de Pablo, Juan J.

We examine the role of molecular connectivity and architecture on the complexation of ionic macromolecules (polyelectrolytes) of finite size. A unified framework is developed and applied to evaluate the electrostatic correlation free energy for point-like, rod-like, and coil-like molecules. That framework is generalized to molecules of variable fractal dimensions, including dendrimers. Analytical expressions for the free energy, correlation length, and osmotic pressure are derived, thereby enabling consideration of the effects of charge connectivity, fractal dimension, and backbone stiffness on the complexation behavior of a wide range of polyelectrolytes. Results are presented for regions in the immediate vicinity of the criticalmore » region and far from it. A transparent and explicit expression for the coexistence curve is derived in order to facilitate analysis of experimentally observed phase diagrams.« less

Analytical Expressions for Deformation from an Arbitrarily Oriented Spheroid in a Half-Space

NASA Astrophysics Data System (ADS)

Cervelli, P. F.

2013-12-01

Deformation from magma chambers can be modeled by an elastic half-space with an embedded cavity subject to uniform pressure change along its interior surface. For a small number of cavity shapes, such as a sphere or a prolate spheroid, closed-form, analytical expressions for deformation have been derived, although these only approximate the uniform-pressure-change boundary condition, with the approximation becoming more accurate as the ratio of source depth to source dimension increases. Using the method of Elshelby [1957] and Yang [1988], which consists of a distribution of double forces and centers of dilatation along the vertical axis, I have derived expressions for displacement from a finite spheroid of arbitrary orientation and aspect ratio that are exact in an infinite elastic medium and approximate in a half-space. The approximation, like those for other cavity shapes, becomes increasingly accurate as the depth to source ratio grows larger, and is accurate to within a few percent in most real-world cases. I have also derived expressions for the deformation-gradient tensor, i.e., the derivatives of each component of displacement with respect to each coordinate direction. These can be transformed easily into the strain and stress tensors. The expressions give deformation both at the surface and at any point within the half-space, and include conditional statements that account for limiting cases that would otherwise prove singular. I have developed MATLAB code for these expressions (and their derivatives), which I use to demonstrate the accuracy of the approximation by showing how well the uniform-pressure-change boundary condition is satisfied in a variety of cases. I also show that a vertical, oblate spheroid with a zero-length vertical axis is equivalent to the penny-shaped crack of Fialko [2001] in an infinite medium and an excellent approximation in a half-space. Finally, because, in many cases, volume change is more tangible than pressure change, I have derived an equation that relates these two quantities for the spheroid: volume change equals pressure change × 2/3 × π/μ × a constant that depends on Poisson's ratio and the spheroid geometry. Eshelby, J. D., The determination of the elastic field of an ellipsoidal inclusion and related problems, Proc. R. Soc. London, Ser. A, 241, 376-396, 1957. Fialko, Y., Khazan, Y, and Simons, M. Deformation due to a pressurized horizontal circular crack in an elastic half-space, with applications to volcano geodesy, Geophys. J. Int., no. 146, 181-190, 2001. Yang, X., Davis, P. M., and Dieterich, J.H, Deformation from inflation of a dipping finite prolate spheroid in an Elastic Half-Space as a model for volcanic stressing, J. Geophys. Res., vol. 93, no. B5, 4249-4257, 1988.

Analytical modeling of the temporal evolution of hot spot temperatures in silicon solar cells

NASA Astrophysics Data System (ADS)

Wasmer, Sven; Rajsrima, Narong; Geisemeyer, Ino; Fertig, Fabian; Greulich, Johannes Michael; Rein, Stefan

2018-03-01

We present an approach to predict the equilibrium temperature of hot spots in crystalline silicon solar cells based on the analysis of their temporal evolution right after turning on a reverse bias. For this end, we derive an analytical expression for the time-dependent heat diffusion of a breakdown channel that is assumed to be cylindrical. We validate this by means of thermography imaging of hot spots right after turning on a reverse bias. The expression allows to be used to extract hot spot powers and radii from short-term measurements, targeting application in inline solar cell characterization. The extracted hot spot powers are validated at the hands of long-term dark lock-in thermography imaging. Using a look-up table of expected equilibrium temperatures determined by numerical and analytical simulations, we utilize the determined hot spot properties to predict the equilibrium temperatures of about 100 industrial aluminum back-surface field solar cells and achieve a high correlation coefficient of 0.86 and a mean absolute error of only 3.3 K.

Morse oscillator propagator in the high temperature limit II: Quantum dynamics and spectroscopy

NASA Astrophysics Data System (ADS)

Toutounji, Mohamad

2018-04-01

This paper is a continuation of Paper I (Toutounji, 2017) of which motivation was testing the applicability of Morse oscillator propagator whose analytical form was derived by Duru (1983). This is because the Morse oscillator propagator was reported (Duru, 1983) in a triple-integral form of a functional of modified Bessel function of the first kind, which considerably limits its applicability. For this reason, I was prompted to find a regime under which Morse oscillator propagator may be simplified and hence be expressed in a closed-form. This was well accomplished in Paper I. Because Morse oscillator is of central importance and widely used in modelling vibrations, its propagator applicability will be extended to applications in quantum dynamics and spectroscopy as will be reported in this paper using the off-diagonal propagator of Morse oscillator whose analytical form is derived.

Joint multifractal analysis based on the partition function approach: analytical analysis, numerical simulation and empirical application

NASA Astrophysics Data System (ADS)

Xie, Wen-Jie; Jiang, Zhi-Qiang; Gu, Gao-Feng; Xiong, Xiong; Zhou, Wei-Xing

2015-10-01

Many complex systems generate multifractal time series which are long-range cross-correlated. Numerous methods have been proposed to characterize the multifractal nature of these long-range cross correlations. However, several important issues about these methods are not well understood and most methods consider only one moment order. We study the joint multifractal analysis based on partition function with two moment orders, which was initially invented to investigate fluid fields, and derive analytically several important properties. We apply the method numerically to binomial measures with multifractal cross correlations and bivariate fractional Brownian motions without multifractal cross correlations. For binomial multifractal measures, the explicit expressions of mass function, singularity strength and multifractal spectrum of the cross correlations are derived, which agree excellently with the numerical results. We also apply the method to stock market indexes and unveil intriguing multifractality in the cross correlations of index volatilities.

Two-population dynamics in a growing network model

NASA Astrophysics Data System (ADS)

Ivanova, Kristinka; Iordanov, Ivan

2012-02-01

We introduce a growing network evolution model with nodal attributes. The model describes the interactions between potentially violent V and non-violent N agents who have different affinities in establishing connections within their own population versus between the populations. The model is able to generate all stable triads observed in real social systems. In the framework of rate equations theory, we employ the mean-field approximation to derive analytical expressions of the degree distribution and the local clustering coefficient for each type of nodes. Analytical derivations agree well with numerical simulation results. The assortativity of the potentially violent network qualitatively resembles the connectivity pattern in terrorist networks that was recently reported. The assortativity of the network driven by aggression shows clearly different behavior than the assortativity of the networks with connections of non-aggressive nature in agreement with recent empirical results of an online social system.

Liquid Drop Model for Charged Spherical Metal Clusters

NASA Astrophysics Data System (ADS)

Seidl, M.; Brack, M.

1996-02-01

The average ground-state energy of a charged spherical metal cluster withNatoms andzexcessive valence electrons, i.e., with net chargeQ=-ezand radiusR=rsN1/3, is presented in the liquid drop model (LDM) expansionE(N, z)=avN+asN2/3+acN1/3+a0(z)+a-1(z) N-1/3+O(N-2/3). We derive analytical expressions for the leading LDM coefficientsav,as,ac, and, in particular, for the charge dependence of the further LDM coefficientsa0anda-1, using the jellium model and density functional theory in the local density approximation. We obtain for the ionization energyI(R)=W+α(e2/R)+O(R-2), with the bulk work functionW=[Φ(+∞)-Φ(0)]-eb, given first by Mahan and Schaich in terms of the electrostatic potentialΦand the bulk energy per electroneb, and a new analytical expression for the dimensionless coefficientα. We demonstrate that within classical theoryα={1}/{2} but, in agreement with experimental information,αtends to ∼0.4 if quantum-mechanical contributions are included. In order to test and confirm our analytical expressions, we discuss the numerical results of semiclassical density variational calculations in the extended Thomas-Fermi model.

Design of permanent magnet eddy current brake for a small scaled electromagnetic launch model

NASA Astrophysics Data System (ADS)

Zhou, Shigui; Yu, Haitao; Hu, Minqiang; Huang, Lei

2012-04-01

A variable pole-pitch double-sided permanent magnet (PM) linear eddy current brake (LECB) is proposed for a small scaled electromagnetic launch model. A two-dimensional (2D) analytical steady state model is presented for the double-sided PM-LECB, and the expression for the braking force is derived. Based on the analytical model, the material and eddy current skin effect of the conducting plate are analyzed. Moreover, a variable pole-pitch double-sided PM-LECB is proposed for the effective braking of the moving plate. In addition, the braking force is predicted by finite element (FE) analysis, and the simulated results are in good agreement with the analytical model. Finally, a prototype is presented to test the braking profile for validation of the proposed design.

Quantitative phase imaging method based on an analytical nonparaxial partially coherent phase optical transfer function.

PubMed

Bao, Yijun; Gaylord, Thomas K

2016-11-01

Multifilter phase imaging with partially coherent light (MFPI-PC) is a promising new quantitative phase imaging method. However, the existing MFPI-PC method is based on the paraxial approximation. In the present work, an analytical nonparaxial partially coherent phase optical transfer function is derived. This enables the MFPI-PC to be extended to the realistic nonparaxial case. Simulations over a wide range of test phase objects as well as experimental measurements on a microlens array verify higher levels of imaging accuracy compared to the paraxial method. Unlike the paraxial version, the nonparaxial MFPI-PC with obliquity factor correction exhibits no systematic error. In addition, due to its analytical expression, the increase in computation time compared to the paraxial version is negligible.

Scaling laws for gas breakdown for nanoscale to microscale gaps at atmospheric pressure

DOE Office of Scientific and Technical Information (OSTI.GOV)

Loveless, Amanda M.; Garner, Allen L., E-mail: algarner@purdue.edu

2016-06-06

Electronics miniaturization motivates gas breakdown predictions for microscale and smaller gaps, since traditional breakdown theory fails when gap size, d, is smaller than ∼15 μm at atmospheric pressure, p{sub atm}. We perform a matched asymptotic analysis to derive analytic expressions for breakdown voltage, V{sub b}, at p{sub atm} for 1 nm ≤ d ≤ 35 μm. We obtain excellent agreement between numerical, analytic, and particle-in-cell simulations for argon, and show V{sub b} decreasing as d → 0, instead of increasing as predicted by Paschen's law. This work provides an analytic framework for determining V{sub b} at atmospheric pressure for various gap distances that may be extended tomore » other gases.« less

Analytical model of radiation-induced precipitation at the surface of dilute binary alloy

NASA Astrophysics Data System (ADS)

Pechenkin, V. A.; Stepanov, I. A.; Konobeev, Yu. V.

2002-12-01

Growth of precipitate layer at the foil surface of an undersaturated binary alloy under uniform irradiation is treated analytically. Analytical expressions for the layer growth rate, layer thickness limit and final component concentrations in the matrix are derived for coherent and incoherent precipitate-matrix interfaces. It is shown that the high temperature limit of radiation-induced precipitation is the same for both types of interfaces, whereas layer thickness limits are different. A parabolic law of the layer growth predicted for both types of interfaces is in agreement with experimental data on γ '-phase precipitation at the surface of Ni-Si dilute alloys under ion irradiation. Effect of sputtering on the precipitation rate and on the low temperature limit of precipitation under ion irradiation is discussed.

Highly Accurate Analytical Approximate Solution to a Nonlinear Pseudo-Oscillator

NASA Astrophysics Data System (ADS)

Wu, Baisheng; Liu, Weijia; Lim, C. W.

2017-07-01

A second-order Newton method is presented to construct analytical approximate solutions to a nonlinear pseudo-oscillator in which the restoring force is inversely proportional to the dependent variable. The nonlinear equation is first expressed in a specific form, and it is then solved in two steps, a predictor and a corrector step. In each step, the harmonic balance method is used in an appropriate manner to obtain a set of linear algebraic equations. With only one simple second-order Newton iteration step, a short, explicit, and highly accurate analytical approximate solution can be derived. The approximate solutions are valid for all amplitudes of the pseudo-oscillator. Furthermore, the method incorporates second-order Taylor expansion in a natural way, and it is of significant faster convergence rate.

Estimating Tsunami Runup with Fault Plane Parameters

NASA Astrophysics Data System (ADS)

Sepulveda, I.; Liu, P. L. F.

2016-12-01

The forecasting of tsunami runup has often been done by solving numerical models. The execution times, however, make them unsuitable for the purpose of warning. We offer an alternative method that provides analytical relationship between the runup height, the fault plane parameters and the characteristic of coastal bathymetry. The method uses the model of Okada (1985) to estimate the coseismic deformation and the corresponding sea surface displacement (η(x,0)). Once the tsunami waves are generated, Carrier & Greenspan (1958) solution (C&G) is adopted to yield analytical expressions for the shoreline elevation and velocity. Two types of problems are investigated. In the first, the bathymetry is modeled as a constant slope that is connected to a constant depth region, where a seismic event occurs. This is a boundary value problem (BVP). In the second, the bathymetry is further simplified as a constant slope, on which a seismic event occurs. This is an initial value problem (IVP). Both problems are depicted in Figure 1. We derive runup solutions in terms of the fault parameters. The earthquake is associated with vertical coseismic seafloor displacements by using Okada's elastic model. In addition to the simplifications considered in Okada's model, we further assume (1) a strike parallel to the shoreline, (2) a very long rupture area and (3) a fast earthquake so surface elevation mimics the seafloor displacements. Then the tsunami origin is modeled in terms of the fault depth (d), fault width (W), fault slip (s) and dip angle (δ). We describe the solution for the BVP. Madsen & Schaeffer (2010) utilized C&G to derive solutions for the shoreline elevation of sinusoidal waves imposed in the offshore boundary. A linear superposition of this solution represents any arbitrary incident wave. Furthermore, we can prescribe the boundary condition at the toe of sloping beach by adopting the linear shallow wave equations in the constant depth area. By means of a dimensional analysis, the runup R is determined by Eq.1. Kanoglu (2004) derived a non-dimensional expression for long wave runup originated over a sloping beach. In our work we determine an analytical expression for a sinusoidal initial condition. Following the same procedure as the BVP, the expression for the runup R in the IVP is given by Eq.2. The curves F1 and F2 are plotted in Figure 2.

Phonon dispersion on Ag (100) surface: A modified analytic embedded atom method study

NASA Astrophysics Data System (ADS)

Xiao-Jun, Zhang; Chang-Le, Chen

2016-01-01

Within the harmonic approximation, the analytic expression of the dynamical matrix is derived based on the modified analytic embedded atom method (MAEAM) and the dynamics theory of surface lattice. The surface phonon dispersions along three major symmetry directions , and X¯M¯ are calculated for the clean Ag (100) surface by using our derived formulas. We then discuss the polarization and localization of surface modes at points X¯ and M¯ by plotting the squared polarization vectors as a function of the layer index. The phonon frequencies of the surface modes calculated by MAEAM are compared with the available experimental and other theoretical data. It is found that the present results are generally in agreement with the referenced experimental or theoretical results, with a maximum deviation of 10.4%. The agreement shows that the modified analytic embedded atom method is a reasonable many-body potential model to quickly describe the surface lattice vibration. It also lays a significant foundation for studying the surface lattice vibration in other metals. Project supported by the National Natural Science Foundation of China (Grant Nos. 61471301 and 61078057), the Scientific Research Program Funded by Shaanxi Provincial Education Department, China (Grant No. 14JK1301), and the Specialized Research Fund for the Doctoral Program of Higher Education, China (Grant No. 20126102110045).

Derivation of a closed form analytical expression for fluorescence recovery after photo bleaching in the case of continuous bleaching during read out

NASA Astrophysics Data System (ADS)

Endress, E.; Weigelt, S.; Reents, G.; Bayerl, T. M.

2005-01-01

Measurements of very slow diffusive processes in membranes, like the diffusion of integral membrane proteins, by fluorescence recovery after photo bleaching (FRAP) are hampered by bleaching of the probe during the read out of the fluorescence recovery. In the limit of long observation time (very slow diffusion as in the case of large membrane proteins), this bleaching may cause errors to the recovery function and thus provides error-prone diffusion coefficients. In this work we present a new approach to a two-dimensional closed form analytical solution of the reaction-diffusion equation, based on the addition of a dissipative term to the conventional diffusion equation. The calculation was done assuming (i) a Gaussian laser beam profile for bleaching the spot and (ii) that the fluorescence intensity profile emerging from the spot can be approximated by a two-dimensional Gaussian. The detection scheme derived from the analytical solution allows for diffusion measurements without the constraint of observation bleaching. Recovery curves of experimental FRAP data obtained under non-negligible read-out bleaching for native membranes (rabbit endoplasmic reticulum) on a planar solid support showed excellent agreement with the analytical solution and allowed the calculation of the lipid diffusion coefficient.

Phase transformations in an ascending adiabatic mixed-phase cloud volume

NASA Astrophysics Data System (ADS)

Pinsky, M.; Khain, A.; Korolev, A.

2015-04-01

Regimes of liquid-ice coexistence that may form in an adiabatic parcel ascending at constant velocity at freezing temperatures are investigated. Four zones with different microphysical structures succeeding one another along the vertical direction have been established. On the basis of a novel balance equation, analytical expressions are derived to determine the conditions specific for each of these zones. In particular, the necessary and sufficient conditions for formation of liquid water phase within an ascending parcel containing only ice particles are determined. The results are compared to findings reported in earlier studies. The role of the Wegener-Bergeron-Findeisen mechanism in the phase transformation is analyzed. The dependence of the phase relaxation time on height in the four zones is investigated on the basis of a novel analytical expression. The results obtained in the study can be instrumental for analysis and interpretation of observed mixed-phase clouds.

QCD corrections to decay-lepton polar and azimuthal angular distributions in e+e- --> tt(bar) in the soft-gluon approximation

NASA Astrophysics Data System (ADS)

Rindani, Saurabh D.

2002-04-01

QCD corrections to order as in the soft-gluon approximation to angular distributions of decay charged leptons in the process e+e- --> t t(bar), followed by semileptonic decay of t or t(bar), are obtained in the e+e- centre-of-mass frame. As compared to distributions in the top rest frame, these have the advantage that they would allow direct comparison with experiment without the need to reconstruct the top rest frame. The results also do not depend on the choice of a spin quantization axis for t or t (bar). Analytic expression for the triple distribution in the polar angle of t and polar and azimuthal angles of the lepton is obtained. Analytic expression is also derived for the distribution in the charged-lepton polar angle. Numerical values are discussed for (s) 1/2 = 400, 800 and 1500 GeV.

Band-to-band tunneling distance analysis in the heterogate electron–hole bilayer tunnel field-effect transistor

DOE Office of Scientific and Technical Information (OSTI.GOV)

Padilla, J. L., E-mail: jose.padilladelatorre@epfl.ch; Departamento de Electrónica y Tecnología de los Computadores, Universidad de Granada, Avda. Fuentenueva s/n, 18071 Granada; Palomares, A.

In this work, we analyze the behavior of the band-to-band tunneling distance between electron and hole subbands resulting from field-induced quantum confinement in the heterogate electron–hole bilayer tunnel field-effect transistor. We show that, analogously to the explicit formula for the tunneling distance that can be easily obtained in the semiclassical framework where the conduction and valence band edges are allowed states, an equivalent analytical expression can be derived in the presence of field-induced quantum confinement for describing the dependence of the tunneling distance on the body thickness and material properties of the channel. This explicit expression accounting for quantum confinementmore » holds valid provided that the potential wells for electrons and holes at the top and bottom of the channel can be approximated by triangular profiles. Analytical predictions are compared to simulation results showing very accurate agreement.« less

Bayesian estimation of the discrete coefficient of determination.

PubMed

Chen, Ting; Braga-Neto, Ulisses M

2016-12-01

The discrete coefficient of determination (CoD) measures the nonlinear interaction between discrete predictor and target variables and has had far-reaching applications in Genomic Signal Processing. Previous work has addressed the inference of the discrete CoD using classical parametric and nonparametric approaches. In this paper, we introduce a Bayesian framework for the inference of the discrete CoD. We derive analytically the optimal minimum mean-square error (MMSE) CoD estimator, as well as a CoD estimator based on the Optimal Bayesian Predictor (OBP). For the latter estimator, exact expressions for its bias, variance, and root-mean-square (RMS) are given. The accuracy of both Bayesian CoD estimators with non-informative and informative priors, under fixed or random parameters, is studied via analytical and numerical approaches. We also demonstrate the application of the proposed Bayesian approach in the inference of gene regulatory networks, using gene-expression data from a previously published study on metastatic melanoma.

Critical power for self-focusing of optical beam in absorbing media

NASA Astrophysics Data System (ADS)

Qi, Pengfei; Zhang, Lin; Lin, Lie; Zhang, Nan; Wang, Yan; Liu, Weiwei

2018-04-01

Self-focusing effects are of central importance for most nonlinear optical effects. The critical power for self-focusing is commonly investigated theoretically without considering a material’s absorption. Although this is practicable for various materials, investigating the critical power for self-focusing in media with non-negligible absorption is also necessary, because this is the situation usually met in practice. In this paper, the simple analytical expressions describing the relationships among incident power, absorption coefficient and focal position are provided by a simple physical model based on the Fermat principle. Expressions for the absorption dependent critical power are also derived; these can play important roles in experimental and applied research on self-focusing-related nonlinear optical phenomena in absorbing media. Numerical results, based on the nonlinear wave equation—and which can predict experimental results perfectly—are also presented, and agree quantitatively with the analytical results proposed in this paper.

Interaction of a conductive crack and of an electrode at a piezoelectric bimaterial interface

NASA Astrophysics Data System (ADS)

Onopriienko, Oleg; Loboda, Volodymyr; Sheveleva, Alla; Lapusta, Yuri

2018-06-01

The interaction of a conductive crack and an electrode at a piezoelectric bi-material interface is studied. The bimaterial is subjected to an in-plane electrical field parallel to the interface and an anti-plane mechanical loading. The problem is formulated and reduced, via the application of sectionally analytic vector functions, to a combined Dirichlet-Riemann boundary value problem. Simple analytical expressions for the stress, the electric field, and their intensity factors as well as for the crack faces' displacement jump are derived. Our numerical results illustrate the proposed approach and permit to draw some conclusions on the crack-electrode interaction.

Analytical model for electromagnetic radiation from a wakefield excited by intense short laser pulses in an unmagnetized plasma

NASA Astrophysics Data System (ADS)

Chen, Zi-Yu; Chen, Shi; Dan, Jia-Kun; Li, Jian-Feng; Peng, Qi-Xian

2011-10-01

A simple one-dimensional analytical model for electromagnetic emission from an unmagnetized wakefield excited by an intense short-pulse laser in the nonlinear regime has been developed in this paper. The expressions for the spectral and angular distributions of the radiation have been derived. The model suggests that the origin of the radiation can be attributed to the violent sudden acceleration of plasma electrons experiencing the accelerating potential of the laser wakefield. The radiation process could help to provide a qualitative interpretation of existing experimental results, and offers useful information for future laser wakefield experiments.

A closed form of a kurtosis parameter of a hypergeometric-Gaussian type-II beam

NASA Astrophysics Data System (ADS)

F, Khannous; A, A. A. Ebrahim; A, Belafhal

2016-04-01

Based on the irradiance moment definition and the analytical expression of waveform propagation for hypergeometric-Gaussian type-II beams passing through an ABCD system, the kurtosis parameter is derived analytically and illustrated numerically. The kurtosis parameters of the Gaussian beam, modified Bessel modulated Gaussian beam with quadrature radial and elegant Laguerre-Gaussian beams are obtained by treating them as special cases of the present treatment. The obtained results show that the kurtosis parameter depends on the change of the beam order m and the hollowness parameter p, such as its decrease with increasing m and increase with increasing p.

Limits on estimating the width of thin tubular structures in 3D images.

PubMed

Wörz, Stefan; Rohr, Karl

2006-01-01

This work studies limits on estimating the width of thin tubular structures in 3D images. Based on nonlinear estimation theory we analyze the minimal stochastic error of estimating the width. Given a 3D analytic model of the image intensities of tubular structures, we derive a closed-form expression for the Cramér-Rao bound of the width estimate under image noise. We use the derived lower bound as a benchmark and compare it with three previously proposed accuracy limits for vessel width estimation. Moreover, by experimental investigations we demonstrate that the derived lower bound can be achieved by fitting a 3D parametric intensity model directly to the image data.

Laugh and Smile upon the Holy Quran: The Study of Analytical Objectivities

ERIC Educational Resources Information Center

al-Domi, Mohammad Mahmoud

2015-01-01

This study aims to examine the positive impact of The Holy Quran based on the laugh and smile. This kind of derivatives in which context of praise, expression the feeling of happiness and joyful in the positive senses. Everyone needs to relieve his heart so that happiness and joy on their faces can be seen. Laughter also are some of attribute…

Electron-cyclotron absorption in high-temperature plasmas: quasi-exact analytical evaluation and comparative numerical analysis

NASA Astrophysics Data System (ADS)

Albajar, F.; Bertelli, N.; Bornatici, M.; Engelmann, F.

2007-01-01

On the basis of the electromagnetic energy balance equation, a quasi-exact analytical evaluation of the electron-cyclotron (EC) absorption coefficient is performed for arbitrary propagation (with respect to the magnetic field) in a (Maxwellian) magneto-plasma for the temperature range of interest for fusion reactors (in which EC radiation losses tend to be important in the plasma power balance). The calculation makes use of Bateman's expansion for the product of two Bessel functions, retaining the lowest-order contribution. The integration over electron momentum can then be carried out analytically, fully accounting for finite Larmor radius effects in this approximation. On the basis of the analytical expressions for the EC absorption coefficients of both the extraordinary and ordinary modes thus obtained, (i) for the case of perpendicular propagation simple formulae are derived for both modes and (ii) a numerical analysis of the angular distribution of EC absorption is carried out. An assessment of the accuracy of asymptotic expressions that have been given earlier is also performed, showing that these approximations can be usefully applied for calculating EC power losses from reactor-grade plasmas. Presented in part at the 14th Joint Workshop on Electron Cyclotron Emission and Electron Cyclotron Resonance Heating, Santorini, Greece, 9-12 May 2006.

Analytic second derivative of the energy for density functional theory based on the three-body fragment molecular orbital method

NASA Astrophysics Data System (ADS)

Nakata, Hiroya; Fedorov, Dmitri G.; Zahariev, Federico; Schmidt, Michael W.; Kitaura, Kazuo; Gordon, Mark S.; Nakamura, Shinichiro

2015-03-01

Analytic second derivatives of the energy with respect to nuclear coordinates have been developed for spin restricted density functional theory (DFT) based on the fragment molecular orbital method (FMO). The derivations were carried out for the three-body expansion (FMO3), and the two-body expressions can be obtained by neglecting the three-body corrections. Also, the restricted Hartree-Fock (RHF) Hessian for FMO3 can be obtained by neglecting the density-functional related terms. In both the FMO-RHF and FMO-DFT Hessians, certain terms with small magnitudes are neglected for computational efficiency. The accuracy of the FMO-DFT Hessian in terms of the Gibbs free energy is evaluated for a set of polypeptides and water clusters and found to be within 1 kcal/mol of the corresponding full (non-fragmented) ab initio calculation. The FMO-DFT method is also applied to transition states in SN2 reactions and for the computation of the IR and Raman spectra of a small Trp-cage protein (PDB: 1L2Y). Some computational timing analysis is also presented.

Constant-roll tachyon inflation and observational constraints

NASA Astrophysics Data System (ADS)

Gao, Qing; Gong, Yungui; Fei, Qin

2018-05-01

For the constant-roll tachyon inflation, we derive the analytical expressions for the scalar and tensor power spectra, the scalar and tensor spectral tilts and the tensor to scalar ratio to the first order of epsilon1 by using the method of Bessel function approximation. The derived ns-r results are compared with the observations, we find that only the constant-roll inflation with ηH being a constant is consistent with the observations and observations constrain the constant-roll inflation to be slow-roll inflation. The tachyon potential is also reconstructed for the constant-roll inflation which is consistent with the observations.

Ion transfer through solvent polymeric membranes driven by an exponential current flux.

PubMed

Molina, A; Torralba, E; González, J; Serna, C; Ortuño, J A

2011-03-21

General analytical equations which govern ion transfer through liquid membranes with one and two polarized interfaces driven by an exponential current flux are derived. Expressions for the transient and stationary E-t, dt/dE-E and dI/dE-E curves are obtained, and the evolution from transient to steady behaviour has been analyzed in depth. We have also shown mathematically that the voltammetric and stationary chronopotentiometric I(N)-E curves are identical (with E being the applied potential for voltammetric techniques and the measured potential for chronopotentiometric techniques), and hence, their derivatives provide identical information.

Approximate method of variational Bayesian matrix factorization/completion with sparse prior

NASA Astrophysics Data System (ADS)

Kawasumi, Ryota; Takeda, Koujin

2018-05-01

We derive the analytical expression of a matrix factorization/completion solution by the variational Bayes method, under the assumption that the observed matrix is originally the product of low-rank, dense and sparse matrices with additive noise. We assume the prior of a sparse matrix is a Laplace distribution by taking matrix sparsity into consideration. Then we use several approximations for the derivation of a matrix factorization/completion solution. By our solution, we also numerically evaluate the performance of a sparse matrix reconstruction in matrix factorization, and completion of a missing matrix element in matrix completion.

Development of ultrasound transducer diffractive field theory for nonlinear propagation-based imaging

NASA Astrophysics Data System (ADS)

Kharin, Nikolay A.

2000-04-01

In nonlinear ultrasound imaging the images are formed using the second harmonic energy generated due to the nonlinear nature of finite amplitude propagation. This propagation can be modeled using the KZK wave equation. This paper presents further development of nonlinear diffractive field theory based on the KZK equation and its solution by means of the slowly changing profile method for moderate nonlinearity. The analytical expression for amplitudes and phases of sum frequency wave are obtained in addition to the second harmonic wave. Also, the analytical expression for the relative curvature of the wave fronts of fundamental and second harmonic signals are derived. The media with different nonlinear properties and absorption coefficients were investigated to characterize the diffractive field of the transducer at medical frequencies. All expressions demonstrate good agreement with experimental results. The expressions are novel and provide an easy way for prediction of amplitude and phase structure of nonlinearly distorted field of a transducer. The sum frequency signal technique could be implemented as well as second harmonic technique to improve the quality of biomedical images. The results obtained are of importance for medical diagnostic ultrasound equipment design.

On the connection coefficients and recurrence relations arising from expansions in series of Laguerre polynomials

NASA Astrophysics Data System (ADS)

Doha, E. H.

2003-05-01

A formula expressing the Laguerre coefficients of a general-order derivative of an infinitely differentiable function in terms of its original coefficients is proved, and a formula expressing explicitly the derivatives of Laguerre polynomials of any degree and for any order as a linear combination of suitable Laguerre polynomials is deduced. A formula for the Laguerre coefficients of the moments of one single Laguerre polynomial of certain degree is given. Formulae for the Laguerre coefficients of the moments of a general-order derivative of an infinitely differentiable function in terms of its Laguerre coefficients are also obtained. A simple approach in order to build and solve recursively for the connection coefficients between Jacobi-Laguerre and Hermite-Laguerre polynomials is described. An explicit formula for these coefficients between Jacobi and Laguerre polynomials is given, of which the ultra-spherical polynomials of the first and second kinds and Legendre polynomials are important special cases. An analytical formula for the connection coefficients between Hermite and Laguerre polynomials is also obtained.

Identification of FOPDT and SOPDT process dynamics using closed loop test.

PubMed

Bajarangbali, Raghunath; Majhi, Somanath; Pandey, Saurabh

2014-07-01

In this paper, identification of stable and unstable first order, second order overdamped and underdamped process dynamics with time delay is presented. Relay with hysteresis is used to induce a limit cycle output and using this information, unknown process model parameters are estimated. State space based generalized analytical expressions are derived to achieve accurate results. To show the performance of the proposed method expressions are also derived for systems with a zero. In real time systems, measurement noise is an important issue during identification of process dynamics. A relay with hysteresis reduces the effect of measurement noise, in addition a new multiloop control strategy is proposed to recover the original limit cycle. Simulation results are included to validate the effectiveness of the proposed method. Copyright © 2014 ISA. Published by Elsevier Ltd. All rights reserved.

On the changes of IAU 2000 nutation theory stemming from IAU 2006 precession theory

NASA Astrophysics Data System (ADS)

Escapa, A.; Getino, J.; Ferrándiz, J. M.; Baenas, T.

2014-12-01

The adoption of IAU 2006 precession theory (Capitaine et al. 2003) introduced some small changes in IAU 2000A nutation theory, relevant at the mircroarcsecond level. These adjustments were derived in Capitaine et al. (2005) and are currently considered in international standards like, for example, IERS Conventions (2010) or in the Explanatory Supplement to the Astronomical Almanac (2013). We reexamine the issue, working out the induced modifications due to a change in the value of the obliquity of the ecliptic and to the secular variation of the Earth dynamical flattening. In particular, within the framework of the Hamiltonian theory of the rotation of the Earth we derive analytical expressions of those changes for the motion of the figure axis. These expressions and their corresponding numerical contributions will be compared with those obtained in Capitaine et al. (2005).

Turbulent equipartition pinch of toroidal momentum in spherical torus

NASA Astrophysics Data System (ADS)

Hahm, T. S.; Lee, J.; Wang, W. X.; Diamond, P. H.; Choi, G. J.; Na, D. H.; Na, Y. S.; Chung, K. J.; Hwang, Y. S.

2014-12-01

We present a new analytic expression for turbulent equipartition (TEP) pinch of toroidal angular momentum originating from magnetic field inhomogeneity of spherical torus (ST) plasmas. Starting from a conservative modern nonlinear gyrokinetic equation (Hahm et al 1988 Phys. Fluids 31 2670), we derive an expression for pinch to momentum diffusivity ratio without using a usual tokamak approximation of B ∝ 1/R which has been previously employed for TEP momentum pinch derivation in tokamaks (Hahm et al 2007 Phys. Plasmas 14 072302). Our new formula is evaluated for model equilibria of National Spherical Torus eXperiment (NSTX) (Ono et al 2001 Nucl. Fusion 41 1435) and Versatile Experiment Spherical Torus (VEST) (Chung et al 2013 Plasma Sci. Technol. 15 244) plasmas. Our result predicts stronger inward pinch for both cases, as compared to the prediction based on the tokamak formula.

Generalisation of Gilbert damping and magnetic inertia parameter as a series of higher-order relativistic terms

NASA Astrophysics Data System (ADS)

Mondal, Ritwik; Berritta, Marco; Oppeneer, Peter M.

2018-07-01

The phenomenological Landau–Lifshitz–Gilbert (LLG) equation of motion remains as the cornerstone of contemporary magnetisation dynamics studies, wherein the Gilbert damping parameter has been attributed to first-order relativistic effects. To include magnetic inertial effects the LLG equation has previously been extended with a supplemental inertia term; the arising inertial dynamics has been related to second-order relativistic effects. Here we start from the relativistic Dirac equation and, performing a Foldy–Wouthuysen transformation, derive a generalised Pauli spin Hamiltonian that contains relativistic correction terms to any higher order. Using the Heisenberg equation of spin motion we derive general relativistic expressions for the tensorial Gilbert damping and magnetic inertia parameters, and show that these tensors can be expressed as series of higher-order relativistic correction terms. We further show that, in the case of a harmonic external driving field, these series can be summed and we provide closed analytical expressions for the Gilbert and inertial parameters that are functions of the frequency of the driving field.

Generalisation of Gilbert damping and magnetic inertia parameter as a series of higher-order relativistic terms.

PubMed

Mondal, Ritwik; Berritta, Marco; Oppeneer, Peter M

2018-05-17

The phenomenological Landau-Lifshitz-Gilbert (LLG) equation of motion remains as the cornerstone of contemporary magnetisation dynamics studies, wherein the Gilbert damping parameter has been attributed to first-order relativistic effects. To include magnetic inertial effects the LLG equation has previously been extended with a supplemental inertia term; the arising inertial dynamics has been related to second-order relativistic effects. Here we start from the relativistic Dirac equation and, performing a Foldy-Wouthuysen transformation, derive a generalised Pauli spin Hamiltonian that contains relativistic correction terms to any higher order. Using the Heisenberg equation of spin motion we derive general relativistic expressions for the tensorial Gilbert damping and magnetic inertia parameters, and show that these tensors can be expressed as series of higher-order relativistic correction terms. We further show that, in the case of a harmonic external driving field, these series can be summed and we provide closed analytical expressions for the Gilbert and inertial parameters that are functions of the frequency of the driving field.

Target-based coherent beam combining of an optical phased array fed by a broadband laser source

NASA Astrophysics Data System (ADS)

Hyde, Milo W., IV; McCrae, Jack E.; Tyler, Glenn A.

2017-11-01

The target-based phasing of an optical phased array (OPA) fed by a broadband master oscillator laser source is investigated. The specific scenario examined here considers an OPA phasing through atmospheric turbulence on a rough curved object. An analytical expression for the detected or received intensity is derived. Gleaned from this expression are the conditions under which target-based phasing is possible. A detailed OPA wave optics simulation is performed to validate the theoretical findings. Key aspects of the simulation set-up as well as the results are thoroughly discussed.

A Novel Capacity Analysis for Wireless Backhaul Mesh Networks

NASA Astrophysics Data System (ADS)

Chung, Tein-Yaw; Lee, Kuan-Chun; Lee, Hsiao-Chih

This paper derived a closed-form expression for inter-flow capacity of a backhaul wireless mesh network (WMN) with centralized scheduling by employing a ring-based approach. Through the definition of an interference area, we are able to accurately describe a bottleneck collision area for a WMN and calculate the upper bound of inter-flow capacity. The closed-form expression shows that the upper bound is a function of the ratio between transmission range and network radius. Simulations and numerical analysis show that our analytic solution can better estimate the inter-flow capacity of WMNs than that of previous approach.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Williams, S.; Rohlfing, E.A.; Rahn, L.A.

We present analytical signal expressions for each two-color resonant four-wave mixing (TC-RFWM) scheme that can be used for double-resonance molecular spectroscopy in the limit of weak fields (no saturation). The theoretical approach employs time-independent, diagrammatic perturbation theory and a spherical tensor analysis in an extension of recent treatments of degenerate four-wave mixing [S. Williams, R. N. Zare, and L. A. Rahn, J. Chem. Phys. {bold 101}, 1072 (1994)] and TC-RFWM for the specific case of stimulated emission pumping [S. Williams {ital et al.}, J. Chem. Phys. {bold 102}, 8342 (1995)]. Under the assumption that the relaxation of the population, themore » orientation, and the alignment are the same, simple analytic expressions are derived for commonly used experimental configurations. The TC-RFWM signal is found to be a product of a concentration term, a one-photon molecular term, a line shape function, and a laboratory-frame geometric factor. These expressions are intended to facilitate the practical analysis of TC-RFWM spectra by clarifying, for example, the dependence on beam polarizations and rotational branch combinations. {copyright} {ital 1997 American Institute of Physics.}« less

Quasi-analytical treatment of spatially averaged radiation transfer in complex terrain

NASA Astrophysics Data System (ADS)

Löwe, H.; Helbig, N.

2012-04-01

We provide a new quasi-analytical method to compute the topographic influence on the effective albedo of complex topography as required for meteorological, land-surface or climate models. We investigate radiative transfer in complex terrain via the radiosity equation on isotropic Gaussian random fields. Under controlled approximations we derive expressions for domain averages of direct, diffuse and terrain radiation and the sky view factor. Domain averaged quantities are related to a type of level-crossing probability of the random field which is approximated by longstanding results developed for acoustic scattering at ocean boundaries. This allows us to express all non-local horizon effects in terms of a local terrain parameter, namely the mean squared slope. Emerging integrals are computed numerically and fit formulas are given for practical purposes. As an implication of our approach we provide an expression for the effective albedo of complex terrain in terms of the sun elevation angle, mean squared slope, the area averaged surface albedo, and the direct-to-diffuse ratio of solar radiation. As an application, we compute the effective albedo for the Swiss Alps and discuss possible generalizations of the method.

Comprehensive Reactive Receiver Modeling for Diffusive Molecular Communication Systems: Reversible Binding, Molecule Degradation, and Finite Number of Receptors.

PubMed

Ahmadzadeh, Arman; Arjmandi, Hamidreza; Burkovski, Andreas; Schober, Robert

2016-10-01

This paper studies the problem of receiver modeling in molecular communication systems. We consider the diffusive molecular communication channel between a transmitter nano-machine and a receiver nano-machine in a fluid environment. The information molecules released by the transmitter nano-machine into the environment can degrade in the channel via a first-order degradation reaction and those that reach the receiver nano-machine can participate in a reversible bimolecular reaction with receiver receptor proteins. Thereby, we distinguish between two scenarios. In the first scenario, we assume that the entire surface of the receiver is covered by receptor molecules. We derive a closed-form analytical expression for the expected received signal at the receiver, i.e., the expected number of activated receptors on the surface of the receiver. Then, in the second scenario, we consider the case where the number of receptor molecules is finite and the uniformly distributed receptor molecules cover the receiver surface only partially. We show that the expected received signal for this scenario can be accurately approximated by the expected received signal for the first scenario after appropriately modifying the forward reaction rate constant. The accuracy of the derived analytical results is verified by Brownian motion particle-based simulations of the considered environment, where we also show the impact of the effect of receptor occupancy on the derived analytical results.

Symmetry analysis of strain, electric and magnetic fields in the Bi2Se3-class of topological insulators

NASA Astrophysics Data System (ADS)

Rosdahl Brems, Mathias; Paaske, Jens; Lunde, Anders Mathias; Willatzen, Morten

2018-05-01

Based on group theoretical arguments we derive the most general Hamiltonian for the Bi2Se3-class of materials including terms to third order in the wave vector, first order in electric and magnetic fields, first order in strain and first order in both strain and wave vector. We determine analytically the effects of strain on the electronic structure of Bi2Se3. For the most experimentally relevant surface termination we analytically derive the surface state (SS) spectrum, revealing an anisotropic Dirac cone with elliptical constant energy contours giving rise to a direction-dependent group velocity. The spin-momentum locking of strained Bi2Se3 is shown to be modified. Hence, strain control can be used to manipulate the spin degree of freedom via the spin–orbit coupling. We show that for a thin film of Bi2Se3 the SS band gap induced by coupling between the opposite surfaces changes opposite to the bulk band gap under strain. Tuning the SS band gap by strain, gives new possibilities for the experimental investigation of the thickness dependent gap and optimization of optical properties relevant for, e.g., photodetector and energy harvesting applications. We finally derive analytical expressions for the effective mass tensor of the Bi2Se3 class of materials as a function of strain and electric field.

Nonlinear pulse propagation and phase velocity of laser-driven plasma waves

NASA Astrophysics Data System (ADS)

Benedetti, Carlo; Rossi, Francesco; Schroeder, Carl; Esarey, Eric; Leemans, Wim

2014-10-01

We investigate and characterize the laser evolution and plasma wave excitation by a relativistically intense, short-pulse laser propagating in a preformed parabolic plasma channel, including the effects of pulse steepening, frequency redshifting, and energy depletion. We derived in 3D, and in the weakly relativistic intensity regime, analytical expressions for the laser energy depletion, the pulse self-steepening rate, the laser intensity centroid velocity, and the phase velocity of the plasma wave. Analytical results have been validated numerically using the 2D-cylindrical, ponderomotive code INF&RNO. We also discuss the extension of these results to the nonlinear regime, where an analytical theory of the nonlinear wake phase velocity is lacking. Work supported by the Office of Science, Office of High Energy Physics, of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.

Photon statistics of a two-mode squeezed vacuum

NASA Technical Reports Server (NTRS)

Schrade, Guenter; Akulin, V. M.; Schleich, W. P.; Manko, Vladimir I.

1994-01-01

We investigate the general case of the photon distribution of a two-mode squeezed vacuum and show that the distribution of photons among the two modes depends on four parameters: two squeezing parameters, the relative phase between the two oscillators and their spatial orientation. The distribution of the total number of photons depends only on the two squeezing parameters. We derive analytical expressions and present pictures for both distributions.

INTERNATIONAL CONFERENCE ON SEMICONDUCTOR INJECTION LASERS SELCO-87: Simple formula for the thermal conductivity of a quaternary solid solution

NASA Astrophysics Data System (ADS)

Nakwaski, W.

1988-11-01

An analysis is made of the thermal conductivity of quaternary solid solutions (alloys) allowing for their disordered structure on the basis of a phenomenological analysis proposed by Abeles. This method is applied to a quaternary solid solution In1 - xGaxAsyP1 - y. A simple analytic expression is derived for the thermal conductivity of this material.

Closed-form recursive formula for an optimal tracker with terminal constraints

NASA Technical Reports Server (NTRS)

Juang, J. N.; Turner, J. D.; Chun, H. M.

1986-01-01

Feedback control laws are derived for a class of optimal finite time tracking problems with terminal constraints. Analytical solutions are obtained for the feedback gain and the closed-loop response trajectory. Such formulations are expressed in recursive forms so that a real-time computer implementation becomes feasible. An example involving the feedback slewing of a flexible spacecraft is given to illustrate the validity and usefulness of the formulations.

Explicit formulas for effective piezoelectric coefficients of ferroelectric 0-3 composites based on effective medium theory

NASA Astrophysics Data System (ADS)

Wong, C. K.; Poon, Y. M.; Shin, F. G.

2003-01-01

Explicit formulas were derived for the effective piezoelectric stress coefficients of a 0-3 composite of ferroelectric spherical particles in a ferroelectric matrix which were then combined to give the more commonly used strain coefficients. Assuming that the elastic stiffness of the inclusion phase is sufficiently larger than that of the matrix phase, the previously derived explicit expressions for the case of a low volume concentration of inclusion particles [C. K. Wong, Y. M. Poon, and F. G. Shin, Ferroelectrics 264, 39 (2001); J. Appl. Phys. 90, 4690 (2001)] were "transformed" analytically by an effective medium theory (EMT) with appropriate approximations, to suit the case of a more concentrated suspension. Predictions of the EMT expressions were compared with the experimental values of composites of lead zirconate titanate ceramic particles dispersed in polyvinylidene fluoride and polyvinylidene fluoride-trifluoroethylene copolymer, reported by Furukawa [IEEE Trans. Electr. Insul. 24, 375 (1989)] and by Ng et al. [IEEE Trans. Ultrason. Ferroelectr. Freq. Control 47, 1308 (2000)] respectively. Fairly good agreement was obtained. Comparisons with other predictions, including the predictions given by numerically solving the EMT scheme, were also made. It was found that the analytic and numeric EMT schemes agreed with each other very well for an inclusion of volume fraction not exceeding 60%.

Field-driven chiral bubble dynamics analysed by a semi-analytical approach

NASA Astrophysics Data System (ADS)

Vandermeulen, J.; Leliaert, J.; Dupré, L.; Van Waeyenberge, B.

2017-12-01

Nowadays, field-driven chiral bubble dynamics in the presence of the Dzyaloshinskii-Moriya interaction are a topic of thorough investigation. In this paper, a semi-analytical approach is used to derive equations of motion that express the bubble wall (BW) velocity and the change in in-plane magnetization angle as function of the micromagnetic parameters of the involved interactions, thereby taking into account the two-dimensional nature of the bubble wall. It is demonstrated that the equations of motion enable an accurate description of the expanding and shrinking convex bubble dynamics and an expression for the transition field between shrinkage and expansion is derived. In addition, these equations of motion show that the BW velocity is not only dependent on the driving force, but also on the BW curvature. The absolute BW velocity increases for both a shrinking and an expanding bubble, but for different reasons: for expanding bubbles, it is due to the increasing importance of the driving force, while for shrinking bubbles, it is due to the increasing importance of contributions related to the BW curvature. Finally, using this approach we show how the recently proposed magnetic bubblecade memory can operate in the flow regime in the presence of a tilted sinusoidal magnetic field and at greatly reduced bubble sizes compared to the original device prototype.

The Analytical Solution of the Transient Radial Diffusion Equation with a Nonuniform Loss Term.

NASA Astrophysics Data System (ADS)

Loridan, V.; Ripoll, J. F.; De Vuyst, F.

2017-12-01

Many works have been done during the past 40 years to perform the analytical solution of the radial diffusion equation that models the transport and loss of electrons in the magnetosphere, considering a diffusion coefficient proportional to a power law in shell and a constant loss term. Here, we propose an original analytical method to address this challenge with a nonuniform loss term. The strategy is to match any L-dependent electron losses with a piecewise constant function on M subintervals, i.e., dealing with a constant lifetime on each subinterval. Applying an eigenfunction expansion method, the eigenvalue problem becomes presently a Sturm-Liouville problem with M interfaces. Assuming the continuity of both the distribution function and its first spatial derivatives, we are able to deal with a well-posed problem and to find the full analytical solution. We further show an excellent agreement between both the analytical solutions and the solutions obtained directly from numerical simulations for different loss terms of various shapes and with a diffusion coefficient DLL L6. We also give two expressions for the required number of eigenmodes N to get an accurate snapshot of the analytical solution, highlighting that N is proportional to 1/√t0, where t0 is a time of interest, and that N increases with the diffusion power. Finally, the equilibrium time, defined as the time to nearly reach the steady solution, is estimated by a closed-form expression and discussed. Applications to Earth and also Jupiter and Saturn are discussed.

Generation of dark hollow beams by using a fractional radial Hilbert transform system

NASA Astrophysics Data System (ADS)

Xie, Qiansen; Zhao, Daomu

2007-07-01

The radial Hilbert transform has been extend to the fractional field, which could be called the fractional radial Hilbert transform (FRHT). Using edge-enhancement characteristics of this transform, we convert a Gaussian light beam into a variety of dark hollow beams (DHBs). Based on the fact that a hard-edged aperture can be expanded approximately as a finite sum of complex Gaussian functions, the analytical expression of a Gaussian beam passing through a FRHT system has been derived. As a numerical example, the properties of the DHBs with different fractional orders are illustrated graphically. The calculation results obtained by use of the analytical method and the integral method are also compared.

Exact Local Correlations and Full Counting Statistics for Arbitrary States of the One-Dimensional Interacting Bose Gas

NASA Astrophysics Data System (ADS)

Bastianello, Alvise; Piroli, Lorenzo; Calabrese, Pasquale

2018-05-01

We derive exact analytic expressions for the n -body local correlations in the one-dimensional Bose gas with contact repulsive interactions (Lieb-Liniger model) in the thermodynamic limit. Our results are valid for arbitrary states of the model, including ground and thermal states, stationary states after a quantum quench, and nonequilibrium steady states arising in transport settings. Calculations for these states are explicitly presented and physical consequences are critically discussed. We also show that the n -body local correlations are directly related to the full counting statistics for the particle-number fluctuations in a short interval, for which we provide an explicit analytic result.

Estimation of the curvature of the solid liquid interface during Bridgman crystal growth

NASA Astrophysics Data System (ADS)

Barat, Catherine; Duffar, Thierry; Garandet, Jean-Paul

1998-11-01

An approximate solution for the solid/liquid interface curvature due to the crucible effect in crystal growth is derived from simple heat flux considerations. The numerical modelling of the problem carried out with the help of the finite element code FIDAP supports the predictions of our analytical expression and allows to identify its range of validity. Experimental interface curvatures, measured in gallium antimonide samples grown by the vertical Bridgman method, are seen to compare satisfactorily to analytical and numerical results. Other literature data are also in fair agreement with the predictions of our models in the case where the amount of heat carried by the crucible is small compared to the overall heat flux.

Optical absorption in planar graphene superlattice: The role of structural parameters

NASA Astrophysics Data System (ADS)

Azadi, L.; Shojaei, S.

2018-04-01

We theoretically studied the optically driven interband transitions in a planar graphene superlattices (PGSL) formed by patterning graphene sheet on laterally hetrostructured substrate as Sio2/hBN. A tunable optical transitions between minibands is observed based on engineering structural parameters. We derive analytically expression for optical absorption from two-band model. Considerable optical absorption is obtained for different ratios between widths of heterostructured substrate and is explained analytically from the view point of wavefunction engineering and miniband dispersion, in details. The role of different statuses of polarization as circular and linear are considered. Our study paves a way toward the control of optical properties of PGSLs to be implemented in optoelectronics devices.

Pressure and wall shear stress in blood hammer - Analytical theory.

PubMed

Mei, Chiang C; Jing, Haixiao

2016-10-01

We describe an analytical theory of blood hammer in a long and stiffened artery due to sudden blockage. Based on the model of a viscous fluid in laminar flow, we derive explicit expressions of oscillatory pressure and wall shear stress. To examine the effects on local plaque formation we also allow the blood vessel radius to be slightly nonuniform. Without resorting to discrete computation, the asymptotic method of multiple scales is utilized to deal with the sharp contrast of time scales. The effects of plaque and blocking time on blood pressure and wall shear stress are studied. The theory is validated by comparison with existing water hammer experiments. Copyright © 2016. Published by Elsevier Inc.

Low-frequency quadrupole impedance of undulators and wigglers

DOE PAGES

Blednykh, A.; Bassi, G.; Hidaka, Y.; ...

2016-10-25

An analytical expression of the low-frequency quadrupole impedance for undulators and wigglers is derived and benchmarked against beam-based impedance measurements done at the 3 GeV NSLS-II storage ring. The adopted theoretical model, valid for an arbitrary number of electromagnetic layers with parallel geometry, allows to calculate the quadrupole impedance for arbitrary values of the magnetic permeability μ r. Here, in the comparison of the analytical results with the measurements for variable magnet gaps, two limit cases of the permeability have been studied: the case of perfect magnets (μ r → ∞), and the case in which the magnets are fullymore » saturated (μ r = 1).« less

Analytic model for ultrasound energy receivers and their optimal electric loads

NASA Astrophysics Data System (ADS)

Gorostiaga, M.; Wapler, M. C.; Wallrabe, U.

2017-08-01

In this paper, we present an analytic model for thickness resonating plate ultrasound energy receivers, which we have derived from the piezoelectric and the wave equations and, in which we have included dielectric, viscosity and acoustic attenuation losses. Afterwards, we explore the optimal electric load predictions by the zero reflection and power maximization approaches present in the literature with different acoustic boundary conditions, and discuss their limitations. To validate our model, we compared our expressions with the KLM model solved numerically with very good agreement. Finally, we discuss the differences between the zero reflection and power maximization optimal electric loads, which start to differ as losses in the receiver increase.

A general statistical test for correlations in a finite-length time series.

PubMed

Hanson, Jeffery A; Yang, Haw

2008-06-07

The statistical properties of the autocorrelation function from a time series composed of independently and identically distributed stochastic variables has been studied. Analytical expressions for the autocorrelation function's variance have been derived. It has been found that two common ways of calculating the autocorrelation, moving-average and Fourier transform, exhibit different uncertainty characteristics. For periodic time series, the Fourier transform method is preferred because it gives smaller uncertainties that are uniform through all time lags. Based on these analytical results, a statistically robust method has been proposed to test the existence of correlations in a time series. The statistical test is verified by computer simulations and an application to single-molecule fluorescence spectroscopy is discussed.

GRB afterglows in the nonrelativistic phase

NASA Astrophysics Data System (ADS)

Huang, Y. F.; Lu, T.

2008-10-01

When discussing the afterglows of gamma-ray bursts analytically, it is usually assumed that the external shock is always ultra-relativisitc, with the bulk Lorentz factor much larger than 1. However, we show that the deceleration of the external shock is actually very quick. The afterglow may typically enter the nonrelativistic phase in several days to teens of days, and may even enter the deep Newtonian phase in tens of days to several months. One thus should be careful in using those familiar analytical expressions that are derived only under the ultra-relativistic assumption. To explain the observed afterglows that typically last for a few weeks to several months, we need to consider the dynamics and radiation in the nonrelativisitic phase.

Modeling and analysis of a magnetically levitated synchronous permanent magnet planar motor

NASA Astrophysics Data System (ADS)

Kou, Baoquan; Zhang, Lu; Li, Liyi; Zhang, Hailin

2012-04-01

In this paper, a new magnetically levitated synchronous permanent magnet planar motor (MLSPMPM) driven by composite-current is proposed, of which the mover is made of a copper coil array and the stator are magnets and magnetic conductor. The coil pitch τt and permanent magnet pole pitch τp satisfy the following relationship 3nτt = (3n ± 1)τp. Firstly, an analytical model of the planar motor is established, flux density distribution of the two-dimensional magnet array is obtained by solving the equations of the scalar magnetic potential. Secondly, the expressions of the electromagnetic forces induced by magnetic field and composite current are derived. To verify the analytical model and the electromagnetic forces, finite element method (FEM) is used for calculating the flux density and electromagnetic forces of the MLSPMPM. And the results from FEM are in good agreement with the results from the analytical equations. This indicates that the analytical model is reasonable.

Analytical Solution for the Anisotropic Rabi Model: Effects of Counter-Rotating Terms

NASA Astrophysics Data System (ADS)

Zhang, Guofeng; Zhu, Hanjie

2015-03-01

The anisotropic Rabi model, which was proposed recently, differs from the original Rabi model: the rotating and counter-rotating terms are governed by two different coupling constants. This feature allows us to vary the counter-rotating interaction independently and explore the effects of it on some quantum properties. In this paper, we eliminate the counter-rotating terms approximately and obtain the analytical energy spectrums and wavefunctions. These analytical results agree well with the numerical calculations in a wide range of the parameters including the ultrastrong coupling regime. In the weak counter-rotating coupling limit we find out that the counter-rotating terms can be considered as the shifts to the parameters of the Jaynes-Cummings model. This modification shows the validness of the rotating-wave approximation on the assumption of near-resonance and relatively weak coupling. Moreover, the analytical expressions of several physics quantities are also derived, and the results show the break-down of the U(1)-symmetry and the deviation from the Jaynes-Cummings model.

Analytical solution for the anisotropic Rabi model: effects of counter-rotating terms.

PubMed

Zhang, Guofeng; Zhu, Hanjie

2015-03-04

The anisotropic Rabi model, which was proposed recently, differs from the original Rabi model: the rotating and counter-rotating terms are governed by two different coupling constants. This feature allows us to vary the counter-rotating interaction independently and explore the effects of it on some quantum properties. In this paper, we eliminate the counter-rotating terms approximately and obtain the analytical energy spectrums and wavefunctions. These analytical results agree well with the numerical calculations in a wide range of the parameters including the ultrastrong coupling regime. In the weak counter-rotating coupling limit we find out that the counter-rotating terms can be considered as the shifts to the parameters of the Jaynes-Cummings model. This modification shows the validness of the rotating-wave approximation on the assumption of near-resonance and relatively weak coupling. Moreover, the analytical expressions of several physics quantities are also derived, and the results show the break-down of the U(1)-symmetry and the deviation from the Jaynes-Cummings model.

Analytical Solution for the Anisotropic Rabi Model: Effects of Counter-Rotating Terms

PubMed Central

Zhang, Guofeng; Zhu, Hanjie

2015-01-01

The anisotropic Rabi model, which was proposed recently, differs from the original Rabi model: the rotating and counter-rotating terms are governed by two different coupling constants. This feature allows us to vary the counter-rotating interaction independently and explore the effects of it on some quantum properties. In this paper, we eliminate the counter-rotating terms approximately and obtain the analytical energy spectrums and wavefunctions. These analytical results agree well with the numerical calculations in a wide range of the parameters including the ultrastrong coupling regime. In the weak counter-rotating coupling limit we find out that the counter-rotating terms can be considered as the shifts to the parameters of the Jaynes-Cummings model. This modification shows the validness of the rotating-wave approximation on the assumption of near-resonance and relatively weak coupling. Moreover, the analytical expressions of several physics quantities are also derived, and the results show the break-down of the U(1)-symmetry and the deviation from the Jaynes-Cummings model. PMID:25736827

Dual metal gate tunneling field effect transistors based on MOSFETs: A 2-D analytical approach

NASA Astrophysics Data System (ADS)

Ramezani, Zeinab; Orouji, Ali A.

2018-01-01

A novel 2-D analytical drain current model of novel Dual Metal Gate Tunnel Field Effect Transistors Based on MOSFETs (DMG-TFET) is presented in this paper. The proposed Tunneling FET is extracted from a MOSFET structure by employing an additional electrode in the source region with an appropriate work function to induce holes in the N+ source region and hence makes it as a P+ source region. The electric field is derived which is utilized to extract the expression of the drain current by analytically integrating the band to band tunneling generation rate in the tunneling region based on the potential profile by solving the Poisson's equation. Through this model, the effects of the thin film thickness and gate voltage on the potential, the electric field, and the effects of the thin film thickness on the tunneling current can be studied. To validate our present model we use SILVACO ATLAS device simulator and the analytical results have been compared with it and found a good agreement.

A one-step method for modelling longitudinal data with differential equations.

PubMed

Hu, Yueqin; Treinen, Raymond

2018-04-06

Differential equation models are frequently used to describe non-linear trajectories of longitudinal data. This study proposes a new approach to estimate the parameters in differential equation models. Instead of estimating derivatives from the observed data first and then fitting a differential equation to the derivatives, our new approach directly fits the analytic solution of a differential equation to the observed data, and therefore simplifies the procedure and avoids bias from derivative estimations. A simulation study indicates that the analytic solutions of differential equations (ASDE) approach obtains unbiased estimates of parameters and their standard errors. Compared with other approaches that estimate derivatives first, ASDE has smaller standard error, larger statistical power and accurate Type I error. Although ASDE obtains biased estimation when the system has sudden phase change, the bias is not serious and a solution is also provided to solve the phase problem. The ASDE method is illustrated and applied to a two-week study on consumers' shopping behaviour after a sale promotion, and to a set of public data tracking participants' grammatical facial expression in sign language. R codes for ASDE, recommendations for sample size and starting values are provided. Limitations and several possible expansions of ASDE are also discussed. © 2018 The British Psychological Society.

A generalized preferential attachment model for business firms growth rates. II. Mathematical treatment

NASA Astrophysics Data System (ADS)

Buldyrev, S. V.; Pammolli, F.; Riccaboni, M.; Yamasaki, K.; Fu, D.-F.; Matia, K.; Stanley, H. E.

2007-05-01

We present a preferential attachment growth model to obtain the distribution P(K) of number of units K in the classes which may represent business firms or other socio-economic entities. We found that P(K) is described in its central part by a power law with an exponent ϕ = 2+b/(1-b) which depends on the probability of entry of new classes, b. In a particular problem of city population this distribution is equivalent to the well known Zipf law. In the absence of the new classes entry, the distribution P(K) is exponential. Using analytical form of P(K) and assuming proportional growth for units, we derive P(g), the distribution of business firm growth rates. The model predicts that P(g) has a Laplacian cusp in the central part and asymptotic power-law tails with an exponent ζ = 3. We test the analytical expressions derived using heuristic arguments by simulations. The model might also explain the size-variance relationship of the firm growth rates.

Hollow vortex Gaussian beams

NASA Astrophysics Data System (ADS)

Zhou, GuoQuan; Cai, YangJian; Dai, ChaoQing

2013-05-01

A kind of hollow vortex Gaussian beam is introduced. Based on the Collins integral, an analytical propagation formula of a hollow vortex Gaussian beam through a paraxial ABCD optical system is derived. Due to the special distribution of the optical field, which is caused by the initial vortex phase, the dark region of a hollow vortex Gaussian beam will not disappear upon propagation. The analytical expressions for the beam propagation factor, the kurtosis parameter, and the orbital angular momentum density of a hollow vortex Gaussian beam passing through a paraxial ABCD optical system are also derived, respectively. The beam propagation factor is determined by the beam order and the topological charge. The kurtosis parameter and the orbital angular momentum density depend on beam order n, topological charge m, parameter γ, and transfer matrix elements A and D. As a numerical example, the propagation properties of a hollow vortex Gaussian beam in free space are demonstrated. The hollow vortex Gaussian beam has eminent propagation stability and has crucial application prospects in optical micromanipulation.

Padé Approximant and Minimax Rational Approximation in Standard Cosmology

NASA Astrophysics Data System (ADS)

Zaninetti, Lorenzo

2016-02-01

The luminosity distance in the standard cosmology as given by $\\Lambda$CDM and consequently the distance modulus for supernovae can be defined by the Pad\\'e approximant. A comparison with a known analytical solution shows that the Pad\\'e approximant for the luminosity distance has an error of $4\\%$ at redshift $= 10$. A similar procedure for the Taylor expansion of the luminosity distance gives an error of $4\\%$ at redshift $=0.7 $; this means that for the luminosity distance, the Pad\\'e approximation is superior to the Taylor series. The availability of an analytical expression for the distance modulus allows applying the Levenberg--Marquardt method to derive the fundamental parameters from the available compilations for supernovae. A new luminosity function for galaxies derived from the truncated gamma probability density function models the observed luminosity function for galaxies when the observed range in absolute magnitude is modeled by the Pad\\'e approximant. A comparison of $\\Lambda$CDM with other cosmologies is done adopting a statistical point of view.

Modal sound transmission loss of a single leaf panel: Effects of inter-modal coupling.

PubMed

Wang, Chong

2015-06-01

Sound transmission through a single leaf panel has mostly been discussed and explained by using the approaching wave concept, from which the well-known mass law can be derived. In this paper, the modal behavior in sound transmission coefficients is explored, and it is shown that the mutual modal radiation impedances in modal sound transmission coefficients may not be ignored even for a panel immersed in a light fluid. By introducing the equivalent modal impedance which incorporates the inter-modal coupling effect, an analytical expression for the modal sound transmission coefficient is derived, and the overall sound transmission coefficient is simply a modal superposition of modal sound transmission coefficients. A good correlation is obtained between analytical calculation and boundary element method. In addition, it is found that inter-modal coupling has noticeable effects in modal sound transmission coefficients in the subsonic region but may be ignored as modes become supersonic. It is also shown that the well-known mass law performance is attributed to all the supersonic modes.

Scaling analysis and instantons for thermally assisted tunneling and quantum Monte Carlo simulations

NASA Astrophysics Data System (ADS)

Jiang, Zhang; Smelyanskiy, Vadim N.; Isakov, Sergei V.; Boixo, Sergio; Mazzola, Guglielmo; Troyer, Matthias; Neven, Hartmut

2017-01-01

We develop an instantonic calculus to derive an analytical expression for the thermally assisted tunneling decay rate of a metastable state in a fully connected quantum spin model. The tunneling decay problem can be mapped onto the Kramers escape problem of a classical random dynamical field. This dynamical field is simulated efficiently by path-integral quantum Monte Carlo (QMC). We show analytically that the exponential scaling with the number of spins of the thermally assisted quantum tunneling rate and the escape rate of the QMC process are identical. We relate this effect to the existence of a dominant instantonic tunneling path. The instanton trajectory is described by nonlinear dynamical mean-field theory equations for a single-site magnetization vector, which we solve exactly. Finally, we derive scaling relations for the "spiky" barrier shape when the spin tunneling and QMC rates scale polynomially with the number of spins N while a purely classical over-the-barrier activation rate scales exponentially with N .

The analysis of non-linear dynamic behavior (including snap-through) of postbuckled plates by simple analytical solution

NASA Technical Reports Server (NTRS)

Ng, C. F.

1988-01-01

Static postbuckling and nonlinear dynamic analysis of plates are usually accomplished by multimode analyses, although the methods are complicated and do not give straightforward understanding of the nonlinear behavior. Assuming single-mode transverse displacement, a simple formula is derived for the transverse load displacement relationship of a plate under in-plane compression. The formula is used to derive a simple analytical expression for the static postbuckling displacement and nonlinear dynamic responses of postbuckled plates under sinusoidal or random excitation. Regions with softening and hardening spring behavior are identified. Also, the highly nonlinear motion of snap-through and its effects on the overall dynamic response can be easily interpreted using the single-mode formula. Theoretical results are compared with experimental results obtained using a buckled aluminum panel, using discrete frequency and broadband point excitation. Some important effects of the snap-through motion on the dynamic response of the postbuckled plates are found.

Concise expression of a classical radiation spectrum

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wang, C.

1993-06-01

In this paper we present a concise expression of the classical electromagnetic radiation spectrum of a moving charge. It is shown to be equivalent to the often used and much more complicated form derived from the Lienard-Wiechert potentials when the observation distance [ital R] satisfies the condition [ital R][much gt][gamma][lambda]. The expression reveals a relationship between the radiation spectrum and the motion of the radiation source. It also forms the basis of an efficient computing approach, which is of practical value in numerical calculations of the spectral output of accelerated charges. The advantages of this approach for analytical and numericalmore » applications are discussed and the bending-magnet synchrotron radiation spectrum is calculated according to the approach.« less

SER Analysis of MPPM-Coded MIMO-FSO System over Uncorrelated and Correlated Gamma-Gamma Atmospheric Turbulence Channels

NASA Astrophysics Data System (ADS)

Khallaf, Haitham S.; Garrido-Balsells, José M.; Shalaby, Hossam M. H.; Sampei, Seiichi

2015-12-01

The performance of multiple-input multiple-output free space optical (MIMO-FSO) communication systems, that adopt multipulse pulse position modulation (MPPM) techniques, is analyzed. Both exact and approximate symbol-error rates (SERs) are derived for both cases of uncorrelated and correlated channels. The effects of background noise, receiver shot-noise, and atmospheric turbulence are taken into consideration in our analysis. The random fluctuations of the received optical irradiance, produced by the atmospheric turbulence, is modeled by the widely used gamma-gamma statistical distribution. Uncorrelated MIMO channels are modeled by the α-μ distribution. A closed-form expression for the probability density function of the optical received irradiance is derived for the case of correlated MIMO channels. Using our analytical expressions, the degradation of the system performance with the increment of the correlation coefficients between MIMO channels is corroborated.

Perihelion precession from power law central force and magnetic-like force

NASA Astrophysics Data System (ADS)

Xu, Feng

2011-04-01

By the Laplace-Runge-Lenz (LRL) vector, we analyzed perihelion precessions of orbit with arbitrary eccentricity from perturbations of 1) power law central force and 2) fThusmagnetic-like force. Exact and analytically closed expressions for the precession rate are derived in both cases. In the central force case, we give a further expansion expression of precession rate in orders of eccentricity, and a rule judging pro- or retrograde precession is also given. We applied the result of central force to precessions of a planet in 1) Schwarzschild space-time, for which the formula for the Mercury’s 43”/century is reproduced, and 2) spherically distributed dark matter, for which we find a formula that is a generalization of the result derived by others for circular orbit. In the magnetic case, the use of the LRL vector proves to be simple and efficient. An example of magnetic-like perturbation is also discussed.

Perihelion precession from power law central force and magnetic-like force

DOE Office of Scientific and Technical Information (OSTI.GOV)

Xu Feng

2011-04-15

By the Laplace-Runge-Lenz (LRL) vector, we analyzed perihelion precessions of orbit with arbitrary eccentricity from perturbations of 1) power law central force and 2) fThusmagnetic-like force. Exact and analytically closed expressions for the precession rate are derived in both cases. In the central force case, we give a further expansion expression of precession rate in orders of eccentricity, and a rule judging pro- or retrograde precession is also given. We applied the result of central force to precessions of a planet in 1) Schwarzschild space-time, for which the formula for the Mercury's 43''/century is reproduced, and 2) spherically distributed darkmore » matter, for which we find a formula that is a generalization of the result derived by others for circular orbit. In the magnetic case, the use of the LRL vector proves to be simple and efficient. An example of magnetic-like perturbation is also discussed.« less

Covariance Bell inequalities

NASA Astrophysics Data System (ADS)

Pozsgay, Victor; Hirsch, Flavien; Branciard, Cyril; Brunner, Nicolas

2017-12-01

We introduce Bell inequalities based on covariance, one of the most common measures of correlation. Explicit examples are discussed, and violations in quantum theory are demonstrated. A crucial feature of these covariance Bell inequalities is their nonlinearity; this has nontrivial consequences for the derivation of their local bound, which is not reached by deterministic local correlations. For our simplest inequality, we derive analytically tight bounds for both local and quantum correlations. An interesting application of covariance Bell inequalities is that they can act as "shared randomness witnesses": specifically, the value of the Bell expression gives device-independent lower bounds on both the dimension and the entropy of the shared random variable in a local model.

Scattering Manipulation and Camouflage of Electrically Small Objects through Metasurfaces

NASA Astrophysics Data System (ADS)

Vellucci, S.; Monti, A.; Toscano, A.; Bilotti, F.

2017-03-01

In this paper, we discuss the intriguing possibility of tailoring the scattering response of an electrically small object for camouflage and illusion applications using metasurfaces. As a significant example, we focus our attention on the cylindrical geometry and derive the analytical conditions needed to camouflage the geometrical and electrical characteristics of dielectric and metallic cylinders coated with ideal metasurfaces. A closed-form expression of the camouflaging metasurface depending on the cylinder's characteristics is derived. Furthermore, the frequency behavior and the limitations of this technique are discussed with the aid of relevant examples. In order to overcome these limitations, a solution based on the use of lossy metasurfaces is proposed.

Expressions for tidal conversion at seafloor topography using physical space integrals

NASA Astrophysics Data System (ADS)

Schorghofer, Norbert

2010-12-01

The barotropic tide interacts with seafloor topography to generate internal gravity waves. Equations for streamfunction and power conversion are derived in terms of integrals over the topography in spatial coordinates. The slope of the topography does not need to be small. Explicit equations are derived up to second order in slope for general topography, and conversion by a bell-shaped topography is calculated analytically to this order. A concise formalism using Hilbert transforms is developed, the minimally converting topographic shape is discussed, and a numerical scheme for the evaluation of power conversion is designed that robustly deals with the singular integrand.

Stress analysis and failure of an internally pressurized composite-jacketed steel cylinder

NASA Technical Reports Server (NTRS)

Chen, Peter C. T.

1992-01-01

This paper presents a nonlinear stress analysis of a thick-walled compound tube subjected to internal pressure. The compound tube is constructed of a steel liner and a graphite-bismaleimide outer shell. Analytical expressions for the stresses, strains, and displacements are derived for all loading ranges up to failure. Numerical results for the stresses and the maximum value that the compound tube can contain without failure are presented.

Radiative transfer in falling snow: A two-stream approximation

NASA Astrophysics Data System (ADS)

Koh, Gary

1989-04-01

Light transmission measurements through falling snow have produced results unexplainable by single scattering arguments. A two-stream approximation to radiative transfer is used to derive an analytical expression that describes the effects of multiple scattering as a function of the snow optical depth and the snow asymmetry parameter. The approximate solution is simple and it may be as accurate as the exact solution for describing the transmission measurements within the limits of experimental uncertainties.

Equilibrium Free Energies from Nonequilibrium Metadynamics

NASA Astrophysics Data System (ADS)

Bussi, Giovanni; Laio, Alessandro; Parrinello, Michele

2006-03-01

In this Letter we propose a new formalism to map history-dependent metadynamics in a Markovian process. We apply this formalism to model Langevin dynamics and determine the equilibrium distribution of a collection of simulations. We demonstrate that the reconstructed free energy is an unbiased estimate of the underlying free energy and analytically derive an expression for the error. The present results can be applied to other history-dependent stochastic processes, such as Wang-Landau sampling.

Analytical expressions for the closure probability of a stiff wormlike chain for finite capture radius.

PubMed

Guérin, T

2017-08-01

Estimating the probability that two monomers of the same polymer chain are close together is a key ingredient to characterize intramolecular reactions and polymer looping. In the case of stiff wormlike polymers (rigid fluctuating elastic rods), for which end-to-end encounters are rare events, we derive an explicit analytical formula for the probability η(r_{c}) that the distance between the chain extremities is smaller than some capture radius r_{c}. The formula is asymptotically exact in the limit of stiff chains, and it leads to the identification of two distinct scaling regimes for the closure factor, originating from a strong variation of the fluctuations of the chain orientation at closure. Our theory is compatible with existing analytical results from the literature that cover the cases of a vanishing capture radius and of nearly fully extended chains.

High-temperature ratchets with sawtooth potentials

NASA Astrophysics Data System (ADS)

Rozenbaum, Viktor M.; Shapochkina, Irina V.; Sheu, Sheh-Yi; Yang, Dah-Yen; Lin, Sheng Hsien

2016-11-01

The concept of the effective potential is suggested as an efficient instrument to get a uniform analytical description of stochastic high-temperature on-off flashing and rocking ratchets. The analytical representation for the average particle velocity, obtained within this technique, allows description of ratchets with sharp potentials (and potentials with jumps in particular). For sawtooth potentials, the explicit analytical expressions for the average velocity of on-off flashing and rocking ratchets valid for arbitrary frequencies of potential energy fluctuations are derived; the difference in their high-frequency asymptotics is explored for the smooth and cusped profiles, and profiles with jumps. The origin of the difference as well as the appearance of the jump behavior in ratchet characteristics are interpreted in terms of self-similar universal solutions which give the continuous description of the effect. It is shown how the jump behavior in motor characteristics arises from the competition between the characteristic times of the system.

Analysis of Mathematical Modelling on Potentiometric Biosensors

PubMed Central

Mehala, N.; Rajendran, L.

2014-01-01

A mathematical model of potentiometric enzyme electrodes for a nonsteady condition has been developed. The model is based on the system of two coupled nonlinear time-dependent reaction diffusion equations for Michaelis-Menten formalism that describes the concentrations of substrate and product within the enzymatic layer. Analytical expressions for the concentration of substrate and product and the corresponding flux response have been derived for all values of parameters using the new homotopy perturbation method. Furthermore, the complex inversion formula is employed in this work to solve the boundary value problem. The analytical solutions obtained allow a full description of the response curves for only two kinetic parameters (unsaturation/saturation parameter and reaction/diffusion parameter). Theoretical descriptions are given for the two limiting cases (zero and first order kinetics) and relatively simple approaches for general cases are presented. All the analytical results are compared with simulation results using Scilab/Matlab program. The numerical results agree with the appropriate theories. PMID:25969765

Analysis of mathematical modelling on potentiometric biosensors.

PubMed

Mehala, N; Rajendran, L

2014-01-01

A mathematical model of potentiometric enzyme electrodes for a nonsteady condition has been developed. The model is based on the system of two coupled nonlinear time-dependent reaction diffusion equations for Michaelis-Menten formalism that describes the concentrations of substrate and product within the enzymatic layer. Analytical expressions for the concentration of substrate and product and the corresponding flux response have been derived for all values of parameters using the new homotopy perturbation method. Furthermore, the complex inversion formula is employed in this work to solve the boundary value problem. The analytical solutions obtained allow a full description of the response curves for only two kinetic parameters (unsaturation/saturation parameter and reaction/diffusion parameter). Theoretical descriptions are given for the two limiting cases (zero and first order kinetics) and relatively simple approaches for general cases are presented. All the analytical results are compared with simulation results using Scilab/Matlab program. The numerical results agree with the appropriate theories.

The exact solutions and approximate analytic solutions of the (2 + 1)-dimensional KP equation based on symmetry method.

PubMed

Gai, Litao; Bilige, Sudao; Jie, Yingmo

2016-01-01

In this paper, we successfully obtained the exact solutions and the approximate analytic solutions of the (2 + 1)-dimensional KP equation based on the Lie symmetry, the extended tanh method and the homotopy perturbation method. In first part, we obtained the symmetries of the (2 + 1)-dimensional KP equation based on the Wu-differential characteristic set algorithm and reduced it. In the second part, we constructed the abundant exact travelling wave solutions by using the extended tanh method. These solutions are expressed by the hyperbolic functions, the trigonometric functions and the rational functions respectively. It should be noted that when the parameters are taken as special values, some solitary wave solutions are derived from the hyperbolic function solutions. Finally, we apply the homotopy perturbation method to obtain the approximate analytic solutions based on four kinds of initial conditions.

Excitation of low-frequency residual currents at combination frequencies of an ionising two-colour laser pulse

NASA Astrophysics Data System (ADS)

Vvedenskii, N. V.; Kostin, V. A.; Laryushin, I. D.; Silaev, A. A.

2016-05-01

We have studied the processes of excitation of low-frequency residual currents in a plasma produced through ionisation of gases by two-colour laser pulses in laser-plasma schemes for THz generation. We have developed an analytical approach that allows one to find residual currents in the case when one of the components of a two-colour pulse is weak enough. The derived analytical expressions show that the effective generation of the residual current (and hence the effective THz generation) is possible if the ratio of the frequencies in the two-colour laser pulse is close to a rational fraction with a not very big odd sum of the numerator and denominator. The results of numerical calculations (including those based on the solution of the three-dimensional time-dependent Schrödinger equation) agree well with the analytical results.

Quantitative Characterization of Tissue Microstructure with Temporal Diffusion Spectroscopy

PubMed Central

Xu, Junzhong; Does, Mark D.; Gore, John C.

2009-01-01

The signals recorded by diffusion-weighted magnetic resonance imaging (DWI) are dependent on the micro-structural properties of biological tissues, so it is possible to obtain quantitative structural information non-invasively from such measurements. Oscillating gradient spin echo (OGSE) methods have the ability to probe the behavior of water diffusion over different time scales and the potential to detect variations in intracellular structure. To assist in the interpretation of OGSE data, analytical expressions have been derived for diffusion-weighted signals with OGSE methods for restricted diffusion in some typical structures, including parallel planes, cylinders and spheres, using the theory of temporal diffusion spectroscopy. These analytical predictions have been confirmed with computer simulations. These expressions suggest how OGSE signals from biological tissues should be analyzed to characterize tissue microstructure, including how to estimate cell nuclear sizes. This approach provides a model to interpret diffusion data obtained from OGSE measurements that can be used for applications such as monitoring tumor response to treatment in vivo. PMID:19616979

Variable Refractive Index Effects on Radiation in Semitransparent Scattering Multilayered Regions

NASA Technical Reports Server (NTRS)

Siegel, R.; Spuckler, C. M.

1993-01-01

A simple set of equations is derived for predicting the temperature distribution and radiative energy flow in a semitransparent layer consisting of an arbitrary number of laminated sublayers that absorb, emit, and scatter radiation. Each sublayer can have a different refractive index and optical thickness. The plane composite region is heated on each exterior side by a different amount of incident radiation. The results are for the limiting case where heat conduction within the layers is very small relative to radiative transfer, and is neglected. The interfaces are assumed diffuse, and all interface reflections are included in the analysis. The thermal behavior is readily calculated from the analytical expressions that are obtained. By using many sublayers, the analytical expressions provide the temperature distribution and heat flow for a diffusing medium with a continuously varying refractive index, including internal reflection effects caused by refractive index gradients. Temperature and heat flux results are given to show the effect of variations in refractive index and optical thickness through the multilayer laminate.

Analytical description of high-aperture STED resolution with 0–2π vortex phase modulation

PubMed Central

Xie, Hao; Liu, Yujia; Jin, Dayong; Santangelo, Philip J.; Xi, Peng

2014-01-01

Stimulated emission depletion (STED) can achieve optical superresolution, with the optical diffraction limit broken by the suppression on the periphery of the fluorescent focal spot. Previously, it is generally experimentally accepted that there exists an inverse square root relationship with the STED power and the resolution, but with arbitrary coefficients in expression. In this paper, we have removed the arbitrary coefficients by exploring the relationship between the STED power and the achievable resolution from vector optical theory for the widely used 0–2π vortex phase modulation. Electromagnetic fields of the focal region of a high numerical aperture objective are calculated and approximated into polynomials of radius in the focal plane, and analytical expression of resolution as a function of the STED intensity has been derived. As a result, the resolution can be estimated directly from the measurement of the saturation power of the dye and the STED power applied in the region of high STED power. PMID:24323224

The spectrum of a vertex model and related spin one chain sitting in a genus five curve

NASA Astrophysics Data System (ADS)

Martins, M. J.

2017-11-01

We derive the transfer matrix eigenvalues of a three-state vertex model whose weights are based on a R-matrix not of difference form with spectral parameters lying on a genus five curve. We have shown that the basic building blocks for both the transfer matrix eigenvalues and Bethe equations can be expressed in terms of meromorphic functions on an elliptic curve. We discuss the properties of an underlying spin one chain originated from a particular choice of the R-matrix second spectral parameter. We present numerical and analytical evidences that the respective low-energy excitations can be gapped or massless depending on the strength of the interaction coupling. In the massive phase we provide analytical and numerical evidences in favor of an exact expression for the lowest energy gap. We point out that the critical point separating these two distinct physical regimes coincides with the one in which the weights geometry degenerate into union of genus one curves.

Sex genes for genomic analysis in human brain: internal controls for comparison of probe level data extraction.

PubMed Central

Galfalvy, Hanga C; Erraji-Benchekroun, Loubna; Smyrniotopoulos, Peggy; Pavlidis, Paul; Ellis, Steven P; Mann, J John; Sibille, Etienne; Arango, Victoria

2003-01-01

Background Genomic studies of complex tissues pose unique analytical challenges for assessment of data quality, performance of statistical methods used for data extraction, and detection of differentially expressed genes. Ideally, to assess the accuracy of gene expression analysis methods, one needs a set of genes which are known to be differentially expressed in the samples and which can be used as a "gold standard". We introduce the idea of using sex-chromosome genes as an alternative to spiked-in control genes or simulations for assessment of microarray data and analysis methods. Results Expression of sex-chromosome genes were used as true internal biological controls to compare alternate probe-level data extraction algorithms (Microarray Suite 5.0 [MAS5.0], Model Based Expression Index [MBEI] and Robust Multi-array Average [RMA]), to assess microarray data quality and to establish some statistical guidelines for analyzing large-scale gene expression. These approaches were implemented on a large new dataset of human brain samples. RMA-generated gene expression values were markedly less variable and more reliable than MAS5.0 and MBEI-derived values. A statistical technique controlling the false discovery rate was applied to adjust for multiple testing, as an alternative to the Bonferroni method, and showed no evidence of false negative results. Fourteen probesets, representing nine Y- and two X-chromosome linked genes, displayed significant sex differences in brain prefrontal cortex gene expression. Conclusion In this study, we have demonstrated the use of sex genes as true biological internal controls for genomic analysis of complex tissues, and suggested analytical guidelines for testing alternate oligonucleotide microarray data extraction protocols and for adjusting multiple statistical analysis of differentially expressed genes. Our results also provided evidence for sex differences in gene expression in the brain prefrontal cortex, supporting the notion of a putative direct role of sex-chromosome genes in differentiation and maintenance of sexual dimorphism of the central nervous system. Importantly, these analytical approaches are applicable to all microarray studies that include male and female human or animal subjects. PMID:12962547

Sex genes for genomic analysis in human brain: internal controls for comparison of probe level data extraction.

PubMed

Galfalvy, Hanga C; Erraji-Benchekroun, Loubna; Smyrniotopoulos, Peggy; Pavlidis, Paul; Ellis, Steven P; Mann, J John; Sibille, Etienne; Arango, Victoria

2003-09-08

Genomic studies of complex tissues pose unique analytical challenges for assessment of data quality, performance of statistical methods used for data extraction, and detection of differentially expressed genes. Ideally, to assess the accuracy of gene expression analysis methods, one needs a set of genes which are known to be differentially expressed in the samples and which can be used as a "gold standard". We introduce the idea of using sex-chromosome genes as an alternative to spiked-in control genes or simulations for assessment of microarray data and analysis methods. Expression of sex-chromosome genes were used as true internal biological controls to compare alternate probe-level data extraction algorithms (Microarray Suite 5.0 [MAS5.0], Model Based Expression Index [MBEI] and Robust Multi-array Average [RMA]), to assess microarray data quality and to establish some statistical guidelines for analyzing large-scale gene expression. These approaches were implemented on a large new dataset of human brain samples. RMA-generated gene expression values were markedly less variable and more reliable than MAS5.0 and MBEI-derived values. A statistical technique controlling the false discovery rate was applied to adjust for multiple testing, as an alternative to the Bonferroni method, and showed no evidence of false negative results. Fourteen probesets, representing nine Y- and two X-chromosome linked genes, displayed significant sex differences in brain prefrontal cortex gene expression. In this study, we have demonstrated the use of sex genes as true biological internal controls for genomic analysis of complex tissues, and suggested analytical guidelines for testing alternate oligonucleotide microarray data extraction protocols and for adjusting multiple statistical analysis of differentially expressed genes. Our results also provided evidence for sex differences in gene expression in the brain prefrontal cortex, supporting the notion of a putative direct role of sex-chromosome genes in differentiation and maintenance of sexual dimorphism of the central nervous system. Importantly, these analytical approaches are applicable to all microarray studies that include male and female human or animal subjects.

Analytic second derivative of the energy for density functional theory based on the three-body fragment molecular orbital method

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nakata, Hiroya, E-mail: nakata.h.ab@m.titech.ac.jp; RIKEN, Research Cluster for Innovation, Nakamura Lab, 2-1 Hirosawa, Wako, Saitama 351-0198; Japan Society for the Promotion of Science, Kojimachi Business Center Building, 5-3-1 Kojimachi, Chiyoda-ku, Tokyo 102-0083

2015-03-28

Analytic second derivatives of the energy with respect to nuclear coordinates have been developed for spin restricted density functional theory (DFT) based on the fragment molecular orbital method (FMO). The derivations were carried out for the three-body expansion (FMO3), and the two-body expressions can be obtained by neglecting the three-body corrections. Also, the restricted Hartree-Fock (RHF) Hessian for FMO3 can be obtained by neglecting the density-functional related terms. In both the FMO-RHF and FMO-DFT Hessians, certain terms with small magnitudes are neglected for computational efficiency. The accuracy of the FMO-DFT Hessian in terms of the Gibbs free energy is evaluatedmore » for a set of polypeptides and water clusters and found to be within 1 kcal/mol of the corresponding full (non-fragmented) ab initio calculation. The FMO-DFT method is also applied to transition states in S{sub N}2 reactions and for the computation of the IR and Raman spectra of a small Trp-cage protein (PDB: 1L2Y). Some computational timing analysis is also presented.« less

Towards a model of pion generalized parton distributions from Dyson-Schwinger equations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Moutarde, H.

2015-04-10

We compute the pion quark Generalized Parton Distribution H{sup q} and Double Distributions F{sup q} and G{sup q} in a coupled Bethe-Salpeter and Dyson-Schwinger approach. We use simple algebraic expressions inspired by the numerical resolution of Dyson-Schwinger and Bethe-Salpeter equations. We explicitly check the support and polynomiality properties, and the behavior under charge conjugation or time invariance of our model. We derive analytic expressions for the pion Double Distributions and Generalized Parton Distribution at vanishing pion momentum transfer at a low scale. Our model compares very well to experimental pion form factor or parton distribution function data.

An analytical model for inductively coupled implantable biomedical devices with ferrite rods.

PubMed

Theilmann, P T; Asbeck, P M

2009-02-01

Using approximations applicable to near field coupled implants simplified expressions for the complex mutual inductance of coaxial aligned coils with and without a cylindrical ferrite rod are derived. Experimental results for ferrite rods of various sizes and permeabilities are presented to verify the accuracy of this expression. An equivalent circuit model for the inductive link between an implant and power coil is then presented and used to investigate how ferrite size, permeability and loss affect the power available to the implant device. Enhancements in coupling provided by high frequency, low permeability nickel zinc rods are compared with low frequency high permeability manganese zinc rods.

Scintillation index and performance analysis of wireless optical links over non-Kolmogorov weak turbulence based on generalized atmospheric spectral model.

PubMed

Cang, Ji; Liu, Xu

2011-09-26

Based on the generalized spectral model for non-Kolmogorov atmospheric turbulence, analytic expressions of the scintillation index (SI) are derived for plane, spherical optical waves and a partially coherent Gaussian beam propagating through non-Kolmogorov turbulence horizontally in the weak fluctuation regime. The new expressions relate the SI to the finite turbulence inner and outer scales, spatial coherence of the source and spectral power-law and then used to analyze the effects of atmospheric condition and link length on the performance of wireless optical communication links. © 2011 Optical Society of America

Analysis in temporal regime of dispersive invisible structures designed from transformation optics

NASA Astrophysics Data System (ADS)

Gralak, B.; Arismendi, G.; Avril, B.; Diatta, A.; Guenneau, S.

2016-03-01

A simple invisible structure made of two anisotropic homogeneous layers is analyzed theoretically in temporal regime. The frequency dispersion is introduced and analytic expression of the transient part of the field is derived for large times when the structure is illuminated by a causal excitation. This expression shows that the limiting amplitude principle applies with transient fields decaying as the power -3 /4 of the time. The quality of the cloak is then reduced at short times and remains preserved at large times. The one-dimensional theoretical analysis is supplemented with full-wave numerical simulations in two-dimensional situations which confirm the effect of dispersion.

Analytical results for a stochastic model of gene expression with arbitrary partitioning of proteins

NASA Astrophysics Data System (ADS)

Tschirhart, Hugo; Platini, Thierry

2018-05-01

In biophysics, the search for analytical solutions of stochastic models of cellular processes is often a challenging task. In recent work on models of gene expression, it was shown that a mapping based on partitioning of Poisson arrivals (PPA-mapping) can lead to exact solutions for previously unsolved problems. While the approach can be used in general when the model involves Poisson processes corresponding to creation or degradation, current applications of the method and new results derived using it have been limited to date. In this paper, we present the exact solution of a variation of the two-stage model of gene expression (with time dependent transition rates) describing the arbitrary partitioning of proteins. The methodology proposed makes full use of the PPA-mapping by transforming the original problem into a new process describing the evolution of three biological switches. Based on a succession of transformations, the method leads to a hierarchy of reduced models. We give an integral expression of the time dependent generating function as well as explicit results for the mean, variance, and correlation function. Finally, we discuss how results for time dependent parameters can be extended to the three-stage model and used to make inferences about models with parameter fluctuations induced by hidden stochastic variables.

High-Speed Rotor Analytical Dynamics on Flexible Foundation Subjected to Internal and External Excitation

NASA Astrophysics Data System (ADS)

Jivkov, Venelin S.; Zahariev, Evtim V.

2016-12-01

The paper presents a geometrical approach to dynamics simulation of a rigid and flexible system, compiled of high speed rotating machine with eccentricity and considerable inertia and mass. The machine is mounted on a vertical flexible pillar with considerable height. The stiffness and damping of the column, as well as, of the rotor bearings and the shaft are taken into account. Non-stationary vibrations and transitional processes are analyzed. The major frequency and modal mode of the flexible column are used for analytical reduction of its mass, stiffness and damping properties. The rotor and the foundation are modelled as rigid bodies, while the flexibility of the bearings is estimated by experiments and the requirements of the manufacturer. The transition effects as a result of limited power are analyzed by asymptotic methods of averaging. Analytical expressions for the amplitudes and unstable vibrations throughout resonance are derived by quasi-static approach increasing and decreasing of the exciting frequency. Analytical functions give the possibility to analyze the influence of the design parameter of many structure applications as wind power generators, gas turbines, turbo-generators, and etc. A numerical procedure is applied to verify the effectiveness and precision of the simulation process. Nonlinear and transitional effects are analyzed and compared to the analytical results. External excitations, as wave propagation and earthquakes, are discussed. Finite elements in relative and absolute coordinates are applied to model the flexible column and the high speed rotating machine. Generalized Newton - Euler dynamics equations are used to derive the precise dynamics equations. Examples of simulation of the system vibrations and nonstationary behaviour are presented.

An Analytical Solution for the Impact of Vegetation Changes on Hydrological Partitioning Within the Budyko Framework

NASA Astrophysics Data System (ADS)

Zhang, Shulei; Yang, Yuting; McVicar, Tim R.; Yang, Dawen

2018-01-01

Vegetation change is a critical factor that profoundly affects the terrestrial water cycle. Here we derive an analytical solution for the impact of vegetation changes on hydrological partitioning within the Budyko framework. This is achieved by deriving an analytical expression between leaf area index (LAI) change and the Budyko land surface parameter (n) change, through the combination of a steady state ecohydrological model with an analytical carbon cost-benefit model for plant rooting depth. Using China where vegetation coverage has experienced dramatic changes over the past two decades as a study case, we quantify the impact of LAI changes on the hydrological partitioning during 1982-2010 and predict the future influence of these changes for the 21st century using climate model projections. Results show that LAI change exhibits an increasing importance on altering hydrological partitioning as climate becomes drier. In semiarid and arid China, increased LAI has led to substantial streamflow reductions over the past three decades (on average -8.5% in 1990s and -11.7% in 2000s compared to the 1980s baseline), and this decreasing trend in streamflow is projected to continue toward the end of this century due to predicted LAI increases. Our result calls for caution regarding the large-scale revegetation activities currently being implemented in arid and semiarid China, which may result in serious future water scarcity issues here. The analytical model developed here is physically based and suitable for simultaneously assessing both vegetation changes and climate change induced changes to streamflow globally.

Propagation of error from parameter constraints in quantitative MRI: Example application of multiple spin echo T2 mapping.

PubMed

Lankford, Christopher L; Does, Mark D

2018-02-01

Quantitative MRI may require correcting for nuisance parameters which can or must be constrained to independently measured or assumed values. The noise and/or bias in these constraints propagate to fitted parameters. For example, the case of refocusing pulse flip angle constraint in multiple spin echo T 2 mapping is explored. An analytical expression for the mean-squared error of a parameter of interest was derived as a function of the accuracy and precision of an independent estimate of a nuisance parameter. The expression was validated by simulations and then used to evaluate the effects of flip angle (θ) constraint on the accuracy and precision of T⁁2 for a variety of multi-echo T 2 mapping protocols. Constraining θ improved T⁁2 precision when the θ-map signal-to-noise ratio was greater than approximately one-half that of the first spin echo image. For many practical scenarios, constrained fitting was calculated to reduce not just the variance but the full mean-squared error of T⁁2, for bias in θ⁁≲6%. The analytical expression derived in this work can be applied to inform experimental design in quantitative MRI. The example application to T 2 mapping provided specific cases, depending on θ⁁ accuracy and precision, in which θ⁁ measurement and constraint would be beneficial to T⁁2 variance or mean-squared error. Magn Reson Med 79:673-682, 2018. © 2017 International Society for Magnetic Resonance in Medicine. © 2017 International Society for Magnetic Resonance in Medicine.

High Energy Laser Beam Propagation in the Atmosphere: The Integral Invariants of the Nonlinear Parabolic Equation and the Method of Moments

NASA Technical Reports Server (NTRS)

Manning, Robert M.

2012-01-01

The method of moments is used to define and derive expressions for laser beam deflection and beam radius broadening for high-energy propagation through the Earth s atmosphere. These expressions are augmented with the integral invariants of the corresponding nonlinear parabolic equation that describes the electric field of high-energy laser beam to propagation to yield universal equations for the aforementioned quantities; the beam deflection is a linear function of the propagation distance whereas the beam broadening is a quadratic function of distance. The coefficients of these expressions are then derived from a thin screen approximation solution of the nonlinear parabolic equation to give corresponding analytical expressions for a target located outside the Earth s atmospheric layer. These equations, which are graphically presented for a host of propagation scenarios, as well as the thin screen model, are easily amenable to the phase expansions of the wave front for the specification and design of adaptive optics algorithms to correct for the inherent phase aberrations. This work finds application in, for example, the analysis of beamed energy propulsion for space-based vehicles.

A computer code for calculations in the algebraic collective model of the atomic nucleus

NASA Astrophysics Data System (ADS)

Welsh, T. A.; Rowe, D. J.

2016-03-01

A Maple code is presented for algebraic collective model (ACM) calculations. The ACM is an algebraic version of the Bohr model of the atomic nucleus, in which all required matrix elements are derived by exploiting the model's SU(1 , 1) × SO(5) dynamical group. This paper reviews the mathematical formulation of the ACM, and serves as a manual for the code. The code enables a wide range of model Hamiltonians to be analysed. This range includes essentially all Hamiltonians that are rational functions of the model's quadrupole moments qˆM and are at most quadratic in the corresponding conjugate momenta πˆN (- 2 ≤ M , N ≤ 2). The code makes use of expressions for matrix elements derived elsewhere and newly derived matrix elements of the operators [ π ˆ ⊗ q ˆ ⊗ π ˆ ] 0 and [ π ˆ ⊗ π ˆ ] LM. The code is made efficient by use of an analytical expression for the needed SO(5)-reduced matrix elements, and use of SO(5) ⊃ SO(3) Clebsch-Gordan coefficients obtained from precomputed data files provided with the code.

Outage Probability of MRC for κ-μ Shadowed Fading Channels under Co-Channel Interference.

PubMed

Chen, Changfang; Shu, Minglei; Wang, Yinglong; Yang, Ming; Zhang, Chongqing

2016-01-01

In this paper, exact closed-form expressions are derived for the outage probability (OP) of the maximal ratio combining (MRC) scheme in the κ-μ shadowed fading channels, in which both the independent and correlated shadowing components are considered. The scenario assumes the received desired signals are corrupted by the independent Rayleigh-faded co-channel interference (CCI) and background white Gaussian noise. To this end, first, the probability density function (PDF) of the κ-μ shadowed fading distribution is obtained in the form of a power series. Then the incomplete generalized moment-generating function (IG-MGF) of the received signal-to-interference-plus-noise ratio (SINR) is derived in the closed form. By using the IG-MGF results, closed-form expressions for the OP of MRC scheme are obtained over the κ-μ shadowed fading channels. Simulation results are included to validate the correctness of the analytical derivations. These new statistical results can be applied to the modeling and analysis of several wireless communication systems, such as body centric communications.

Outage Probability of MRC for κ-μ Shadowed Fading Channels under Co-Channel Interference

PubMed Central

Chen, Changfang; Shu, Minglei; Wang, Yinglong; Yang, Ming; Zhang, Chongqing

2016-01-01

In this paper, exact closed-form expressions are derived for the outage probability (OP) of the maximal ratio combining (MRC) scheme in the κ-μ shadowed fading channels, in which both the independent and correlated shadowing components are considered. The scenario assumes the received desired signals are corrupted by the independent Rayleigh-faded co-channel interference (CCI) and background white Gaussian noise. To this end, first, the probability density function (PDF) of the κ-μ shadowed fading distribution is obtained in the form of a power series. Then the incomplete generalized moment-generating function (IG-MGF) of the received signal-to-interference-plus-noise ratio (SINR) is derived in the closed form. By using the IG-MGF results, closed-form expressions for the OP of MRC scheme are obtained over the κ-μ shadowed fading channels. Simulation results are included to validate the correctness of the analytical derivations. These new statistical results can be applied to the modeling and analysis of several wireless communication systems, such as body centric communications. PMID:27851817

Analytical model for the radio-frequency sheath

NASA Astrophysics Data System (ADS)

Czarnetzki, Uwe

2013-12-01

A simple analytical model for the planar radio-frequency (rf) sheath in capacitive discharges is developed that is based on the assumptions of a step profile for the electron front, charge exchange collisions with constant cross sections, negligible ionization within the sheath, and negligible ion dynamics. The continuity, momentum conservation, and Poisson equations are combined in a single integro-differential equation for the square of the ion drift velocity, the so called sheath equation. Starting from the kinetic Boltzmann equation, special attention is paid to the derivation and the validity of the approximate fluid equation for momentum balance. The integrals in the sheath equation appear in the screening function which considers the relative contribution of the temporal mean of the electron density to the space charge in the sheath. It is shown that the screening function is quite insensitive to variations of the effective sheath parameters. The two parameters defining the solution are the ratios of the maximum sheath extension to the ion mean free path and the Debye length, respectively. A simple general analytic expression for the screening function is introduced. By means of this expression approximate analytical solutions are obtained for the collisionless as well as the highly collisional case that compare well with the exact numerical solution. A simple transition formula allows application to all degrees of collisionality. In addition, the solutions are used to calculate all static and dynamic quantities of the sheath, e.g., the ion density, fields, and currents. Further, the rf Child-Langmuir laws for the collisionless as well as the collisional case are derived. An essential part of the model is the a priori knowledge of the wave form of the sheath voltage. This wave form is derived on the basis of a cubic charge-voltage relation for individual sheaths, considering both sheaths and the self-consistent self-bias in a discharge with arbitrary symmetry. The externally applied rf voltage is assumed to be sinusoidal, although the model can be extended to arbitrary wave forms, e.g., for dual-frequency discharges. The model calculates explicitly the cubic correction parameter in the charge-voltage relation for the case of highly asymmetric discharges. It is shown that the cubic correction is generally moderate but more pronounced in the collisionless case. The analytical results are compared to experimental data from the literature obtained by laser electric field measurements of the mean and dynamic fields in the capacitive sheath for various gases and pressures. Very good agreement is found throughout.

Analytical model for the radio-frequency sheath.

PubMed

Czarnetzki, Uwe

2013-12-01

A simple analytical model for the planar radio-frequency (rf) sheath in capacitive discharges is developed that is based on the assumptions of a step profile for the electron front, charge exchange collisions with constant cross sections, negligible ionization within the sheath, and negligible ion dynamics. The continuity, momentum conservation, and Poisson equations are combined in a single integro-differential equation for the square of the ion drift velocity, the so called sheath equation. Starting from the kinetic Boltzmann equation, special attention is paid to the derivation and the validity of the approximate fluid equation for momentum balance. The integrals in the sheath equation appear in the screening function which considers the relative contribution of the temporal mean of the electron density to the space charge in the sheath. It is shown that the screening function is quite insensitive to variations of the effective sheath parameters. The two parameters defining the solution are the ratios of the maximum sheath extension to the ion mean free path and the Debye length, respectively. A simple general analytic expression for the screening function is introduced. By means of this expression approximate analytical solutions are obtained for the collisionless as well as the highly collisional case that compare well with the exact numerical solution. A simple transition formula allows application to all degrees of collisionality. In addition, the solutions are used to calculate all static and dynamic quantities of the sheath, e.g., the ion density, fields, and currents. Further, the rf Child-Langmuir laws for the collisionless as well as the collisional case are derived. An essential part of the model is the a priori knowledge of the wave form of the sheath voltage. This wave form is derived on the basis of a cubic charge-voltage relation for individual sheaths, considering both sheaths and the self-consistent self-bias in a discharge with arbitrary symmetry. The externally applied rf voltage is assumed to be sinusoidal, although the model can be extended to arbitrary wave forms, e.g., for dual-frequency discharges. The model calculates explicitly the cubic correction parameter in the charge-voltage relation for the case of highly asymmetric discharges. It is shown that the cubic correction is generally moderate but more pronounced in the collisionless case. The analytical results are compared to experimental data from the literature obtained by laser electric field measurements of the mean and dynamic fields in the capacitive sheath for various gases and pressures. Very good agreement is found throughout.

Kubo formulas for dispersion in heterogeneous periodic nonequilibrium systems.

PubMed

Guérin, T; Dean, D S

2015-12-01

We consider the dispersion properties of tracer particles moving in nonequilibrium heterogeneous periodic media. The tracer motion is described by a Fokker-Planck equation with arbitrary spatially periodic (but constant in time) local diffusion tensors and drifts, eventually with the presence of obstacles. We derive a Kubo-like formula for the time-dependent effective diffusion tensor valid in any dimension. From this general formula, we derive expressions for the late time effective diffusion tensor and drift in these systems. In addition, we find an explicit formula for the late finite-time corrections to these transport coefficients. In one dimension, we give a closed analytical formula for the transport coefficients. The formulas derived here are very general and provide a straightforward method to compute the dispersion properties in arbitrary nonequilibrium periodic advection-diffusion systems.

Edge Fracture in Complex Fluids.

PubMed

Hemingway, Ewan J; Kusumaatmaja, Halim; Fielding, Suzanne M

2017-07-14

We study theoretically the edge fracture instability in sheared complex fluids, by means of linear stability analysis and direct nonlinear simulations. We derive an exact analytical expression for the onset of edge fracture in terms of the shear-rate derivative of the fluid's second normal stress difference, the shear-rate derivative of the shear stress, the jump in shear stress across the interface between the fluid and the outside medium (usually air), the surface tension of that interface, and the rheometer gap size. We provide a full mechanistic understanding of the edge fracture instability, carefully validated against our simulations. These findings, which are robust with respect to choice of rheological constitutive model, also suggest a possible route to mitigating edge fracture, potentially allowing experimentalists to achieve and accurately measure flows stronger than hitherto possible.

Computer simulation of the mathematical modeling involved in constitutive equation development: Via symbolic computations

NASA Technical Reports Server (NTRS)

Arnold, S. M.; Tan, H. Q.; Dong, X.

1989-01-01

Development of new material models for describing the high temperature constitutive behavior of real materials represents an important area of research in engineering disciplines. Derivation of mathematical expressions (constitutive equations) which describe this high temperature material behavior can be quite time consuming, involved and error prone; thus intelligent application of symbolic systems to facilitate this tedious process can be of significant benefit. A computerized procedure (SDICE) capable of efficiently deriving potential based constitutive models, in analytical form is presented. This package, running under MACSYMA, has the following features: partial differentiation, tensor computations, automatic grouping and labeling of common factors, expression substitution and simplification, back substitution of invariant and tensorial relations and a relational data base. Also limited aspects of invariant theory were incorporated into SDICE due to the utilization of potentials as a starting point and the desire for these potentials to be frame invariant (objective). Finally not only calculation of flow and/or evolutionary laws were accomplished but also the determination of history independent nonphysical coefficients in terms of physically measurable parameters, e.g., Young's modulus, was achieved. The uniqueness of SDICE resides in its ability to manipulate expressions in a general yet predefined order and simplify expressions so as to limit expression growth. Results are displayed when applicable utilizing index notation.

Low-frequency analogue Hawking radiation: The Korteweg-de Vries model

NASA Astrophysics Data System (ADS)

Coutant, Antonin; Weinfurtner, Silke

2018-01-01

We derive analytic expressions for the low-frequency properties of the analogue Hawking radiation in a general weak-dispersive medium. A thermal low-frequency part of the spectrum is expected even when dispersive effects become significant. We consider the two most common class of weak-dispersive media and investigate all possible anomalous scattering processes due inhomogeneous background flows. We first argue that under minimal assumptions, the scattering processes in near-critical flows are well described by a linearized Korteweg-de Vries equation. Within our theoretical model grey-body factors are neglected, that is, the mode comoving with the flow decouples from the other ones. We also exhibit a flow example with an exact expression for the effective temperature. We see that this temperature coincides with the Hawking one only when the dispersive length scale is much smaller than the flow gradient scale. We apply the same method in inhomogeneous flows without an analogue horizon. In this case, the spectrum coefficients decrease with decreasing frequencies. Our findings are in agreement with previous numerical works, generalizing their findings to arbitrary flow profiles. Our analytical expressions provide estimates to guide ongoing experimental efforts.

Modal expansions in periodic photonic systems with material loss and dispersion

NASA Astrophysics Data System (ADS)

Wolff, Christian; Busch, Kurt; Mortensen, N. Asger

2018-03-01

We study band-structure properties of periodic optical systems composed of lossy and intrinsically dispersive materials. To this end, we develop an analytical framework based on adjoint modes of a lossy periodic electromagnetic system and show how the problem of linearly dependent eigenmodes in the presence of material dispersion can be overcome. We then formulate expressions for the band-structure derivative (∂ ω )/(∂ k ) (complex group velocity) and the local and total density of transverse optical states. Our exact expressions hold for 3D periodic arrays of materials with arbitrary dispersion properties and in general need to be evaluated numerically. They can be generalized to systems with two, one, or no directions of periodicity provided the fields are localized along nonperiodic directions. Possible applications are photonic crystals, metamaterials, metasurfaces composed of highly dispersive materials such as metals or lossless photonic crystals, and metamaterials or metasurfaces strongly coupled to resonant perturbations such as quantum dots or excitons in 2D materials. For illustration purposes, we analytically evaluate our expressions for some simple systems consisting of lossless dielectrics with one sharp Lorentzian material resonance added. By combining several Lorentz poles, this provides an avenue to perturbatively treat quite general material loss bands in photonic crystals.

Analytical study of the reflection and transmission coefficient of the submarine interface

NASA Astrophysics Data System (ADS)

Zhang, Guangli; Hao, Chongtao; Yao, Chen

2018-05-01

The analytical study of the reflection and transmission coefficient of the seafloor interface is essential for the characterization of the ocean bottom in marine seismic exploration. Based on the boundary conditions of the seafloor interface, the analytical expression of the reflection and transmission coefficient at the submarine interface is derived in this study by using the steady-state wave solution of the elastic wave in a homogeneous, isotropic medium. With this analytical expression, the characteristics of the reflection and transmission coefficient at the submarine interface are analysed and discussed using critical angles. The results show that the change in the reflection and transmission coefficient with the incidence angle presents a "segmented" characteristic, in which the critical angle is the dividing point. The amplitude value and phase angle of the coefficient at the submarine interface change dramatically at the critical angle, which is related to the P- and S-wave velocities in the seabed layer. Compared with the stiff seabed, the soft seabed has a larger P-wave critical angle and an absence of the converted S-wave critical angle, owing to the low P- and S-wave velocities in the solid seabed layer. By analysing and discussing the special changes that occur in the coefficient values at the critical angle, the reflection and transmission characteristics of the different incident angles are obtained. Synthetic models of both stiff and soft seafloors are provided in this study to verify the analytical results. Finally, we compared our synthetic results with real data from the Gulf of Mexico, which enabled the validation of our conclusions.

Spherical Pendulum Small Oscillations for Slewing Crane Motion

PubMed Central

Perig, Alexander V.; Stadnik, Alexander N.; Deriglazov, Alexander I.

2014-01-01

The present paper focuses on the Lagrange mechanics-based description of small oscillations of a spherical pendulum with a uniformly rotating suspension center. The analytical solution of the natural frequencies' problem has been derived for the case of uniform rotation of a crane boom. The payload paths have been found in the inertial reference frame fixed on earth and in the noninertial reference frame, which is connected with the rotating crane boom. The numerical amplitude-frequency characteristics of the relative payload motion have been found. The mechanical interpretation of the terms in Lagrange equations has been outlined. The analytical expression and numerical estimation for cable tension force have been proposed. The numerical computational results, which correlate very accurately with the experimental observations, have been shown. PMID:24526891

Exact solution for the Poisson field in a semi-infinite strip.

PubMed

Cohen, Yossi; Rothman, Daniel H

2017-04-01

The Poisson equation is associated with many physical processes. Yet exact analytic solutions for the two-dimensional Poisson field are scarce. Here we derive an analytic solution for the Poisson equation with constant forcing in a semi-infinite strip. We provide a method that can be used to solve the field in other intricate geometries. We show that the Poisson flux reveals an inverse square-root singularity at a tip of a slit, and identify a characteristic length scale in which a small perturbation, in a form of a new slit, is screened by the field. We suggest that this length scale expresses itself as a characteristic spacing between tips in real Poisson networks that grow in response to fluxes at tips.

Interacting steps with finite-range interactions: Analytical approximation and numerical results

NASA Astrophysics Data System (ADS)

Jaramillo, Diego Felipe; Téllez, Gabriel; González, Diego Luis; Einstein, T. L.

2013-05-01

We calculate an analytical expression for the terrace-width distribution P(s) for an interacting step system with nearest- and next-nearest-neighbor interactions. Our model is derived by mapping the step system onto a statistically equivalent one-dimensional system of classical particles. The validity of the model is tested with several numerical simulations and experimental results. We explore the effect of the range of interactions q on the functional form of the terrace-width distribution and pair correlation functions. For physically plausible interactions, we find modest changes when next-nearest neighbor interactions are included and generally negligible changes when more distant interactions are allowed. We discuss methods for extracting from simulated experimental data the characteristic scale-setting terms in assumed potential forms.

Phase-locking of combination-cylinder discharge CO2 laser

NASA Astrophysics Data System (ADS)

Xu, Yonggen

2014-05-01

A new type of laser resonator is presented to obtained good coherent beam and the parameters of the laser beam are calculated. The principle of phase-locking is described based on the injection-locking, the properties of the injected beam in the resonator are studied in detail. The output beam from output mirror is an annular laser beam with zero central intensity. An analytical expression for the annular laser beam through the ABCD optical system is derived. Typical numerical examples are calculated to confirm our analytical results. It is shown that the good coherent beam can be obtained through phase-locking, and the central intensity of annular beam through ABCD optical system will become maximum when the parameters of laser beam are selected reasonably.

Combinatoric analysis of heterogeneous stochastic self-assembly.

PubMed

D'Orsogna, Maria R; Zhao, Bingyu; Berenji, Bijan; Chou, Tom

2013-09-28

We analyze a fully stochastic model of heterogeneous nucleation and self-assembly in a closed system with a fixed total particle number M, and a fixed number of seeds Ns. Each seed can bind a maximum of N particles. A discrete master equation for the probability distribution of the cluster sizes is derived and the corresponding cluster concentrations are found using kinetic Monte-Carlo simulations in terms of the density of seeds, the total mass, and the maximum cluster size. In the limit of slow detachment, we also find new analytic expressions and recursion relations for the cluster densities at intermediate times and at equilibrium. Our analytic and numerical findings are compared with those obtained from classical mass-action equations and the discrepancies between the two approaches analyzed.

Investigation of Cepstrum Analysis for Seismic/Acoustic Signal Sensor Range Determination.

DTIC Science & Technology

1981-01-01

distorted by transmission through a linear system . For example, the effect of multipath and reverberation may be modeled in terms of a signal that is...called the short time averaged cepstrum. To derive some analytical expressions for short time average cepstrums we choose some functions of interest...linear process applied to the time series or any equivalent time function Repiod Period The amount of time required for one cycle of a time series Saphe

Saturation of the Hosing Instability in Quasilinear Plasma Accelerators

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lehe, R.; Schroeder, C. B.; Vay, J. -L.

The beam hosing instability is analyzed theoretically for a witness beam in the quasilinear regime of plasma accelerators. In this regime, the hosing instability saturates, even for a monoenergetic bunch, at a level much less than standard scalings predict. Analytic expressions are derived for the saturation distance and amplitude and are in agreement with numerical results. Saturation is due to the natural head-to-tail variations in the focusing force, including the self-consistent transverse beam loading.

Saturation of the Hosing Instability in Quasilinear Plasma Accelerators

DOE PAGES

Lehe, R.; Schroeder, C. B.; Vay, J. -L.; ...

2017-12-13

The beam hosing instability is analyzed theoretically for a witness beam in the quasilinear regime of plasma accelerators. In this regime, the hosing instability saturates, even for a monoenergetic bunch, at a level much less than standard scalings predict. Analytic expressions are derived for the saturation distance and amplitude and are in agreement with numerical results. Saturation is due to the natural head-to-tail variations in the focusing force, including the self-consistent transverse beam loading.

Multiscale modeling of fluid transport in tumors.

PubMed

Chapman, S Jonathan; Shipley, Rebecca J; Jawad, Rossa

2008-11-01

A model for fluid flow through the leaky neovasculature and porous interstitium of a solid tumor is developed. A network of isolated capillaries is analyzed in the limit of small capillary radius, and analytical expressions for the hydraulic conductivities and fractional leakage coefficients derived. This model is then homogenized to give a continuum description in terms of the vascular density. The resulting equations comprise a double porous medium with coupled Darcy flow through the interstitium and vasculature.

Spectral analysis for weighted tree-like fractals

NASA Astrophysics Data System (ADS)

Dai, Meifeng; Chen, Yufei; Wang, Xiaoqian; Sun, Yu; Su, Weiyi

2018-02-01

Much information about the structural properties and dynamical aspects of a network is measured by the eigenvalues of its normalized Laplacian matrix. In this paper, we aim to present a study on the spectra of the normalized Laplacian of weighted tree-like fractals. We analytically obtain the relationship between the eigenvalues and their multiplicities for two successive generations. As an example of application of these results, we then derive closed-form expressions for their multiplicative Kirchhoff index and Kemeny's constant.

Dispersion of TE modes in slab waveguides with reference to double heterostructure semiconductor lasers

NASA Astrophysics Data System (ADS)

Buus, J.

1980-06-01

The group index for TE modes in an asymmetrical slab waveguide is investigated, and a simple analytical expression is derived. It is shown that the product of the phase and group indices is related to the power fraction in each of the three layers of the waveguide. The results are of interest in the analysis of double heterostructure semiconductor lasers. Theoretical and experimental results for lasers emitting at 1.55 microns are compared.

Eddy currents in a conducting sphere

NASA Technical Reports Server (NTRS)

Bergman, John; Hestenes, David

1986-01-01

This report analyzes the eddy current induced in a solid conducting sphere by a sinusoidal current in a circular loop. Analytical expressions for the eddy currents are derived as a power series in the vectorial displacement of the center of the sphere from the axis of the loop. These are used for first order calculations of the power dissipated in the sphere and the force and torque exerted on the sphere by the electromagnetic field of the loop.

A generalized theory for the design of contraction cones and other low speed ducts

NASA Technical Reports Server (NTRS)

Barger, R. L.; Bowen, J. T.

1972-01-01

A generalization of the Tsien method of contraction cone design is described. The design velocity distribution is expressed in such a form that the required high order derivatives can be obtained by recursion rather than by numerical or analytic differentiation. The method is applicable to the design of diffusers and converging-diverging ducts as well as contraction cones. The computer program is described and a FORTRAN listing of the program is provided.

Accurate monoenergetic electron parameters of laser wakefield in a bubble model

NASA Astrophysics Data System (ADS)

Raheli, A.; Rahmatallahpur, S. H.

2012-11-01

A reliable analytical expression for the potential of plasma waves with phase velocities near the speed of light is derived. The presented spheroid cavity model is more consistent than the previous spherical and ellipsoidal model and it explains the mono-energetic electron trajectory more accurately, especially at the relativistic region. As a result, the quasi-mono-energetic electrons output beam interacting with the laser plasma can be more appropriately described with this model.

A scaling law for the critical current of Nb3Sn stands based on strong-coupling theory of superconductivity

NASA Astrophysics Data System (ADS)

Oh, Sangjun; Kim, Keeman

2006-02-01

We study the transition temperature Tc, the thermodynamic critical field Bc, and the upper critical field Bc2 of Nb3Sn with Eliashberg theory of strongly coupled superconductors using the Einstein spectrum α2(ω)F(ω)=λ<ω2>1/2δ(ω-<ω2>1/2). The strain dependences of λ(ɛ) and <ω2>1/2(V) are introduced from the empirical strain dependence of Tc(V) for three model cases. It is found that the empirical relation Tc(V)/Tc(0)=[Bc2(4.2 K,V)/Bc2(4.2 K,0)]1/w (w~3) is mainly due to the low-energy-phonon mode softening. We derive analytic expressions for the strain and temperature dependences of Bc(T,V) and Bc2(T,V) and the Ginzburg-Landau parameter κ(T,V) from the numerical calculation results. The Summers refinement on the temperature dependence of κ(T) shows deviation from our calculation results. We propose a unified scaling law of flux pinning in Nb3Sn strands in the form of the Kramer model with the analytic expressions of Bc2(T,V) and κ(T,V) derived in this work. It is shown that the proposed scaling law gives a reasonable fit to the reported data with only eight fitting parameters.

Linear and Order Statistics Combiners for Pattern Classification

NASA Technical Reports Server (NTRS)

Tumer, Kagan; Ghosh, Joydeep; Lau, Sonie (Technical Monitor)

2001-01-01

Several researchers have experimentally shown that substantial improvements can be obtained in difficult pattern recognition problems by combining or integrating the outputs of multiple classifiers. This chapter provides an analytical framework to quantify the improvements in classification results due to combining. The results apply to both linear combiners and order statistics combiners. We first show that to a first order approximation, the error rate obtained over and above the Bayes error rate, is directly proportional to the variance of the actual decision boundaries around the Bayes optimum boundary. Combining classifiers in output space reduces this variance, and hence reduces the 'added' error. If N unbiased classifiers are combined by simple averaging. the added error rate can be reduced by a factor of N if the individual errors in approximating the decision boundaries are uncorrelated. Expressions are then derived for linear combiners which are biased or correlated, and the effect of output correlations on ensemble performance is quantified. For order statistics based non-linear combiners, we derive expressions that indicate how much the median, the maximum and in general the i-th order statistic can improve classifier performance. The analysis presented here facilitates the understanding of the relationships among error rates, classifier boundary distributions, and combining in output space. Experimental results on several public domain data sets are provided to illustrate the benefits of combining and to support the analytical results.

Analytical results for the time-dependent current density distribution of expanding ultracold gases after a sudden change of the confining potential

NASA Astrophysics Data System (ADS)

Boumaza, R.; Bencheikh, K.

2017-12-01

Using the so-called operator product expansion to lowest order, we extend the work in Campbell et al (2015 Phys. Rev. Lett 114 125302) by deriving a simple analytical expression for the long-time asymptotic one-body reduced density matrix during free expansion for a one-dimensional system of bosons with large atom number interacting through a repulsive delta potential initially confined by a potential well. This density matrix allows direct access to the momentum distribution and also to the mass current density. For initially confining power-law potentials we give explicit expressions, in the limits of very weak and very strong interaction, for the current density distributions during the free expansion. In the second part of the work we consider the expansion of ultracold gas from a confining harmonic trap to another harmonic trap with a different frequency. For the case of a quantum impenetrable gas of bosons (a Tonks-Girardeau gas) with a given atom number, we present an exact analytical expression for the mass current distribution (mass transport) after release from one harmonic trap to another harmonic trap. It is shown that, for a harmonically quenched Tonks-Girardeau gas, the current distribution is a suitable collective observable and under the weak quench regime, it exhibits oscillations at the same frequencies as those recently predicted for the peak momentum distribution in the breathing mode. The analysis is extended to other possible quenched systems.

Analytic gradient for second order Møller-Plesset perturbation theory with the polarizable continuum model based on the fragment molecular orbital method.

PubMed

Nagata, Takeshi; Fedorov, Dmitri G; Li, Hui; Kitaura, Kazuo

2012-05-28

A new energy expression is proposed for the fragment molecular orbital method interfaced with the polarizable continuum model (FMO/PCM). The solvation free energy is shown to be more accurate on a set of representative polypeptides with neutral and charged residues, in comparison to the original formulation at the same level of the many-body expansion of the electrostatic potential determining the apparent surface charges. The analytic first derivative of the energy with respect to nuclear coordinates is formulated at the second-order Møller-Plesset (MP2) perturbation theory level combined with PCM, for which we derived coupled perturbed Hartree-Fock equations. The accuracy of the analytic gradient is demonstrated on test calculations in comparison to numeric gradient. Geometry optimization of the small Trp-cage protein (PDB: 1L2Y) is performed with FMO/PCM/6-31(+)G(d) at the MP2 and restricted Hartree-Fock with empirical dispersion (RHF/D). The root mean square deviations between the FMO optimized and NMR experimental structure are found to be 0.414 and 0.426 Å for RHF/D and MP2, respectively. The details of the hydrogen bond network in the Trp-cage protein are revealed.

Analytic gradient for second order Møller-Plesset perturbation theory with the polarizable continuum model based on the fragment molecular orbital method

NASA Astrophysics Data System (ADS)

Nagata, Takeshi; Fedorov, Dmitri G.; Li, Hui; Kitaura, Kazuo

2012-05-01

A new energy expression is proposed for the fragment molecular orbital method interfaced with the polarizable continuum model (FMO/PCM). The solvation free energy is shown to be more accurate on a set of representative polypeptides with neutral and charged residues, in comparison to the original formulation at the same level of the many-body expansion of the electrostatic potential determining the apparent surface charges. The analytic first derivative of the energy with respect to nuclear coordinates is formulated at the second-order Møller-Plesset (MP2) perturbation theory level combined with PCM, for which we derived coupled perturbed Hartree-Fock equations. The accuracy of the analytic gradient is demonstrated on test calculations in comparison to numeric gradient. Geometry optimization of the small Trp-cage protein (PDB: 1L2Y) is performed with FMO/PCM/6-31(+)G(d) at the MP2 and restricted Hartree-Fock with empirical dispersion (RHF/D). The root mean square deviations between the FMO optimized and NMR experimental structure are found to be 0.414 and 0.426 Å for RHF/D and MP2, respectively. The details of the hydrogen bond network in the Trp-cage protein are revealed.

How classical gluon fields generate odd azimuthal harmonics for the two-gluon correlation function in high-energy collisions

NASA Astrophysics Data System (ADS)

Kovchegov, Yuri V.; Skokov, Vladimir V.

2018-05-01

We show that, in the saturation/color glass condensate framework, odd azimuthal harmonics of the two-gluon correlation function with a long-range separation in rapidity are generated by the higher-order saturation corrections in the interactions with the projectile and the target. At the very least, the odd harmonics require three scatterings in the projectile and three scatterings in the target. We derive the leading-order expression for the two-gluon production cross section which generates odd harmonics: the expression includes all-order interactions with the target and three interactions with the projectile. We evaluate the obtained expression both analytically and numerically, confirming that the odd-harmonics contribution to the two-gluon production in the saturation framework is nonzero.

Complementarity of SOMAscan to LC-MS/MS and RNA-seq for quantitative profiling of human embryonic and mesenchymal stem cells.

PubMed

Billing, Anja M; Ben Hamidane, Hisham; Bhagwat, Aditya M; Cotton, Richard J; Dib, Shaima S; Kumar, Pankaj; Hayat, Shahina; Goswami, Neha; Suhre, Karsten; Rafii, Arash; Graumann, Johannes

2017-01-06

Dynamic range limitations are challenging to proteomics, particularly in clinical samples. Affinity proteomics partially overcomes this, yet suffers from dependence on reagent quality. SOMAscan, an aptamer-based platform for over 1000 proteins, avoids that issue using nucleic acid binders. Targets include low expressed proteins not easily accessible by other approaches. Here we report on the potential of SOMAscan for the study of differently sourced mesenchymal stem cells (MSC) in comparison to LC-MS/MS and RNA sequencing. While targeting fewer analytes, SOMAscan displays high precision and dynamic range coverage, allowing quantification of proteins not measured by the other platforms. Expression between cell types (ESC and MSC) was compared across techniques and uncovered the expected large differences. Sourcing was investigated by comparing subtypes: bone marrow-derived, standard in clinical studies, and ESC-derived MSC, thought to hold similar potential but devoid of inter-donor variability and proliferating faster in vitro. We confirmed subtype-equivalency, as well as vesicle and extracellular matrix related processes in MSC. In contrast, the proliferative nature of ESC was captured less by SOMAscan, where nuclear proteins are underrepresented. The complementary of SOMAscan allowed the comprehensive exploration of CD markers and signaling molecules, not readily accessible otherwise and offering unprecedented potential in subtype characterization. Mesenchymal stem cells (MSC) represent promising stem cell-derived therapeutics as indicated by their application in >500 clinical trials currently registered with the NIH. Tissue-derived MSC require invasive harvesting and imply donor-to-donor differences, to which embryonic stem cell (ESC)-derived MSC may provide an alternative and thus warrant thorough characterization. In continuation of our previous study where we compared in depth embryonic stem cells (ESC) and MSC from two sources (bone marrow and ESC-derived), we included the aptamer-based SOMAscan assay, complementing LC-MS/MS and RNA-seq data. Furthermore, SOMAscan, a targeted proteomics platform developed for analyzing clinical samples, has been benchmarked against established analytical platforms (LC-MS/MS and RNA-seq) using stem cell comparisons as a model. Copyright © 2016 The Authors. Published by Elsevier B.V. All rights reserved.

Pattern Storage, Bifurcations, and Groupwise Correlation Structure of an Exactly Solvable Asymmetric Neural Network Model.

PubMed

Fasoli, Diego; Cattani, Anna; Panzeri, Stefano

2018-05-01

Despite their biological plausibility, neural network models with asymmetric weights are rarely solved analytically, and closed-form solutions are available only in some limiting cases or in some mean-field approximations. We found exact analytical solutions of an asymmetric spin model of neural networks with arbitrary size without resorting to any approximation, and we comprehensively studied its dynamical and statistical properties. The network had discrete time evolution equations and binary firing rates, and it could be driven by noise with any distribution. We found analytical expressions of the conditional and stationary joint probability distributions of the membrane potentials and the firing rates. By manipulating the conditional probability distribution of the firing rates, we extend to stochastic networks the associating learning rule previously introduced by Personnaz and coworkers. The new learning rule allowed the safe storage, under the presence of noise, of point and cyclic attractors, with useful implications for content-addressable memories. Furthermore, we studied the bifurcation structure of the network dynamics in the zero-noise limit. We analytically derived examples of the codimension 1 and codimension 2 bifurcation diagrams of the network, which describe how the neuronal dynamics changes with the external stimuli. This showed that the network may undergo transitions among multistable regimes, oscillatory behavior elicited by asymmetric synaptic connections, and various forms of spontaneous symmetry breaking. We also calculated analytically groupwise correlations of neural activity in the network in the stationary regime. This revealed neuronal regimes where, statistically, the membrane potentials and the firing rates are either synchronous or asynchronous. Our results are valid for networks with any number of neurons, although our equations can be realistically solved only for small networks. For completeness, we also derived the network equations in the thermodynamic limit of infinite network size and we analytically studied their local bifurcations. All the analytical results were extensively validated by numerical simulations.

The effect of newly synthesized progesterone derivatives on apoptotic and angiogenic pathway in MCF-7 breast cancer cells.

PubMed

Yahya, Shaymaa M M; Abdelhamid, Abdou O; Abd-Elhalim, Mervat M; Elsayed, Ghada H; Eskander, Emad F

2017-10-01

Due to its high potency and selectivity, anticancer agents consisting of combined molecules have gained great interests. The current study introduces newly synthesized progesterone derivatives of promising anticancer effect. Moreover, the pro-apoptotic and anti-angiogenic effects of these compounds were studied extensively. Several thiazole, pyridine, pyrazole, thiazolopyridine and pyrazolopyridine progesterone derivatives were synthesized. The structure of the novel progesterone derivatives was elucidated and confirmed using the analytical and spectral data. This novel derivatives were tested for their cytotoxic effect against human breast cancer cells (MCF-7) using neutral red uptake assay. Tested compounds showed anticancer activity against MCF-7 cancer cell line in the descending order of 7>2>3>8>6>9>4. The expression levels of Bcl-2, survivin, CCND1, CDC2, P53 and P21, VEGF, Hif-1α, MMP-2, MMP-9, Ang-1, Ang-2, and FGF-1 genes were investigated using QRT-PCR (Quantitative real time-polymerase chain reaction). The study clarified that compounds 2, 3, 4, 6, 7, 8 and 9 showed significant pro-apoptotic effect through the down regulation of Bcl-2., besides, survivin and CCND1 expression levels were down regulated by compounds 3, 4, 6, 7, 8, 9. However, Compound 4 may exert this pro-apoptotic effect through the up-regulation of P53 gene expression. On the other hand, the anti-angiogenic effect of these newly synthesized derivatives was due to their down regulation of VEGF, Ang-2, MMP-9 and FGF-1; and the up-regulation of HIF-1α and ang-1. This study recommended promising pro-apoptotic and anti-angiogenic anticancer agents acting through the regulation of key regulators of apoptosis, cell cycle genes, and pro-angiogenic genes. Copyright © 2017 Elsevier Inc. All rights reserved.

On the specta of X-ray bursters: Expansion and contraction stages

NASA Technical Reports Server (NTRS)

Titarchuk, Lev

1994-01-01

The theory of spectral formation during the explosion and contraction stages of X-ray bursters, which include the effects of Computonization and free-free absorption and emission, is described. Analytical expressions are provided for color ratios, and the spectral shape is given as a function of input parameters, elemental abundance, neutron star mass and radius, and Eddington ratio. An Eulerian calculation is used to determine the photospheric evolution accurately during the Eddington luminosity phase. The developed analytical theory for hydrodynamics of the expansion takes into account the dependence of Compton scattering opacity on electron temperature. An analytical expression is derived from the sonic point position and the value of the sonic velcoity. Using this value as a boundary condition at the sonic point, the velocity, density, and temperature profile are calculated throughout the whole photosphere. It is shown that the atmopsphere radiates spectra having a low-energy power-law shape and blackbody-like hard tail. In the expansion stage the spectra depend strongly on the temperature of the helium-burning zone at the neutron star surface. The X-ray photosheric radius increases to approximately 100 km or more, depending on the condition of the nuclear burning on the surface of the neutron star in the course of the expansion.

Optical soliton solutions, periodic wave solutions and complexitons of the cubic Schrödinger equation with a bounded potential

NASA Astrophysics Data System (ADS)

Yan, Xue-Wei; Tian, Shou-Fu; Dong, Min-Jie; Zou, Li

2018-01-01

In this paper, we consider the cubic Schrödinger equation with a bounded potential, which describes the propagation properties of optical soliton solutions. By employing an ansatz method, we precisely derive the bright and dark soliton solutions of the equation. Moreover, we obtain three classes of analytic periodic wave solutions expressed in terms of the Jacobi's elliptic functions including cn ,sn and dn functions. Finally, by using a tanh function method, its complexitons solutions are derived in a very natural way. It is hoped that our results can enrich the nonlinear dynamical behaviors of the cubic Schrödinger equation with a bounded potential.

Direct S -matrix calculation for diffractive structures and metasurfaces

NASA Astrophysics Data System (ADS)

Shcherbakov, Alexey A.; Stebunov, Yury V.; Baidin, Denis F.; Kämpfe, Thomas; Jourlin, Yves

2018-06-01

The paper presents a derivation of analytical components of S matrices for arbitrary planar diffractive structures and metasurfaces in the Fourier domain. The attained general formulas for S -matrix components can be applied within both formulations in the Cartesian and curvilinear metric. A numerical method based on these results can benefit from all previous improvements of the Fourier domain methods. In addition, we provide expressions for S -matrix calculation in the case of periodically corrugated layers of two-dimensional materials, which are valid for arbitrary corrugation depth-to-period ratios. As an example, the derived equations are used to simulate resonant grating excitation of graphene plasmons and the impact of a silica interlayer on corresponding reflection curves.

Laboratory study of low-β forces in arched, line-tied magnetic flux ropes

NASA Astrophysics Data System (ADS)

Myers, C. E.; Yamada, M.; Ji, H.; Yoo, J.; Jara-Almonte, J.; Fox, W.

2016-11-01

The loss-of-equilibrium is a solar eruption mechanism whereby a sudden breakdown of the magnetohydrodynamic force balance in the Sun's corona ejects a massive burst of particles and energy into the heliosphere. Predicting a loss-of-equilibrium, which has more recently been formulated as the torus instability, relies on a detailed understanding of the various forces that hold the pre-eruption magnetic flux rope in equilibrium. Traditionally, idealized analytical force expressions are used to derive simplified eruption criteria that can be compared to solar observations and modeling. What is missing, however, is a validation that these idealized analytical force expressions can be applied to the line-tied, low-aspect-ratio conditions of the corona. In this paper, we address this shortcoming by using a laboratory experiment to study the forces that act on long-lived, arched, line-tied magnetic flux ropes. Three key force terms are evaluated over a wide range of experimental conditions: (1) the upward hoop force; (2) the downward strapping force; and (3) the downward toroidal field tension force. First, the laboratory force measurements show that, on average, the three aforementioned force terms cancel to produce a balanced line-tied equilibrium. This finding validates the laboratory force measurement techniques developed here, which were recently used to identify a dynamic toroidal field tension force that can prevent flux rope eruptions [Myers et al., Nature 528, 526 (2015)]. The verification of magnetic force balance also confirms the low-β assumption that the plasma thermal pressure is negligible in these experiments. Next, the measured force terms are directly compared to corresponding analytical expressions. While the measured and analytical forces are found to be well correlated, the low-aspect-ratio, line-tied conditions in the experiment are found to both reduce the measured hoop force and increase the measured tension force with respect to analytical expectations. These two co-directed effects combine to generate laboratory flux rope equilibria at lower altitudes than are predicted analytically. Such considerations are expected to modify the loss-of-equilibrium eruption criteria for analogous flux ropes in the solar corona.

Laboratory study of low- β forces in arched, line-tied magnetic flux ropes

DOE PAGES

Myers, C. E.; Yamada, M.; Ji, H.; ...

2016-11-04

Here, the loss-of-equilibrium is a solar eruption mechanism whereby a sudden breakdown of the magnetohydrodynamic force balance in the Sun's corona ejects a massive burst of particles and energy into the heliosphere. Predicting a loss-of-equilibrium, which has more recently been formulated as the torus instability, relies on a detailed understanding of the various forces that hold the pre-eruption magnetic flux rope in equilibrium. Traditionally, idealized analytical force expressions are used to derive simplified eruption criteria that can be compared to solar observations and modeling. What is missing, however, is a validation that these idealized analytical force expressions can be appliedmore » to the line-tied, low-aspect-ratio conditions of the corona. In this paper, we address this shortcoming by using a laboratory experiment to study the forces that act on long-lived, arched, line-tied magnetic flux ropes. Three key force terms are evaluated over a wide range of experimental conditions: (1) the upward hoop force; (2) the downward strapping force; and (3) the downward toroidal field tension force. First, the laboratory force measurements show that, on average, the three aforementioned force terms cancel to produce a balanced line-tied equilibrium. This finding validates the laboratory force measurement techniques developed here, which were recently used to identify a dynamic toroidal field tension force that can prevent flux rope eruption. The verification of magnetic force balance also confirms the low-beta assumption that the plasma thermal pressure is negligible in these experiments. Next, the measured force terms are directly compared to corresponding analytical expressions. While the measured and analytical forces are found to be well correlated, the low-aspect-ratio, line-tied conditions in the experiment are found to both reduce the measured hoop force and increase the measured tension force with respect to analytical expectations. These two co-directed effects combine to generate laboratory flux rope equilibria at lower altitudes than are predicted analytically. Such considerations are expected to modify the loss-of-equilibrium eruption criteria for analogous flux ropes in the solar corona.« less

Laboratory study of low- β forces in arched, line-tied magnetic flux ropes

DOE Office of Scientific and Technical Information (OSTI.GOV)

Myers, C. E.; Yamada, M.; Ji, H.

Here, the loss-of-equilibrium is a solar eruption mechanism whereby a sudden breakdown of the magnetohydrodynamic force balance in the Sun's corona ejects a massive burst of particles and energy into the heliosphere. Predicting a loss-of-equilibrium, which has more recently been formulated as the torus instability, relies on a detailed understanding of the various forces that hold the pre-eruption magnetic flux rope in equilibrium. Traditionally, idealized analytical force expressions are used to derive simplified eruption criteria that can be compared to solar observations and modeling. What is missing, however, is a validation that these idealized analytical force expressions can be appliedmore » to the line-tied, low-aspect-ratio conditions of the corona. In this paper, we address this shortcoming by using a laboratory experiment to study the forces that act on long-lived, arched, line-tied magnetic flux ropes. Three key force terms are evaluated over a wide range of experimental conditions: (1) the upward hoop force; (2) the downward strapping force; and (3) the downward toroidal field tension force. First, the laboratory force measurements show that, on average, the three aforementioned force terms cancel to produce a balanced line-tied equilibrium. This finding validates the laboratory force measurement techniques developed here, which were recently used to identify a dynamic toroidal field tension force that can prevent flux rope eruption. The verification of magnetic force balance also confirms the low-beta assumption that the plasma thermal pressure is negligible in these experiments. Next, the measured force terms are directly compared to corresponding analytical expressions. While the measured and analytical forces are found to be well correlated, the low-aspect-ratio, line-tied conditions in the experiment are found to both reduce the measured hoop force and increase the measured tension force with respect to analytical expectations. These two co-directed effects combine to generate laboratory flux rope equilibria at lower altitudes than are predicted analytically. Such considerations are expected to modify the loss-of-equilibrium eruption criteria for analogous flux ropes in the solar corona.« less

Thinking outside the box: effects of modes larger than the survey on matter power spectrum covariance

DOE Office of Scientific and Technical Information (OSTI.GOV)

Putter, Roland de; Wagner, Christian; Verde, Licia

2012-04-01

Accurate power spectrum (or correlation function) covariance matrices are a crucial requirement for cosmological parameter estimation from large scale structure surveys. In order to minimize reliance on computationally expensive mock catalogs, it is important to have a solid analytic understanding of the different components that make up a covariance matrix. Considering the matter power spectrum covariance matrix, it has recently been found that there is a potentially dominant effect on mildly non-linear scales due to power in modes of size equal to and larger than the survey volume. This beat coupling effect has been derived analytically in perturbation theory andmore » while it has been tested with simulations, some questions remain unanswered. Moreover, there is an additional effect of these large modes, which has so far not been included in analytic studies, namely the effect on the estimated average density which enters the power spectrum estimate. In this article, we work out analytic, perturbation theory based expressions including both the beat coupling and this local average effect and we show that while, when isolated, beat coupling indeed causes large excess covariance in agreement with the literature, in a realistic scenario this is compensated almost entirely by the local average effect, leaving only ∼ 10% of the excess. We test our analytic expressions by comparison to a suite of large N-body simulations, using both full simulation boxes and subboxes thereof to study cases without beat coupling, with beat coupling and with both beat coupling and the local average effect. For the variances, we find excellent agreement with the analytic expressions for k < 0.2 hMpc{sup −1} at z = 0.5, while the correlation coefficients agree to beyond k = 0.4 hMpc{sup −1}. As expected, the range of agreement increases towards higher redshift and decreases slightly towards z = 0. We finish by including the large-mode effects in a full covariance matrix description for arbitrary survey geometry and confirming its validity using simulations. This may be useful as a stepping stone towards building an actual galaxy (or other tracer's) power spectrum covariance matrix.« less

A unifying theory for genetic epidemiological analysis of binary disease data

PubMed Central

2014-01-01

Background Genetic selection for host resistance offers a desirable complement to chemical treatment to control infectious disease in livestock. Quantitative genetics disease data frequently originate from field studies and are often binary. However, current methods to analyse binary disease data fail to take infection dynamics into account. Moreover, genetic analyses tend to focus on host susceptibility, ignoring potential variation in infectiousness, i.e. the ability of a host to transmit the infection. This stands in contrast to epidemiological studies, which reveal that variation in infectiousness plays an important role in the progression and severity of epidemics. In this study, we aim at filling this gap by deriving an expression for the probability of becoming infected that incorporates infection dynamics and is an explicit function of both host susceptibility and infectiousness. We then validate this expression according to epidemiological theory and by simulating epidemiological scenarios, and explore implications of integrating this expression into genetic analyses. Results Our simulations show that the derived expression is valid for a range of stochastic genetic-epidemiological scenarios. In the particular case of variation in susceptibility only, the expression can be incorporated into conventional quantitative genetic analyses using a complementary log-log link function (rather than probit or logit). Similarly, if there is moderate variation in both susceptibility and infectiousness, it is possible to use a logarithmic link function, combined with an indirect genetic effects model. However, in the presence of highly infectious individuals, i.e. super-spreaders, the use of any model that is linear in susceptibility and infectiousness causes biased estimates. Thus, in order to identify super-spreaders, novel analytical methods using our derived expression are required. Conclusions We have derived a genetic-epidemiological function for quantitative genetic analyses of binary infectious disease data, which, unlike current approaches, takes infection dynamics into account and allows for variation in host susceptibility and infectiousness. PMID:24552188

A unifying theory for genetic epidemiological analysis of binary disease data.

PubMed

Lipschutz-Powell, Debby; Woolliams, John A; Doeschl-Wilson, Andrea B

2014-02-19

Genetic selection for host resistance offers a desirable complement to chemical treatment to control infectious disease in livestock. Quantitative genetics disease data frequently originate from field studies and are often binary. However, current methods to analyse binary disease data fail to take infection dynamics into account. Moreover, genetic analyses tend to focus on host susceptibility, ignoring potential variation in infectiousness, i.e. the ability of a host to transmit the infection. This stands in contrast to epidemiological studies, which reveal that variation in infectiousness plays an important role in the progression and severity of epidemics. In this study, we aim at filling this gap by deriving an expression for the probability of becoming infected that incorporates infection dynamics and is an explicit function of both host susceptibility and infectiousness. We then validate this expression according to epidemiological theory and by simulating epidemiological scenarios, and explore implications of integrating this expression into genetic analyses. Our simulations show that the derived expression is valid for a range of stochastic genetic-epidemiological scenarios. In the particular case of variation in susceptibility only, the expression can be incorporated into conventional quantitative genetic analyses using a complementary log-log link function (rather than probit or logit). Similarly, if there is moderate variation in both susceptibility and infectiousness, it is possible to use a logarithmic link function, combined with an indirect genetic effects model. However, in the presence of highly infectious individuals, i.e. super-spreaders, the use of any model that is linear in susceptibility and infectiousness causes biased estimates. Thus, in order to identify super-spreaders, novel analytical methods using our derived expression are required. We have derived a genetic-epidemiological function for quantitative genetic analyses of binary infectious disease data, which, unlike current approaches, takes infection dynamics into account and allows for variation in host susceptibility and infectiousness.

Optimum design of structures subject to general periodic loads

NASA Technical Reports Server (NTRS)

Reiss, Robert; Qian, B.

1989-01-01

A simplified version of Icerman's problem regarding the design of structures subject to a single harmonic load is discussed. The nature of the restrictive conditions that must be placed on the design space in order to ensure an analytic optimum are discussed in detail. Icerman's problem is then extended to include multiple forcing functions with different driving frequencies. And the conditions that now must be placed upon the design space to ensure an analytic optimum are again discussed. An important finding is that all solutions to the optimality condition (analytic stationary design) are local optima, but the global optimum may well be non-analytic. The more general problem of distributing the fixed mass of a linear elastic structure subject to general periodic loads in order to minimize some measure of the steady state deflection is also considered. This response is explicitly expressed in terms of Green's functional and the abstract operators defining the structure. The optimality criterion is derived by differentiating the response with respect to the design parameters. The theory is applicable to finite element as well as distributed parameter models.

Analytic and numeric Green's functions for a two-dimensional electron gas in an orthogonal magnetic field

DOE Office of Scientific and Technical Information (OSTI.GOV)

Cresti, Alessandro; Grosso, Giuseppe; Parravicini, Giuseppe Pastori

2006-05-15

We have derived closed analytic expressions for the Green's function of an electron in a two-dimensional electron gas threaded by a uniform perpendicular magnetic field, also in the presence of a uniform electric field and of a parabolic spatial confinement. A workable and powerful numerical procedure for the calculation of the Green's functions for a large infinitely extended quantum wire is considered exploiting a lattice model for the wire, the tight-binding representation for the corresponding matrix Green's function, and the Peierls phase factor in the Hamiltonian hopping matrix element to account for the magnetic field. The numerical evaluation of themore » Green's function has been performed by means of the decimation-renormalization method, and quite satisfactorily compared with the analytic results worked out in this paper. As an example of the versatility of the numerical and analytic tools here presented, the peculiar semilocal character of the magnetic Green's function is studied in detail because of its basic importance in determining magneto-transport properties in mesoscopic systems.« less

Analytic Reflected Lightcurves for Exoplanets

NASA Astrophysics Data System (ADS)

Haggard, Hal M.; Cowan, Nicolas B.

2018-04-01

The disk-integrated reflected brightness of an exoplanet changes as a function of time due to orbital and rotational motion coupled with an inhomogeneous albedo map. We have previously derived analytic reflected lightcurves for spherical harmonic albedo maps in the special case of a synchronously-rotating planet on an edge-on orbit (Cowan, Fuentes & Haggard 2013). In this letter, we present analytic reflected lightcurves for the general case of a planet on an inclined orbit, with arbitrary spin period and non-zero obliquity. We do so for two different albedo basis maps: bright points (δ-maps), and spherical harmonics (Y_l^m-maps). In particular, we use Wigner D-matrices to express an harmonic lightcurve for an arbitrary viewing geometry as a non-linear combination of harmonic lightcurves for the simpler edge-on, synchronously rotating geometry. These solutions will enable future exploration of the degeneracies and information content of reflected lightcurves, as well as fast calculation of lightcurves for mapping exoplanets based on time-resolved photometry. To these ends we make available Exoplanet Analytic Reflected Lightcurves (EARL), a simple open-source code that allows rapid computation of reflected lightcurves.

Point spread functions and deconvolution of ultrasonic images.

PubMed

Dalitz, Christoph; Pohle-Fröhlich, Regina; Michalk, Thorsten

2015-03-01

This article investigates the restoration of ultrasonic pulse-echo C-scan images by means of deconvolution with a point spread function (PSF). The deconvolution concept from linear system theory (LST) is linked to the wave equation formulation of the imaging process, and an analytic formula for the PSF of planar transducers is derived. For this analytic expression, different numerical and analytic approximation schemes for evaluating the PSF are presented. By comparing simulated images with measured C-scan images, we demonstrate that the assumptions of LST in combination with our formula for the PSF are a good model for the pulse-echo imaging process. To reconstruct the object from a C-scan image, we compare different deconvolution schemes: the Wiener filter, the ForWaRD algorithm, and the Richardson-Lucy algorithm. The best results are obtained with the Richardson-Lucy algorithm with total variation regularization. For distances greater or equal twice the near field distance, our experiments show that the numerically computed PSF can be replaced with a simple closed analytic term based on a far field approximation.

Transient Dynamics of Double Quantum Dots Coupled to Two Reservoirs

NASA Astrophysics Data System (ADS)

Fukadai, Takahisa; Sasamoto, Tomohiro

2018-05-01

We study the time-dependent properties of double quantum dots coupled to two reservoirs using the nonequilibrium Green function method. For an arbitrary time-dependent bias, we derive an expression for the time-dependent electron density of a dot and several currents, including the current between the dots in the wide-band-limit approximation. For the special case of a constant bias, we calculate the electron density and the currents numerically. As a result, we find that these quantities oscillate and that the number of crests in a single period of the current from a dot changes with the bias voltage. We also obtain an analytical expression for the relaxation time, which expresses how fast the system converges to its steady state. From the expression, we find that the relaxation time becomes constant when the coupling strength between the dots is sufficiently large in comparison with the difference of coupling strength between the dots and the reservoirs.

An improved 3D MoF method based on analytical partial derivatives

NASA Astrophysics Data System (ADS)

Chen, Xiang; Zhang, Xiong

2016-12-01

MoF (Moment of Fluid) method is one of the most accurate approaches among various surface reconstruction algorithms. As other second order methods, MoF method needs to solve an implicit optimization problem to obtain the optimal approximate surface. Therefore, the partial derivatives of the objective function have to be involved during the iteration for efficiency and accuracy. However, to the best of our knowledge, the derivatives are currently estimated numerically by finite difference approximation because it is very difficult to obtain the analytical derivatives of the object function for an implicit optimization problem. Employing numerical derivatives in an iteration not only increase the computational cost, but also deteriorate the convergence rate and robustness of the iteration due to their numerical error. In this paper, the analytical first order partial derivatives of the objective function are deduced for 3D problems. The analytical derivatives can be calculated accurately, so they are incorporated into the MoF method to improve its accuracy, efficiency and robustness. Numerical studies show that by using the analytical derivatives the iterations are converged in all mixed cells with the efficiency improvement of 3 to 4 times.

General pulsed-field gradient signal attenuation expression based on a fractional integral modified-Bloch equation

NASA Astrophysics Data System (ADS)

Lin, Guoxing

2018-10-01

Anomalous diffusion has been investigated in many polymer and biological systems. The analysis of PFG anomalous diffusion relies on the ability to obtain the signal attenuation expression. However, the general analytical PFG signal attenuation expression based on the fractional derivative has not been previously reported. Additionally, the reported modified-Bloch equations for PFG anomalous diffusion in the literature yielded different results due to their different forms. Here, a new integral type modified-Bloch equation based on the fractional derivative for PFG anomalous diffusion is proposed, which is significantly different from the conventional differential type modified-Bloch equation. The merit of the integral type modified-Bloch equation is that the original properties of the contributions from linear or nonlinear processes remain unchanged at the instant of the combination. From the modified-Bloch equation, the general solutions are derived, which includes the finite gradient pulse width (FGPW) effect. The numerical evaluation of these PFG signal attenuation expressions can be obtained either by the Adomian decomposition, or a direct integration method that is fast and practicable. The theoretical results agree with the continuous-time random walk (CTRW) simulations performed in this paper. Additionally, the relaxation effect in PFG anomalous diffusion is found to be different from that in PFG normal diffusion. The new modified-Bloch equations and their solutions provide a fundamental tool to analyze PFG anomalous diffusion in nuclear magnetic resonance (NMR) and magnetic resonance imaging (MRI).

Two-dimensional dispersion of magnetostatic volume spin waves

NASA Astrophysics Data System (ADS)

Buijnsters, Frank J.; van Tilburg, Lennert J. A.; Fasolino, Annalisa; Katsnelson, Mikhail I.

2018-06-01

Owing to the dipolar (magnetostatic) interaction, long-wavelength spin waves in in-plane magnetized films show an unusual dispersion behavior, which can be mathematically described by the model of and and refinements thereof. However, solving the two-dimensional dispersion requires the evaluation of a set of coupled transcendental equations and one has to rely on numerics. In this work, we present a systematic perturbative analysis of the spin wave model. An expansion in the in-plane wavevector allows us to obtain explicit closed-form expressions for the dispersion relation and mode profiles in various asymptotic regimes. Moreover, we derive a very accurate semi-analytical expression for the dispersion relation of the lowest-frequency mode that is straightforward to evaluate.

Albedo of an irradiated plane-parallel atmosphere with finite optical depth

NASA Astrophysics Data System (ADS)

Fukue, Jun

2018-03-01

We analytically derive albedo for a plane-parallel atmosphere with finite optical depth, irradiated by an external source, under the local thermodynamic equilibrium approximation. Albedo is expressed as a function of the photon destruction probability ɛ and optical depth τ, with several parameters such as dilution factors of the external source. In the particular case of the infinite optical depth, albedo A is expressed as A=[1 + (1-W_J/W_H)√{3ɛ}/3]/(1+√{3ɛ}), where WJ and WH are the dilution factors for the mean intensity and Eddington flux, respectively. An example of a model atmosphere is also presented under a gray approximation.

Scaling laws for AC gas breakdown and implications for universality

NASA Astrophysics Data System (ADS)

Loveless, Amanda M.; Garner, Allen L.

2017-10-01

The reduced dependence on secondary electron emission and electrode surface properties makes radiofrequency (RF) and microwave (MW) plasmas advantageous over direct current (DC) plasmas for various applications, such as microthrusters. Theoretical models relating molecular constants to alternating current (AC) breakdown often fail due to incomplete understanding of both the constants and the mechanisms involved. This work derives simple analytic expressions for RF and MW breakdown, demonstrating the transition between these regimes at their high and low frequency limits, respectively. We further show that the limiting expressions for DC, RF, and MW breakdown voltage all have the same universal scaling dependence on pressure and gap distance at high pressure, agreeing with experiment.

Theory of nonstationary Hawkes processes

NASA Astrophysics Data System (ADS)

Tannenbaum, Neta Ravid; Burak, Yoram

2017-12-01

We expand the theory of Hawkes processes to the nonstationary case, in which the mutually exciting point processes receive time-dependent inputs. We derive an analytical expression for the time-dependent correlations, which can be applied to networks with arbitrary connectivity, and inputs with arbitrary statistics. The expression shows how the network correlations are determined by the interplay between the network topology, the transfer functions relating units within the network, and the pattern and statistics of the external inputs. We illustrate the correlation structure using several examples in which neural network dynamics are modeled as a Hawkes process. In particular, we focus on the interplay between internally and externally generated oscillations and their signatures in the spike and rate correlation functions.

Electromagnetic plane-wave pulse transmission into a Lorentz half-space.

PubMed

Cartwright, Natalie A

2011-12-01

The propagation of an electromagnetic plane-wave signal obliquely incident upon a Lorentz half-space is studied analytically. Time-domain asymptotic expressions that increase in accuracy with propagation distance are derived by application of uniform saddle point methods on the Fourier-Laplace integral representation of the transmitted field. The results are shown to be continuous in time and comparable with numerical calculations of the field. Arrival times and angles of refraction are given for prominent transient pulse features and the steady-state signal.

Møller-Plesset perturbation theory gradient in the generalized hybrid orbital quantum mechanical and molecular mechanical method

NASA Astrophysics Data System (ADS)

Jung, Jaewoon; Sugita, Yuji; Ten-no, S.

2010-02-01

An analytic gradient expression is formulated and implemented for the second-order Møller-Plesset perturbation theory (MP2) based on the generalized hybrid orbital QM/MM method. The method enables us to obtain an accurate geometry at a reasonable computational cost. The performance of the method is assessed for various isomers of alanine dipepetide. We also compare the optimized structures of fumaramide-derived [2]rotaxane and cAMP-dependent protein kinase with experiment.

Parametric number covariance in quantum chaotic spectra.

PubMed

Vinayak; Kumar, Sandeep; Pandey, Akhilesh

2016-03-01

We study spectral parametric correlations in quantum chaotic systems and introduce the number covariance as a measure of such correlations. We derive analytic results for the classical random matrix ensembles using the binary correlation method and obtain compact expressions for the covariance. We illustrate the universality of this measure by presenting the spectral analysis of the quantum kicked rotors for the time-reversal invariant and time-reversal noninvariant cases. A local version of the parametric number variance introduced earlier is also investigated.

Effects of surface roughness and absorption on light propagation in graded-profile waveguides

DOE Office of Scientific and Technical Information (OSTI.GOV)

Danilenko, S S; Osovitskii, A N

2011-06-30

This paper examines the effects of surface roughness and absorption on laser light propagation in graded-profile waveguiding structures. We derive analytical expressions for the scattering and absorption coefficients of guided waves and analyse these coefficients in relation to parameters of the waveguiding structure and the roughness of its boundary. A new approach is proposed to measuring roughness parameters of precision dielectric surfaces. Experimental evidence is presented which supports the main conclusions of the theory. (integraled-optical waweguides)

Tsunami and acoustic-gravity waves in water of constant depth

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hendin, Gali; Stiassnie, Michael

2013-08-15

A study of wave radiation by a rather general bottom displacement, in a compressible ocean of otherwise constant depth, is carried out within the framework of a three-dimensional linear theory. Simple analytic expressions for the flow field, at large distance from the disturbance, are derived. Realistic numerical examples indicate that the Acoustic-Gravity waves, which significantly precede the Tsunami, are expected to leave a measurable signature on bottom-pressure records that should be considered for early detection of Tsunami.

Improvement of the accuracy of noise measurements by the two-amplifier correlation method.

PubMed

Pellegrini, B; Basso, G; Fiori, G; Macucci, M; Maione, I A; Marconcini, P

2013-10-01

We present a novel method for device noise measurement, based on a two-channel cross-correlation technique and a direct "in situ" measurement of the transimpedance of the device under test (DUT), which allows improved accuracy with respect to what is available in the literature, in particular when the DUT is a nonlinear device. Detailed analytical expressions for the total residual noise are derived, and an experimental investigation of the increased accuracy provided by the method is performed.

Nonlinear propagation of electromagnetic waves in negative-refraction-index composite materials.

PubMed

Kourakis, I; Shukla, P K

2005-07-01

We investigate the nonlinear propagation of electromagnetic waves in left-handed materials. For this purpose, we consider a set of coupled nonlinear Schrödinger (CNLS) equations, which govern the dynamics of coupled electric and magnetic field envelopes. The CNLS equations are used to obtain a nonlinear dispersion, which depicts the modulational stability profile of the coupled plane-wave solutions in left-handed materials. An exact (in)stability criterion for modulational interactions is derived, and analytical expressions for the instability growth rate are obtained.

Thermal convection in a cylindrical enclosure

NASA Astrophysics Data System (ADS)

Shukla, K. N.

2005-02-01

The paper highlights the onset of convection in a fluid layer partially filled in an axisymmetric container. The equilibrium of the fluid is disturbed with the deformation of the interface due to residual acceleration. The general problem of deformable interface involves a dimensionless parameter, the Bond number. An analytical expression for the natural frequencies of the deformable surface is derived in terms of the Bond number, which determines the time period required for the stable location of the fluid for the propellant management of the spacecraft.

Optical response of photopolymer materials for holographic data storage applications.

PubMed

Sheridan, J T; Gleeson, M R; Close, C E; Kelly, J V

2007-01-01

We briefly review the application of photopolymer recording materials in the area of holographic data storage. In particular we discuss the recent development of the Non-local Polymerisation Driven Diffusion model. Applying this model we develop simple first-order analytic expressions describing the spatial frequency response of photopolymer materials. The assumptions made in the derivation of these formulae are described and their ranges of validity are examined. The effects of particular physical parameters of a photopolymer on the material response are discussed.

Hilbert-Schmidt Measure of Pairwise Quantum Discord for Three-Qubit X States

NASA Astrophysics Data System (ADS)

Daoud, M.; Laamara, R. Ahl; Seddik, S.

2015-10-01

The Hilbert-Schmidt distance between a mixed three-qubit state and its closest state is used to quantify the amount of pairwise quantum correlations in a tripartite system. Analytical expressions of geometric quantum discord are derived. A particular attention is devoted to two special classes of three-qubit X states. They include three-qubit states of W, GHZ and Bell type. We also discuss the monogamy property of geometric quantum discord in some mixed three-qubit systems.

Analytic Expressions for the Inner-rim Structure of Passively Heated Protoplanetary Disks

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ueda, Takahiro; Okuzumi, Satoshi; Flock, Mario, E-mail: t_ueda@geo.titech.ac.jp

We analytically derive the expressions for the structure of the inner region of protoplanetary disks based on the results from the recent hydrodynamical simulations. The inner part of a disk can be divided into four regions: a dust-free region with a gas temperature in the optically thin limit, an optically thin dust halo, an optically thick condensation front, and the classical, optically thick region, in order from the innermost to the outermost. We derive the dust-to-gas mass ratio profile in the dust halo using the fact that partial dust condensation regulates the temperature relative to the dust evaporation temperature. Beyondmore » the dust halo, there is an optically thick condensation front where all the available silicate gas condenses out. The curvature of the condensation surface is determined by the condition that the surface temperature must be nearly equal to the characteristic temperature ∼1200 K. We derive the midplane temperature in the outer two regions using the two-layer approximation, with the additional heating by the condensation front for the outermost region. As a result, the overall temperature profile is step-like, with steep gradients at the borders between the outer three regions. The borders might act as planet traps where the inward migration of planets due to gravitational interaction with the gas disk stops. The temperature at the border between the two outermost regions coincides with the temperature needed to activate magnetorotational instability, suggesting that the inner edge of the dead zone must lie at this border. The radius of the dead zone inner edge predicted from our solution is ∼2–3 times larger than that expected from the classical optically thick temperature.« less

Analytic Expressions for the Inner-rim Structure of Passively Heated Protoplanetary Disks

NASA Astrophysics Data System (ADS)

Ueda, Takahiro; Okuzumi, Satoshi; Flock, Mario

2017-07-01

We analytically derive the expressions for the structure of the inner region of protoplanetary disks based on the results from the recent hydrodynamical simulations. The inner part of a disk can be divided into four regions: a dust-free region with a gas temperature in the optically thin limit, an optically thin dust halo, an optically thick condensation front, and the classical, optically thick region, in order from the innermost to the outermost. We derive the dust-to-gas mass ratio profile in the dust halo using the fact that partial dust condensation regulates the temperature relative to the dust evaporation temperature. Beyond the dust halo, there is an optically thick condensation front where all the available silicate gas condenses out. The curvature of the condensation surface is determined by the condition that the surface temperature must be nearly equal to the characteristic temperature ˜1200 K. We derive the midplane temperature in the outer two regions using the two-layer approximation, with the additional heating by the condensation front for the outermost region. As a result, the overall temperature profile is step-like, with steep gradients at the borders between the outer three regions. The borders might act as planet traps where the inward migration of planets due to gravitational interaction with the gas disk stops. The temperature at the border between the two outermost regions coincides with the temperature needed to activate magnetorotational instability, suggesting that the inner edge of the dead zone must lie at this border. The radius of the dead zone inner edge predicted from our solution is ˜2-3 times larger than that expected from the classical optically thick temperature.

Reverie and metaphor. Some thoughts on how I work as a psychoanalyst.

PubMed

Ogden, T

1997-08-01

In this paper, the author presents parts of an ongoing internal dialogue concerning how he works as an analyst. He describes the way in which he attempts to sense what is most alive and most real in each analytic encounter, as well as his use of his own reveries in his effort to locate himself in what is going on at an unconscious level in the analytic relationship. The author views each analytic situation as reflecting, to a large degree, a specific type of unconscious intersubjective construction. Since unconscious experience is by definition outside of conscious awareness, the analyst must make use of indirect (associational) methods such as the scrutiny of his own reverie experience in his efforts to 'catch the drift' (Freud, 1923, p. 239) of the unconscious intersubjective constructions being generated. Reveries (and all other derivatives of the unconscious) are viewed not as glimpses into the unconscious, but as metaphorical expressions of what the unconscious experience is like. In the author's experience, when an analysis is 'a going concern', the analytic dialogue often takes the form of a verbal 'squiggle game' (Winnicott, 1971a, p. 3) in which the analytic pair elaborates and modifies the metaphors that the other has unself-consciously introduced. The analytic use of reverie and of the role of metaphor in the analytic experience is clinically illustrated.

Excitation, response, and fatigue life estimation methods for the structural design of externally blown flaps

NASA Technical Reports Server (NTRS)

Ungar, E. E.; Chandiramani, K. L.; Barger, J. E.

1972-01-01

Means for predicting the fluctuating pressures acting on externally blown flap surfaces are developed on the basis of generalizations derived from non-dimensionalized empirical data. Approaches for estimation of the fatigue lives of skin-stringer and honeycomb-core sandwich flap structures are derived from vibration response analyses and panel fatigue data. Approximate expressions for fluctuating pressures, structural response, and fatigue life are combined to reveal the important parametric dependences. The two-dimensional equations of motion of multi-element flap systems are derived in general form, so that they can be specialized readily for any particular system. An introduction is presented of an approach to characterizing the excitation pressures and structural responses which makes use of space-time spectral concepts and promises to provide useful insights, as well as experimental and analytical savings.

Analytic Expressions for the Gravity Gradient Tensor of 3D Prisms with Depth-Dependent Density

NASA Astrophysics Data System (ADS)

Jiang, Li; Liu, Jie; Zhang, Jianzhong; Feng, Zhibing

2017-12-01

Variable-density sources have been paid more attention in gravity modeling. We conduct the computation of gravity gradient tensor of given mass sources with variable density in this paper. 3D rectangular prisms, as simple building blocks, can be used to approximate well 3D irregular-shaped sources. A polynomial function of depth can represent flexibly the complicated density variations in each prism. Hence, we derive the analytic expressions in closed form for computing all components of the gravity gradient tensor due to a 3D right rectangular prism with an arbitrary-order polynomial density function of depth. The singularity of the expressions is analyzed. The singular points distribute at the corners of the prism or on some of the lines through the edges of the prism in the lower semi-space containing the prism. The expressions are validated, and their numerical stability is also evaluated through numerical tests. The numerical examples with variable-density prism and basin models show that the expressions within their range of numerical stability are superior in computational accuracy and efficiency to the common solution that sums up the effects of a collection of uniform subprisms, and provide an effective method for computing gravity gradient tensor of 3D irregular-shaped sources with complicated density variation. In addition, the tensor computed with variable density is different in magnitude from that with constant density. It demonstrates the importance of the gravity gradient tensor modeling with variable density.

Motion of vortices in inhomogeneous Bose-Einstein condensates

NASA Astrophysics Data System (ADS)

Groszek, Andrew J.; Paganin, David M.; Helmerson, Kristian; Simula, Tapio P.

2018-02-01

We derive a general and exact equation of motion for a quantized vortex in an inhomogeneous two-dimensional Bose-Einstein condensate. This equation expresses the velocity of a vortex as a sum of local ambient density and phase gradients in the vicinity of the vortex. We perform Gross-Pitaevskii simulations of single-vortex dynamics in both harmonic and hard-walled disk-shaped traps, and find excellent agreement in both cases with our analytical prediction. The simulations reveal that, in a harmonic trap, the main contribution to the vortex velocity is an induced ambient phase gradient, a finding that contradicts the commonly quoted result that the local density gradient is the only relevant effect in this scenario. We use our analytical vortex velocity formula to derive a point-vortex model that accounts for both density and phase contributions to the vortex velocity, suitable for use in inhomogeneous condensates. Although good agreement is obtained between Gross-Pitaevskii and point-vortex simulations for specific few-vortex configurations, the effects of nonuniform condensate density are in general highly nontrivial, and are thus difficult to efficiently and accurately model using a simplified point-vortex description.

Vortex Core Size in the Rotor Near-Wake

NASA Technical Reports Server (NTRS)

Young, Larry A.

2003-01-01

Using a kinetic energy conservation approach, a number of simple analytic expressions are derived for estimating the core size of tip vortices in the near-wake of rotors in hover and axial-flow flight. The influence of thrust, induced power losses, advance ratio, and vortex structure on rotor vortex core size is assessed. Experimental data from the literature is compared to the analytical results derived in this paper. In general, three conclusions can be drawn from the work in this paper. First, the greater the rotor thrust, t h e larger the vortex core size in the rotor near-wake. Second, the more efficient a rotor is with respect to induced power losses, the smaller the resulting vortex core size. Third, and lastly, vortex core size initially decreases for low axial-flow advance ratios, but for large advance ratios core size asymptotically increases to a nominal upper limit. Insights gained from this work should enable improved modeling of rotary-wing aerodynamics, as well as provide a framework for improved experimental investigations of rotor a n d advanced propeller wakes.

Scientists' sense making when hypothesizing about disease mechanisms from expression data and their needs for visualization support.

PubMed

Mirel, Barbara; Görg, Carsten

2014-04-26

A common class of biomedical analysis is to explore expression data from high throughput experiments for the purpose of uncovering functional relationships that can lead to a hypothesis about mechanisms of a disease. We call this analysis expression driven, -omics hypothesizing. In it, scientists use interactive data visualizations and read deeply in the research literature. Little is known, however, about the actual flow of reasoning and behaviors (sense making) that scientists enact in this analysis, end-to-end. Understanding this flow is important because if bioinformatics tools are to be truly useful they must support it. Sense making models of visual analytics in other domains have been developed and used to inform the design of useful and usable tools. We believe they would be helpful in bioinformatics. To characterize the sense making involved in expression-driven, -omics hypothesizing, we conducted an in-depth observational study of one scientist as she engaged in this analysis over six months. From findings, we abstracted a preliminary sense making model. Here we describe its stages and suggest guidelines for developing visualization tools that we derived from this case. A single case cannot be generalized. But we offer our findings, sense making model and case-based tool guidelines as a first step toward increasing interest and further research in the bioinformatics field on scientists' analytical workflows and their implications for tool design.

Scientists’ sense making when hypothesizing about disease mechanisms from expression data and their needs for visualization support

PubMed Central

2014-01-01

A common class of biomedical analysis is to explore expression data from high throughput experiments for the purpose of uncovering functional relationships that can lead to a hypothesis about mechanisms of a disease. We call this analysis expression driven, -omics hypothesizing. In it, scientists use interactive data visualizations and read deeply in the research literature. Little is known, however, about the actual flow of reasoning and behaviors (sense making) that scientists enact in this analysis, end-to-end. Understanding this flow is important because if bioinformatics tools are to be truly useful they must support it. Sense making models of visual analytics in other domains have been developed and used to inform the design of useful and usable tools. We believe they would be helpful in bioinformatics. To characterize the sense making involved in expression-driven, -omics hypothesizing, we conducted an in-depth observational study of one scientist as she engaged in this analysis over six months. From findings, we abstracted a preliminary sense making model. Here we describe its stages and suggest guidelines for developing visualization tools that we derived from this case. A single case cannot be generalized. But we offer our findings, sense making model and case-based tool guidelines as a first step toward increasing interest and further research in the bioinformatics field on scientists’ analytical workflows and their implications for tool design. PMID:24766796

Analytic relations for magnifications and time delays in gravitational lenses with fold and cusp configurations

NASA Astrophysics Data System (ADS)

Congdon, Arthur B.; Keeton, Charles R.; Nordgren, C. Erik

2008-09-01

Gravitational lensing provides a unique and powerful probe of the mass distributions of distant galaxies. Four-image lens systems with fold and cusp configurations have two or three bright images near a critical point. Within the framework of singularity theory, we derive analytic relations that are satisfied for a light source that lies a small but finite distance from the astroid caustic of a four-image lens. Using a perturbative expansion of the image positions, we show that the time delay between the close pair of images in a fold lens scales with the cube of the image separation, with a constant of proportionality that depends on a particular third derivative of the lens potential. We also apply our formalism to cusp lenses, where we develop perturbative expressions for the image positions, magnifications and time delays of the images in a cusp triplet. Some of these results were derived previously for a source asymptotically close to a cusp point, but using a simplified form of the lens equation whose validity may be in doubt for sources that lie at astrophysically relevant distances from the caustic. Along with the work of Keeton, Gaudi & Petters, this paper demonstrates that perturbation theory plays an important role in theoretical lensing studies.

Laser pulse propagation in inhomogeneous magnetoplasma channels and wakefield acceleration

NASA Astrophysics Data System (ADS)

Sharma, B. S.; Jain, Archana; Jaiman, N. K.; Gupta, D. N.; Jang, D. G.; Suk, H.; Kulagin, V. V.

2014-02-01

Wakefield excitation in a preformed inhomogeneous parabolic plasma channel by an intense relativistic (≃1019 W/cm2) circularly polarized Gaussian laser pulse is investigated analytically and numerically in the presence of an external longitudinal magnetic field. A three dimensional envelope equation for the evolution of the laser pulse is derived, which includes the effect of the nonparaxial and applied external magnetic field. A relation for the channel radius with the laser spot size is derived and examines numerically to see the external magnetic field effect. It is observed that the channel radius depends on the applied external magnetic field. An analytical expression for the wakefield is derived and validated with the help of a two dimensional particle in cell (2D PIC) simulation code. It is shown that the electromagnetic nature of the wakes in an inhomogeneous plasma channel makes their excitation nonlocal, which results in change of fields with time and external magnetic field due to phase mixing of the plasma oscillations with spatially varying frequencies. The magnetic field effect on perturbation of the plasma density and decreasing length is also analyzed numerically. In addition, it has been shown that the electron energy gain in the inhomogeneous parabolic magnetoplasma channel can be increased significantly compared with the homogeneous plasma channel.

On the effect of the degeneracy among dark energy parameters

NASA Astrophysics Data System (ADS)

Gong, Yungui; Gao, Qing

2014-01-01

The dynamics of scalar fields as dark energy is well approximated by some general relations between the equation of state parameter and the fractional energy density . Based on the approximation, for slowly rolling scalar fields, we derived the analytical expressions of which reduce to the popular Chevallier-Polarski-Linder parametrization with an explicit degeneracy relation between and . The models approximate the dynamics of scalar fields well and help eliminate the degeneracies among , , and . With the explicit degeneracy relations, we test their effects on the constraints of the cosmological parameters. We find that: (1) The analytical relations between and for the two models are consistent with observational data. (2) The degeneracies have little effect on . (3) The error of was reduced about 30 % with the degeneracy relations.

Quantum Effects at a Proton Relaxation at Low Temperatures

NASA Astrophysics Data System (ADS)

Kalytka, V. A.; Korovkin, M. V.

2016-11-01

Quantum effects during migratory polarization in multi-well crystals (including multi-well silicates and crystalline hydrates) are investigated in a variable electric field at low temperatures by direct quantum-mechanical calculations. Based on analytical solution of the quantum Liouville kinetic equation in the linear approximation for the polarizing field, the non-stationary density matrix is calculated for an ensemble of non-interacting protons moving in the field of one-dimensional multi-well crystal potential relief of rectangular shape. An expression for the complex dielectric constant convenient for a comparison with experiment and calculation of relaxer parameters is derived using the nonequilibrium polarization density matrix. The density matrix apparatus can be used for analytical investigation of the quantum mechanism of spontaneous polarization of a ferroelectric material (KDP and DKDP).

Plate and butt-weld stresses beyond elastic limit, material and structural modeling

NASA Technical Reports Server (NTRS)

Verderaime, V.

1991-01-01

Ultimate safety factors of high performance structures depend on stress behavior beyond the elastic limit, a region not too well understood. An analytical modeling approach was developed to gain fundamental insights into inelastic responses of simple structural elements. Nonlinear material properties were expressed in engineering stresses and strains variables and combined with strength of material stress and strain equations similar to numerical piece-wise linear method. Integrations are continuous which allows for more detailed solutions. Included with interesting results are the classical combined axial tension and bending load model and the strain gauge conversion to stress beyond the elastic limit. Material discontinuity stress factors in butt-welds were derived. This is a working-type document with analytical methods and results applicable to all industries of high reliability structures.

Non-collinear magnetism with analytic Bond-Order Potentials

NASA Astrophysics Data System (ADS)

Ford, Michael E.; Pettifor, D. G.; Drautz, Ralf

2015-03-01

The theory of analytic Bond-Order Potentials as applied to non-collinear magnetic structures of transition metals is extended to take into account explicit rotations of Hamiltonian and local moment matrix elements between locally and globally defined spin-coordinate systems. Expressions for the gradients of the energy with respect to the Hamiltonian matrix elements, the interatomic forces and the magnetic torques are derived. The method is applied to simulations of the rotation of magnetic moments in α iron, as well as α and β manganese, based on d-valent orthogonal tight-binding parametrizations of the electronic structure. A new weighted-average terminator is introduced to improve the convergence of the Bond-Order Potential energies and torques with respect to tight-binding reference values, although the general behavior is qualitatively correct for low-moment expansions.

Derivation of Markov processes that violate detailed balance

NASA Astrophysics Data System (ADS)

Lee, Julian

2018-03-01

Time-reversal symmetry of the microscopic laws dictates that the equilibrium distribution of a stochastic process must obey the condition of detailed balance. However, cyclic Markov processes that do not admit equilibrium distributions with detailed balance are often used to model systems driven out of equilibrium by external agents. I show that for a Markov model without detailed balance, an extended Markov model can be constructed, which explicitly includes the degrees of freedom for the driving agent and satisfies the detailed balance condition. The original cyclic Markov model for the driven system is then recovered as an approximation at early times by summing over the degrees of freedom for the driving agent. I also show that the widely accepted expression for the entropy production in a cyclic Markov model is actually a time derivative of an entropy component in the extended model. Further, I present an analytic expression for the entropy component that is hidden in the cyclic Markov model.

Stability of surface nanobubbles

NASA Astrophysics Data System (ADS)

Maheshwari, Shantanu; van der Hoef, Martin; Zhang, Xuehua; Lohse, Detlef

2015-11-01

We have studied the stability and dissolution of surface nanobubbles on the chemical heterogenous surface by performing Molecular Dynamics (MD) simulations of binary mixture consists of Lennard-Jones (LJ) particles. Recently our group has derived the exact expression for equilibrium contact angle of surface nanobubbles as a function of oversaturation of the gas concentration in bulk liquid and the lateral length of bubble. It has been showed that the contact line pinning and the oversaturation of gas concentration in bulk liquid is crucial in the stability of surface nanobubbles. Our simulations showed that how pinning of the three-phase contact line on the chemical heterogenous surface lead to the stability of the nanobubble. We have calculated the equilibrium contact angle by varying the gas concentration in bulk liquid and the lateral length of the bubble. Our results showed that the equilibrium contact angle follows the expression derived analytically by our group. We have also studied the bubble dissolution dynamics and showed the ''stick-jump'' mechanism which was also observed experimentally in case of dissolution of nanodrops.

Acoustic field of a wedge-shaped section of a spherical cap transducer

NASA Astrophysics Data System (ADS)

Ketterling, Jeffrey A.

2003-12-01

The acoustic pressure field at an arbitrary point in space is derived for a wedge-shaped section of a spherical cap transducer using the spatial impulse response (SIR) method. For a spherical surface centered at the origin, a wedge shape is created by taking cuts in the X-Y and X-Z planes and removing the smallest surface component. Analytic expressions are derived for the SIR based on spatial location. The expressions utilize the SIR solutions for a spherical cap transducer [Arditi et al., Ultrason. Imaging 3, 37-61 (1981)] with additional terms added to account for the reduced surface area of the wedge. Results from the numerical model are compared to experimental measurements from a wedge transducer with an 8-cm outer diameter and 9-cm geometric focus. The experimental and theoretical -3-dB beamwidths agreed to within 10%+/-5%. The SIR model for a wedge-shaped transducer is easily extended to other spherically curved transducer geometries that consist of combinations of wedge sections and spherical caps.

Acoustic field of a wedge-shaped section of a spherical cap transducer.

PubMed

Ketterling, Jeffrey A

2003-12-01

The acoustic pressure field at an arbitrary point in space is derived for a wedge-shaped section of a spherical cap transducer using the spatial impulse response (SIR) method. For a spherical surface centered at the origin, a wedge shape is created by taking cuts in the X-Y and X-Z planes and removing the smallest surface component. Analytic expressions are derived for the SIR based on spatial location. The expressions utilize the SIR solutions for a spherical cap transducer [Arditi et al., Ultrason. Imaging 3, 37-61 (1981)] with additional terms added to account for the reduced surface area of the wedge. Results from the numerical model are compared to experimental measurements from a wedge transducer with an 8-cm outer diameter and 9-cm geometric focus. The experimental and theoretical -3-dB beamwidths agreed to within 10% +/- 5%. The SIR model for a wedge-shaped transducer is easily extended to other spherically curved transducer geometries that consist of combinations of wedge sections and spherical caps.

Simple phenomenological modeling of transition-region capacitance of forward-biased p-n junction diodes and transistor diodes

NASA Technical Reports Server (NTRS)

Lindholm, F. A.

1982-01-01

The derivation of a simple expression for the capacitance C(V) associated with the transition region of a p-n junction under a forward bias is derived by phenomenological reasoning. The treatment of C(V) is based on the conventional Shockley equations, and simpler expressions for C(V) result that are in general accord with the previous analytical and numerical results. C(V) consists of two components resulting from changes in majority carrier concentration and from free hole and electron accumulation in the space-charge region. The space-charge region is conceived as the intrinsic region of an n-i-p structure for a space-charge region markedly wider than the extrinsic Debye lengths at its edges. This region is excited in the sense that the forward bias creates hole and electron densities orders of magnitude larger than those in equilibrium. The recent Shirts-Gordon (1979) modeling of the space-charge region using a dielectric response function is contrasted with the more conventional Schottky-Shockley modeling.

Level crossings and excess times due to a superposition of uncorrelated exponential pulses

NASA Astrophysics Data System (ADS)

Theodorsen, A.; Garcia, O. E.

2018-01-01

A well-known stochastic model for intermittent fluctuations in physical systems is investigated. The model is given by a superposition of uncorrelated exponential pulses, and the degree of pulse overlap is interpreted as an intermittency parameter. Expressions for excess time statistics, that is, the rate of level crossings above a given threshold and the average time spent above the threshold, are derived from the joint distribution of the process and its derivative. Limits of both high and low intermittency are investigated and compared to previously known results. In the case of a strongly intermittent process, the distribution of times spent above threshold is obtained analytically. This expression is verified numerically, and the distribution of times above threshold is explored for other intermittency regimes. The numerical simulations compare favorably to known results for the distribution of times above the mean threshold for an Ornstein-Uhlenbeck process. This contribution generalizes the excess time statistics for the stochastic model, which find applications in a wide diversity of natural and technological systems.

Analytical approximations for spatial stochastic gene expression in single cells and tissues

PubMed Central

Smith, Stephen; Cianci, Claudia; Grima, Ramon

2016-01-01

Gene expression occurs in an environment in which both stochastic and diffusive effects are significant. Spatial stochastic simulations are computationally expensive compared with their deterministic counterparts, and hence little is currently known of the significance of intrinsic noise in a spatial setting. Starting from the reaction–diffusion master equation (RDME) describing stochastic reaction–diffusion processes, we here derive expressions for the approximate steady-state mean concentrations which are explicit functions of the dimensionality of space, rate constants and diffusion coefficients. The expressions have a simple closed form when the system consists of one effective species. These formulae show that, even for spatially homogeneous systems, mean concentrations can depend on diffusion coefficients: this contradicts the predictions of deterministic reaction–diffusion processes, thus highlighting the importance of intrinsic noise. We confirm our theory by comparison with stochastic simulations, using the RDME and Brownian dynamics, of two models of stochastic and spatial gene expression in single cells and tissues. PMID:27146686

Distributed parameter modeling to prevent charge cancellation for discrete thickness piezoelectric energy harvester

NASA Astrophysics Data System (ADS)

Krishnasamy, M.; Qian, Feng; Zuo, Lei; Lenka, T. R.

2018-03-01

The charge cancellation due to the change of strain along single continuous piezoelectric layer can remarkably affect the performance of a cantilever based harvester. In this paper, analytical models using distributed parameters are developed with some extent of averting the charge cancellation in cantilever piezoelectric transducer where the piezoelectric layers are segmented at strain nodes of concerned vibration mode. The electrode of piezoelectric segments are parallelly connected with a single external resistive load in the 1st model (Model 1). While each bimorph piezoelectric layers are connected in parallel to a resistor to form an independent circuit in the 2nd model (Model 2). The analytical expressions of the closed-form electromechanical coupling responses in frequency domain under harmonic base excitation are derived based on the Euler-Bernoulli beam assumption for both models. The developed analytical models are validated by COMSOL and experimental results. The results demonstrate that the energy harvesting performance of the developed segmented piezoelectric layer models is better than the traditional model of continuous piezoelectric layer.

Highly accurate analytic formulae for projectile motion subjected to quadratic drag

NASA Astrophysics Data System (ADS)

Turkyilmazoglu, Mustafa

2016-05-01

The classical phenomenon of motion of a projectile fired (thrown) into the horizon through resistive air charging a quadratic drag onto the object is revisited in this paper. No exact solution is known that describes the full physical event under such an exerted resistance force. Finding elegant analytical approximations for the most interesting engineering features of dynamical behavior of the projectile is the principal target. Within this purpose, some analytical explicit expressions are derived that accurately predict the maximum height, its arrival time as well as the flight range of the projectile at the highest ascent. The most significant property of the proposed formulas is that they are not restricted to the initial speed and firing angle of the object, nor to the drag coefficient of the medium. In combination with the available approximations in the literature, it is possible to gain information about the flight and complete the picture of a trajectory with high precision, without having to numerically simulate the full governing equations of motion.

Electron transfer statistics and thermal fluctuations in molecular junctions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Goswami, Himangshu Prabal; Harbola, Upendra

2015-02-28

We derive analytical expressions for probability distribution function (PDF) for electron transport in a simple model of quantum junction in presence of thermal fluctuations. Our approach is based on the large deviation theory combined with the generating function method. For large number of electrons transferred, the PDF is found to decay exponentially in the tails with different rates due to applied bias. This asymmetry in the PDF is related to the fluctuation theorem. Statistics of fluctuations are analyzed in terms of the Fano factor. Thermal fluctuations play a quantitative role in determining the statistics of electron transfer; they tend tomore » suppress the average current while enhancing the fluctuations in particle transfer. This gives rise to both bunching and antibunching phenomena as determined by the Fano factor. The thermal fluctuations and shot noise compete with each other and determine the net (effective) statistics of particle transfer. Exact analytical expression is obtained for delay time distribution. The optimal values of the delay time between successive electron transfers can be lowered below the corresponding shot noise values by tuning the thermal effects.« less

Path-integral formalism for stochastic resetting: Exactly solved examples and shortcuts to confinement

NASA Astrophysics Data System (ADS)

Roldán, Édgar; Gupta, Shamik

2017-08-01

We study the dynamics of overdamped Brownian particles diffusing in conservative force fields and undergoing stochastic resetting to a given location at a generic space-dependent rate of resetting. We present a systematic approach involving path integrals and elements of renewal theory that allows us to derive analytical expressions for a variety of statistics of the dynamics such as (i) the propagator prior to first reset, (ii) the distribution of the first-reset time, and (iii) the spatial distribution of the particle at long times. We apply our approach to several representative and hitherto unexplored examples of resetting dynamics. A particularly interesting example for which we find analytical expressions for the statistics of resetting is that of a Brownian particle trapped in a harmonic potential with a rate of resetting that depends on the instantaneous energy of the particle. We find that using energy-dependent resetting processes is more effective in achieving spatial confinement of Brownian particles on a faster time scale than performing quenches of parameters of the harmonic potential.

Reflection and Transmission of P-Waves in an Intermediate Layer Lying Between Two Semi-infinite Media

NASA Astrophysics Data System (ADS)

Singh, Pooja; Chattopadhyay, Amares; Srivastava, Akanksha; Singh, Abhishek Kumar

2018-05-01

With a motivation to gain physical insight of reflection as well as transmission phenomena in frozen (river/ocean) situation for example in Antarctica and other coldest place on Earth, the present article undertakes the analysis of reflection and transmission of a plane wave at the interfaces of layered structured comprised of a water layer of finite thickness sandwiched between an upper half-space constituted of ice and a lower isotropic elastic half-space, which may be useful in geophysical exploration in such conditions. A closed form expression of reflection/transmission coefficients of reflected and transmitted waves has been derived in terms of angles of incidence, propagation vector, displacement vector and elastic constants of the media. Expressions corresponding to the energy partition of various reflected and transmitted waves have also been established analytically. It has been remarkably shown that the law of conservation of energy holds good in the entire reflection and transmission phenomena for different angles of incidence. A numerical examples were performed so to graphically portray the analytical findings. Further the deduced results are validated with the pre-established classical results.

Correlation of basic TL, OSL and IRSL properties of ten K-feldspar samples of various origins

NASA Astrophysics Data System (ADS)

Sfampa, I. K.; Polymeris, G. S.; Pagonis, V.; Theodosoglou, E.; Tsirliganis, N. C.; Kitis, G.

2015-09-01

Feldspars stand among the most widely used minerals in dosimetric methods of dating using thermoluminescence (TL), optically stimulated luminescence (OSL) and infrared stimulated luminescence (IRSL). Having very good dosimetric properties, they can in principle contribute to the dating of every site of archaeological and geological interest. The present work studies basic properties of ten naturally occurring K-feldspar samples belonging to three feldspar species, namely sanidine, orthoclase and microcline. The basic properties studied are (a) the influence of blue light and infrared stimulation on the thermoluminescence glow-curves, (b) the growth of OSL, IRSL, residual TL and TL-loss as a function of OSL and IRSL bleaching time and (c) the correlation between the OSL and IRSL signals and the energy levels responsible for the TL glow-curve. All experimental data were fitted using analytical expressions derived from a recently developed tunneling recombination model. The results show that the analytical expressions provide excellent fits to all experimental results, thus verifying the tunneling recombination mechanism in these materials and providing valuable information about the concentrations of luminescence centers.

Integrate-and-fire neurons driven by asymmetric dichotomous noise.

PubMed

Droste, Felix; Lindner, Benjamin

2014-12-01

We consider a general integrate-and-fire (IF) neuron driven by asymmetric dichotomous noise. In contrast to the Gaussian white noise usually used in the so-called diffusion approximation, this noise is colored, i.e., it exhibits temporal correlations. We give an analytical expression for the stationary voltage distribution of a neuron receiving such noise and derive recursive relations for the moments of the first passage time density, which allow us to calculate the firing rate and the coefficient of variation of interspike intervals. We study how correlations in the input affect the rate and regularity of firing under variation of the model's parameters for leaky and quadratic IF neurons. Further, we consider the limit of small correlation times and find lowest order corrections to the first passage time moments to be proportional to the square root of the correlation time. We show analytically that to this lowest order, correlations always lead to a decrease in firing rate for a leaky IF neuron. All theoretical expressions are compared to simulations of leaky and quadratic IF neurons.

Extracting surface diffusion coefficients from batch adsorption measurement data: application of the classic Langmuir kinetics model.

PubMed

Chu, Khim Hoong

2017-11-09

Surface diffusion coefficients may be estimated by fitting solutions of a diffusion model to batch kinetic data. For non-linear systems, a numerical solution of the diffusion model's governing equations is generally required. We report here the application of the classic Langmuir kinetics model to extract surface diffusion coefficients from batch kinetic data. The use of the Langmuir kinetics model in lieu of the conventional surface diffusion model allows derivation of an analytical expression. The parameter estimation procedure requires determining the Langmuir rate coefficient from which the pertinent surface diffusion coefficient is calculated. Surface diffusion coefficients within the 10 -9 to 10 -6  cm 2 /s range obtained by fitting the Langmuir kinetics model to experimental kinetic data taken from the literature are found to be consistent with the corresponding values obtained from the traditional surface diffusion model. The virtue of this simplified parameter estimation method is that it reduces the computational complexity as the analytical expression involves only an algebraic equation in closed form which is easily evaluated by spreadsheet computation.

Optimal control, optimization and asymptotic analysis of Purcell's microswimmer model

NASA Astrophysics Data System (ADS)

Wiezel, Oren; Or, Yizhar

2016-11-01

Purcell's swimmer (1977) is a classic model of a three-link microswimmer that moves by performing periodic shape changes. Becker et al. (2003) showed that the swimmer's direction of net motion is reversed upon increasing the stroke amplitude of joint angles. Tam and Hosoi (2007) used numerical optimization in order to find optimal gaits for maximizing either net displacement or Lighthill's energetic efficiency. In our work, we analytically derive leading-order expressions as well as next-order corrections for both net displacement and energetic efficiency of Purcell's microswimmer. Using these expressions enables us to explicitly show the reversal in direction of motion, as well as obtaining an estimate for the optimal stroke amplitude. We also find the optimal swimmer's geometry for maximizing either displacement or energetic efficiency. Additionally, the gait optimization problem is revisited and analytically formulated as an optimal control system with only two state variables, which can be solved using Pontryagin's maximum principle. It can be shown that the optimal solution must follow a "singular arc". Numerical solution of the boundary value problem is obtained, which exactly reproduces Tam and Hosoi's optimal gait.

Mixed analytical-stochastic simulation method for the recovery of a Brownian gradient source from probability fluxes to small windows.

PubMed

Dobramysl, U; Holcman, D

2018-02-15

Is it possible to recover the position of a source from the steady-state fluxes of Brownian particles to small absorbing windows located on the boundary of a domain? To address this question, we develop a numerical procedure to avoid tracking Brownian trajectories in the entire infinite space. Instead, we generate particles near the absorbing windows, computed from the analytical expression of the exit probability. When the Brownian particles are generated by a steady-state gradient at a single point, we compute asymptotically the fluxes to small absorbing holes distributed on the boundary of half-space and on a disk in two dimensions, which agree with stochastic simulations. We also derive an expression for the splitting probability between small windows using the matched asymptotic method. Finally, when there are more than two small absorbing windows, we show how to reconstruct the position of the source from the diffusion fluxes. The present approach provides a computational first principle for the mechanism of sensing a gradient of diffusing particles, a ubiquitous problem in cell biology.

Self-similar solutions to isothermal shock problems

NASA Astrophysics Data System (ADS)

Deschner, Stephan C.; Illenseer, Tobias F.; Duschl, Wolfgang J.

We investigate exact solutions for isothermal shock problems in different one-dimensional geometries. These solutions are given as analytical expressions if possible, or are computed using standard numerical methods for solving ordinary differential equations. We test the numerical solutions against the analytical expressions to verify the correctness of all numerical algorithms. We use similarity methods to derive a system of ordinary differential equations (ODE) yielding exact solutions for power law density distributions as initial conditions. Further, the system of ODEs accounts for implosion problems (IP) as well as explosion problems (EP) by changing the initial or boundary conditions, respectively. Taking genuinely isothermal approximations into account leads to additional insights of EPs in contrast to earlier models. We neglect a constant initial energy contribution but introduce a parameter to adjust the initial mass distribution of the system. Moreover, we show that due to this parameter a constant initial density is not allowed for isothermal EPs. Reasonable restrictions for this parameter are given. Both, the (genuinely) isothermal implosion as well as the explosion problem are solved for the first time.

Residual zonal flows in tokamaks and stellarators at arbitrary wavelengths

NASA Astrophysics Data System (ADS)

Monreal, Pedro; Calvo, Iván; Sánchez, Edilberto; Parra, Félix I.; Bustos, Andrés; Könies, Axel; Kleiber, Ralf; Görler, Tobias

2016-04-01

In the linear collisionless limit, a zonal potential perturbation in a toroidal plasma relaxes, in general, to a non-zero residual value. Expressions for the residual value in tokamak and stellarator geometries, and for arbitrary wavelengths, are derived. These expressions involve averages over the lowest order particle trajectories, that typically cannot be evaluated analytically. In this work, an efficient numerical method for the evaluation of such expressions is reported. It is shown that this method is faster than direct gyrokinetic simulations performed with the Gene and EUTERPE codes. Calculations of the residual value in stellarators are provided for much shorter wavelengths than previously available in the literature. Electrons must be treated kinetically in stellarators because, unlike in tokamaks, kinetic electrons modify the residual value even at long wavelengths. This effect, that had already been predicted theoretically, is confirmed by gyrokinetic simulations.

When gluons go odd: how classical gluon fields generate odd azimuthal harmonics for the two-gluon correlation function in high-energy collisions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kovchegov, Yuri V.; Skokov, Vladimir V.

We show that, in the saturation/Color Glass Condensate framework, odd azimuthal harmonics of the two-gluon correlation function with a long-range separation in rapidity are generated by the higher-order saturation corrections in the interactions with the projectile and the target. At the very least, the odd harmonics require three scatterings in the projectile and three scatterings in the target. We derive the leading-order expression for the two-gluon production cross section which generates odd harmonics: the expression includes all-order interactions with the target and three interactions with the projectile. Here, we evaluate the obtained expression both analytically and numerically, confirming that themore » odd-harmonics contribution to the two-gluon production in the saturation framework is non-zero.« less

When gluons go odd: how classical gluon fields generate odd azimuthal harmonics for the two-gluon correlation function in high-energy collisions

DOE PAGES

Kovchegov, Yuri V.; Skokov, Vladimir V.

2018-04-30

We show that, in the saturation/Color Glass Condensate framework, odd azimuthal harmonics of the two-gluon correlation function with a long-range separation in rapidity are generated by the higher-order saturation corrections in the interactions with the projectile and the target. At the very least, the odd harmonics require three scatterings in the projectile and three scatterings in the target. We derive the leading-order expression for the two-gluon production cross section which generates odd harmonics: the expression includes all-order interactions with the target and three interactions with the projectile. Here, we evaluate the obtained expression both analytically and numerically, confirming that themore » odd-harmonics contribution to the two-gluon production in the saturation framework is non-zero.« less

Modeling and evaluating the performance of Brillouin distributed optical fiber sensors.

PubMed

Soto, Marcelo A; Thévenaz, Luc

2013-12-16

A thorough analysis of the key factors impacting on the performance of Brillouin distributed optical fiber sensors is presented. An analytical expression is derived to estimate the error on the determination of the Brillouin peak gain frequency, based for the first time on real experimental conditions. This expression is experimentally validated, and describes how this frequency uncertainty depends on measurement parameters, such as Brillouin gain linewidth, frequency scanning step and signal-to-noise ratio. Based on the model leading to this expression and considering the limitations imposed by nonlinear effects and pump depletion, a figure-of-merit is proposed to fairly compare the performance of Brillouin distributed sensing systems. This figure-of-merit offers to the research community and to potential users the possibility to evaluate with an objective metric the real performance gain resulting from any proposed configuration.

Recurrences and explicit formulae for the expansion and connection coefficients in series of Bessel polynomials

NASA Astrophysics Data System (ADS)

Doha, E. H.; Ahmed, H. M.

2004-08-01

A formula expressing explicitly the derivatives of Bessel polynomials of any degree and for any order in terms of the Bessel polynomials themselves is proved. Another explicit formula, which expresses the Bessel expansion coefficients of a general-order derivative of an infinitely differentiable function in terms of its original Bessel coefficients, is also given. A formula for the Bessel coefficients of the moments of one single Bessel polynomial of certain degree is proved. A formula for the Bessel coefficients of the moments of a general-order derivative of an infinitely differentiable function in terms of its Bessel coefficients is also obtained. Application of these formulae for solving ordinary differential equations with varying coefficients, by reducing them to recurrence relations in the expansion coefficients of the solution, is explained. An algebraic symbolic approach (using Mathematica) in order to build and solve recursively for the connection coefficients between Bessel-Bessel polynomials is described. An explicit formula for these coefficients between Jacobi and Bessel polynomials is given, of which the ultraspherical polynomial and its consequences are important special cases. Two analytical formulae for the connection coefficients between Laguerre-Bessel and Hermite-Bessel are also developed.

Mechanism of formation of the response of a hydrogen gas sensor based on a silicon MOS diode

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gaman, V. I.; Balyuba, V. I.; Gritsyk, V. Yu.

2008-03-15

Experimental data on the dependence of the flat-band voltage and relaxation time for the capacitance of the space-charge region in an MOS diode (Pd-SiO{sub 2}-n-Si) on the hydrogen concentration in a hydrogen/air gaseous mixture are discussed. It is assumed that variation in the flat-band voltage U{sub fb} in an MOS structure with the thickness d = 369 nm subjected to a hydrogen/air gaseous mixture can be accounted for by the formation of dipoles in the Pd-SiO{sub 2} gap due to polarization of hydrogen atoms (H{sub a}). An analytical expression describing the dependence of variation in the flat-band voltage {Delta}U{sub fb}more » on the hydrogen concentration n{sub H{sub 2}} was derived. In MOS structures with d {<=} 4 nm (or MOS diodes), the value of {Delta}U{sub fb} is mainly controlled by passivation of the centers responsible for the presence of the surface acceptor-type centers at the SiO{sub 2}-n-Si interface by hydrogen atoms. Analytical expressions describing the dependences of {Delta}U{sub fb} and the capacitance relaxation time in the space-charge region on n{sub H{sub 2}} are derived. The values of the density of adsorption centers and the adsorption heat for hydrogen atoms at the Pd-SiO{sub 2} and SiO{sub 2}-n-Si interfaces are found.« less

Mechanism of formation of the response of a hydrogen gas sensor based on a silicon MOS diode

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gaman, V. I.; Balyuba, V. I.; Gritsyk, V. Yu.

2008-03-15

Experimental data on the dependence of the flat-band voltage and relaxation time for the capacitance of the space-charge region in an MOS diode (Pd-SiO{sub 2}-n-Si) on the hydrogen concentration in a hydrogen/air gaseous mixture are discussed. It is assumed that variation in the flat-band voltage U{sub fb} in an MOS structure with the thickness d = 369 nm subjected to a hydrogen/air gaseous mixture can be accounted for by the formation of dipoles in the Pd-SiO{sub 2} gap due to polarization of hydrogen atoms (H{sub a}). An analytical expression describing the dependence of variation in the flat-band voltage {delta}U{sub fb}more » on the hydrogen concentration n{sub H2} was derived. In MOS structures with d {<=} 4 nm (or MOS diodes), the value of {delta}U{sub fb} is mainly controlled by passivation of the centers responsible for the presence of the surface acceptor-type centers at the SiO{sub 2}-n-Si interface by hydrogen atoms. Analytical expressions describing the dependences of {delta}U{sub fb} and the capacitance relaxation time in the space-charge region on n{sub H2} are derived. The values of the density of adsorption centers and the adsorption heat for hydrogen atoms at the Pd-SiO{sub 2} and SiO{sub 2}-n-Si interfaces are found.« less

The MGGB equation-of-state for multifield applications: a numerical recipe for analytic expression of sesame EOS data

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kashiwa, B. A.

2010-12-01

Abstract A thermodynamically consistent and fully general equation–of– state (EOS) for multifield applications is described. EOS functions are derived from a Helmholtz free energy expressed as the sum of thermal (fluctuational) and collisional (condensed–phase) contributions; thus the free energy is of the Mie–Gr¨uneisen1 form. The phase–coexistence region is defined using a parameterized saturation curve by extending the form introduced by Guggenheim,2 which scales the curve relative to conditions at the critical point. We use the zero–temperature condensed–phase contribution developed by Barnes,3 which extends the Thomas–Fermi–Dirac equation to zero pressure. Thus, the functional form of the EOS could be called MGGBmore » (for Mie– Gr¨uneisen–Guggenheim–Barnes). Substance–specific parameters are obtained by fitting the low–density energy to data from the Sesame4 library; fitting the zero–temperature pressure to the Sesame cold curve; and fitting the saturation curve and latent heat to laboratory data,5 if available. When suitable coexistence data, or Sesame data, are not available, then we apply the Principle of Corresponding States.2 Thus MGGB can be thought of as a numerical recipe for rendering the tabular Sesame EOS data in an analytic form that includes a proper coexistence region, and which permits the accurate calculation of derivatives associated with compressibility, expansivity, Joule coefficient, and specific heat, all of which are required for multifield applications. 1« less

Is the Jeffreys' scale a reliable tool for Bayesian model comparison in cosmology?

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nesseris, Savvas; García-Bellido, Juan, E-mail: savvas.nesseris@uam.es, E-mail: juan.garciabellido@uam.es

2013-08-01

We are entering an era where progress in cosmology is driven by data, and alternative models will have to be compared and ruled out according to some consistent criterium. The most conservative and widely used approach is Bayesian model comparison. In this paper we explicitly calculate the Bayes factors for all models that are linear with respect to their parameters. We do this in order to test the so called Jeffreys' scale and determine analytically how accurate its predictions are in a simple case where we fully understand and can calculate everything analytically. We also discuss the case of nestedmore » models, e.g. one with M{sub 1} and another with M{sub 2} superset of M{sub 1} parameters and we derive analytic expressions for both the Bayes factor and the figure of Merit, defined as the inverse area of the model parameter's confidence contours. With all this machinery and the use of an explicit example we demonstrate that the threshold nature of Jeffreys' scale is not a ''one size fits all'' reliable tool for model comparison and that it may lead to biased conclusions. Furthermore, we discuss the importance of choosing the right basis in the context of models that are linear with respect to their parameters and how that basis affects the parameter estimation and the derived constraints.« less

Implementing a Matrix-free Analytical Jacobian to Handle Nonlinearities in Models of 3D Lithospheric Deformation

NASA Astrophysics Data System (ADS)

Kaus, B.; Popov, A.

2015-12-01

The analytical expression for the Jacobian is a key component to achieve fast and robust convergence of the nonlinear Newton-Raphson iterative solver. Accomplishing this task in practice often requires a significant algebraic effort. Therefore it is quite common to use a cheap alternative instead, for example by approximating the Jacobian with a finite difference estimation. Despite its simplicity it is a relatively fragile and unreliable technique that is sensitive to the scaling of the residual and unknowns, as well as to the perturbation parameter selection. Unfortunately no universal rule can be applied to provide both a robust scaling and a perturbation. The approach we use here is to derive the analytical Jacobian for the coupled set of momentum, mass, and energy conservation equations together with the elasto-visco-plastic rheology and a marker in cell/staggered finite difference method. The software project LaMEM (Lithosphere and Mantle Evolution Model) is primarily developed for the thermo-mechanically coupled modeling of the 3D lithospheric deformation. The code is based on a staggered grid finite difference discretization in space, and uses customized scalable solvers form PETSc library to efficiently run on the massively parallel machines (such as IBM Blue Gene/Q). Currently LaMEM relies on the Jacobian-Free Newton-Krylov (JFNK) nonlinear solver, which approximates the Jacobian-vector product using a simple finite difference formula. This approach never requires an assembled Jacobian matrix and uses only the residual computation routine. We use an approximate Jacobian (Picard) matrix to precondition the Krylov solver with the Galerkin geometric multigrid. Because of the inherent problems of the finite difference Jacobian estimation, this approach doesn't always result in stable convergence. In this work we present and discuss a matrix-free technique in which the Jacobian-vector product is replaced by analytically-derived expressions and compare results with those obtained with a finite difference approximation of the Jacobian. This project is funded by ERC Starting Grant 258830 and computer facilities were provided by Jülich supercomputer center (Germany).

An analytical solution to assess the SH seismoelectric response of the vadose zone

NASA Astrophysics Data System (ADS)

Monachesi, L. B.; Zyserman, F. I.; Jouniaux, L.

2018-03-01

We derive an analytical solution of the seismoelectric conversions generated in the vadose zone, when this region is crossed by a pure shear horizontal (SH) wave. Seismoelectric conversions are induced by electrokinetic effects linked to relative motions between fluid and porous media. The considered model assumes a one-dimensional soil constituted by a single layer on top of a half space in contact at the water table, and a shearing force located at the earth's surface as the wave source. The water table is an interface expected to induce a seismoelectric interfacial response (IR). The top layer represents a porous rock which porous space is partially saturated by water and air, while the half-space is completely saturated with water, representing the saturated zone. The analytical expressions for the coseismic fields and the interface responses, both electric and magnetic, are derived by solving Pride's equations with proper boundary conditions. An approximate analytical expression of the solution is also obtained, which is very simple and applicable in a fairly broad set of situations. Hypothetical scenarios are proposed to study and analyse the dependence of the electromagnetic fields on various parameters of the medium. An analysis of the approximate solution is also made together with a comparison to the exact solution. The main result of the present analysis is that the amplitude of the interface response generated at the water table is found to be proportional to the jump in the electric current density, which in turn depends on the saturation contrast, poro-mechanical and electrical properties of the medium and on the amplitude of the solid displacement produced by the source. This result is in agreement with the one numerically obtained by the authors, which has been published in a recent work. We also predict the existence of an interface response located at the surface, and that the electric interface response is several orders of magnitude bigger than the electric coseismic field, whereas it is the opposite using compressional waves as shown by theoretical and experimental results. This fact should encourage the performance of field and laboratory tests to check the viability of SHTE seismoelectrics as a near surface prospecting/monitoring tool.

An analytical solution to assess the SH seismoelectric response of the vadose zone

NASA Astrophysics Data System (ADS)

Monachesi, L. B.; Zyserman, F. I.; Jouniaux, L.

2018-06-01

We derive an analytical solution of the seismoelectric conversions generated in the vadose zone, when this region is crossed by a pure shear horizontal (SH) wave. Seismoelectric conversions are induced by electrokinetic effects linked to relative motions between fluid and porous media. The considered model assumes a 1D soil constituted by a single layer on top of a half-space in contact at the water table, and a shearing force located at the earth's surface as the wave source. The water table is an interface expected to induce a seismoelectric interfacial response (IR). The top layer represents a porous rock in which porous space is partially saturated by water and air, while the half-space is completely saturated with water, representing the saturated zone. The analytical expressions for the coseismic fields and the interface responses, both electric and magnetic, are derived by solving Pride's equations with proper boundary conditions. An approximate analytical expression of the solution is also obtained, which is very simple and applicable in a fairly broad set of situations. Hypothetical scenarios are proposed to study and analyse the dependence of the electromagnetic fields on various parameters of the medium. An analysis of the approximate solution is also made together with a comparison to the exact solution. The main result of the present analysis is that the amplitude of the interface response generated at the water table is found to be proportional to the jump in the electric current density, which in turn depends on the saturation contrast, poro-mechanical and electrical properties of the medium and on the amplitude of the solid displacement produced by the source. This result is in agreement with the one numerically obtained by the authors, which has been published in a recent work. We also predict the existence of an interface response located at the surface, and that the electric interface response is several orders of magnitude bigger than the electric coseismic field, whereas it is the opposite using compressional waves as shown by theoretical and experimental results. This fact should encourage the performance of field and laboratory tests to check the viability of SHTE seismoelectrics as a near surface prospecting/monitoring tool.

Frequency response in short thermocouple wires

NASA Technical Reports Server (NTRS)

Forney, L. J.; Meeks, E. L.; Ma, J.; Fralick, G. C.

1992-01-01

Theoretical expressions are derived for the steady state frequency response of a thermocouple wire. In particular, the effects of axial heat conduction are demonstrated for a nonuniform wire with unequal material properties and wire diameters across the junction. The amplitude ratio at low frequency omega approaches 0 agrees with the results of Scadron and Warshawsky (1952) for a steady state temperature distribution. Moreover, the frequency response for a nonuniform wire in the limit of infinite length l approaches infinity is shown to reduce to a simple expression that is analogous to the classic first order solution for a thermocouple wire with uniform properties. Theoretical expressions are also derived for the steady state frequency response of a supported thermocouple wire. In particular, the effects of axial heat conduction are demonstrated for both a supported one material wire and a two material wire with unequal material properties across the junction. For the case of a one material supported wire, an exact solution is derived which compares favorably with an approximate expression that only matches temperatures at the support junction. Moreover, for the case of a two material supported wire, an analytical expression is derived that closely correlates numerical results. Experimental measurements are made for the steady state frequency response of a supported thermocouple wire. In particular, the effects of axial heat conduction are demonstrated for both a supported one material wire (type K) and a two material wire (type T) with unequal material properties across the junction. The data for the amplitude ratio and phase angle are correlated to within 10 pct. with the theoretical predictions of Forney and Fralick (1991). This is accomplished by choosing a natural frequency omega sub n for the wire data to correlate the first order response at large gas temperature frequencies. It is found that a large bead size, however, will increase the amplitude ratio at low frequencies but decrease the natural frequency of the wire. The phase angle data are also distorted for imperfect junctions.

Solvent effects in time-dependent self-consistent field methods. II. Variational formulations and analytical gradients

DOE PAGES

Bjorgaard, J. A.; Velizhanin, K. A.; Tretiak, S.

2015-08-06

This study describes variational energy expressions and analytical excited state energy gradients for time-dependent self-consistent field methods with polarizable solvent effects. Linear response, vertical excitation, and state-specific solventmodels are examined. Enforcing a variational ground stateenergy expression in the state-specific model is found to reduce it to the vertical excitation model. Variational excited state energy expressions are then provided for the linear response and vertical excitation models and analytical gradients are formulated. Using semiempiricalmodel chemistry, the variational expressions are verified by numerical and analytical differentiation with respect to a static external electric field. Lastly, analytical gradients are further tested by performingmore » microcanonical excited state molecular dynamics with p-nitroaniline.« less

Deriving the exact nonadiabatic quantum propagator in the mapping variable representation.

PubMed

Hele, Timothy J H; Ananth, Nandini

2016-12-22

We derive an exact quantum propagator for nonadiabatic dynamics in multi-state systems using the mapping variable representation, where classical-like Cartesian variables are used to represent both continuous nuclear degrees of freedom and discrete electronic states. The resulting Liouvillian is a Moyal series that, when suitably approximated, can allow for the use of classical dynamics to efficiently model large systems. We demonstrate that different truncations of the exact Liouvillian lead to existing approximate semiclassical and mixed quantum-classical methods and we derive an associated error term for each method. Furthermore, by combining the imaginary-time path-integral representation of the Boltzmann operator with the exact Liouvillian, we obtain an analytic expression for thermal quantum real-time correlation functions. These results provide a rigorous theoretical foundation for the development of accurate and efficient classical-like dynamics to compute observables such as electron transfer reaction rates in complex quantized systems.

Nonequilibrium Green's function theory for nonadiabatic effects in quantum electron transport

NASA Astrophysics Data System (ADS)

Kershaw, Vincent F.; Kosov, Daniel S.

2017-12-01

We develop nonequilibrium Green's function-based transport theory, which includes effects of nonadiabatic nuclear motion in the calculation of the electric current in molecular junctions. Our approach is based on the separation of slow and fast time scales in the equations of motion for Green's functions by means of the Wigner representation. Time derivatives with respect to central time serve as a small parameter in the perturbative expansion enabling the computation of nonadiabatic corrections to molecular Green's functions. Consequently, we produce a series of analytic expressions for non-adiabatic electronic Green's functions (up to the second order in the central time derivatives), which depend not solely on the instantaneous molecular geometry but likewise on nuclear velocities and accelerations. An extended formula for electric current is derived which accounts for the non-adiabatic corrections. This theory is concisely illustrated by the calculations on a model molecular junction.

Nonequilibrium Green's function theory for nonadiabatic effects in quantum electron transport.

PubMed

Kershaw, Vincent F; Kosov, Daniel S

2017-12-14

We develop nonequilibrium Green's function-based transport theory, which includes effects of nonadiabatic nuclear motion in the calculation of the electric current in molecular junctions. Our approach is based on the separation of slow and fast time scales in the equations of motion for Green's functions by means of the Wigner representation. Time derivatives with respect to central time serve as a small parameter in the perturbative expansion enabling the computation of nonadiabatic corrections to molecular Green's functions. Consequently, we produce a series of analytic expressions for non-adiabatic electronic Green's functions (up to the second order in the central time derivatives), which depend not solely on the instantaneous molecular geometry but likewise on nuclear velocities and accelerations. An extended formula for electric current is derived which accounts for the non-adiabatic corrections. This theory is concisely illustrated by the calculations on a model molecular junction.

Exact relations between homoclinic and periodic orbit actions in chaotic systems

NASA Astrophysics Data System (ADS)

Li, Jizhou; Tomsovic, Steven

2018-02-01

Homoclinic and unstable periodic orbits in chaotic systems play central roles in various semiclassical sum rules. The interferences between terms are governed by the action functions and Maslov indices. In this article, we identify geometric relations between homoclinic and unstable periodic orbits, and derive exact formulas expressing the periodic orbit classical actions in terms of corresponding homoclinic orbit actions plus certain phase space areas. The exact relations provide a basis for approximations of the periodic orbit actions as action differences between homoclinic orbits with well-estimated errors. This enables an explicit study of relations between periodic orbits, which results in an analytic expression for the action differences between long periodic orbits and their shadowing decomposed orbits in the cycle expansion.

Intrinsic optical conductivity of a {{\\rm{C}}}_{2v} symmetric topological insulator

NASA Astrophysics Data System (ADS)

Sengupta, Parijat; Matsubara, Masahiko; Bellotti, Enrico; Shi, Junxia

2017-07-01

In this work we analytically investigate the longitudinal optical conductivity of the {{{C}}}2v symmetric topological insulator. The conductivity expressions at T = 0 are derived using the Kubo formula and expressed as a function of the ratio of the Dresselhaus and Rashba parameters that characterize the low-energy Hamiltonian. We find that the longitudinal inter-band conductivity vanishes when Dresselhaus and Rashba parameters are equal in strength, also called the persistent spin helix state. The calculations are extended to obtain the frequency-dependent real and imaginary components of the optical conductivity for the topological Kondo insulator SmB6 which exhibits {{{C}}}2v symmetric and anisotropic Dirac cones hosting topological states at \\overline{X} point on the surface Brillouin zone.

Rapid calculation of acoustic fields from arbitrary continuous-wave sources.

PubMed

Treeby, Bradley E; Budisky, Jakub; Wise, Elliott S; Jaros, Jiri; Cox, B T

2018-01-01

A Green's function solution is derived for calculating the acoustic field generated by phased array transducers of arbitrary shape when driven by a single frequency continuous wave excitation with spatially varying amplitude and phase. The solution is based on the Green's function for the homogeneous wave equation expressed in the spatial frequency domain or k-space. The temporal convolution integral is solved analytically, and the remaining integrals are expressed in the form of the spatial Fourier transform. This allows the acoustic pressure for all spatial positions to be calculated in a single step using two fast Fourier transforms. The model is demonstrated through several numerical examples, including single element rectangular and spherically focused bowl transducers, and multi-element linear and hemispherical arrays.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ayala, Alejandro; Hentschinski, Martin; Jalilian-Marian, Jamal

Azimuthal angular correlations between produced hadrons/jets in high energy collisions are a sensitive probe of the dynamics of QCD at small x. Here we derive the triple differential cross section for inclusive production of 3 polarized partons in DIS at small x using the spinor helicity formalism. The target proton or nucleus is described using the Color Glass Condensate (CGC) formalism. The resulting expressions are used to study azimuthal angular correlations between produced partons in order to probe the gluon structure of the target hadron or nucleus. Finally, our analytic expressions can also be used to calculate the real partmore » of the Next to Leading Order (NLO) corrections to di-hadron production in DIS by integrating out one of the three final state partons.« less

Dependence of structure factor and correlation energy on the width of electron wires

NASA Astrophysics Data System (ADS)

Ashokan, Vinod; Bala, Renu; Morawetz, Klaus; Pathak, Kare Narain

2018-02-01

The structure factor and correlation energy of a quantum wire of thickness b ≪ a B are studied in random phase approximation (RPA) and for the less investigated region r s < 1. Using the single-loop approximation, analytical expressions of the structure factor are obtained. The exact expressions for the exchange energy are also derived for a cylindrical and harmonic wire. The correlation energy in RPA is found to be represented by ɛ c ( b, r s ) = α( r s )/ b + β( r s ) ln( b) + η( r s ), for small b and high densities. For a pragmatic width of the wire, the correlation energy is in agreement with the quantum Monte Carlo simulation data.

Modeling the macrodynamics of international migration: determinants of emigration from Cyprus, 1946-85.

PubMed

Diamantides, N D; Constantinou, S T

1989-07-01

"A model is presented of international migration that is based on the concept of a pool of potential emigrants at the origin created by push-pull forces and by the establishment of information feedback between origin and destination. The forces can be economic, political, or both, and are analytically expressed by the 'mediating factor'. The model is macrodynamic in nature and provides both for the main secular component of the migratory flow and for transient components caused by extraordinary events. The model is expressed in a Bernoulli-type differential equation through which quantitative weights can be derived for each of the operating causes. Out-migration from the Republic of Cyprus is used to test the tenets of the model." excerpt

Exact Analytic Result of Contact Value for the Density in a Modified Poisson-Boltzmann Theory of an Electrical Double Layer.

PubMed

Lou, Ping; Lee, Jin Yong

2009-04-14

For a simple modified Poisson-Boltzmann (SMPB) theory, taking into account the finite ionic size, we have derived the exact analytic expression for the contact values of the difference profile of the counterion and co-ion, as well as of the sum (density) and product profiles, near a charged planar electrode that is immersed in a binary symmetric electrolyte. In the zero ionic size or dilute limit, these contact values reduce to the contact values of the Poisson-Boltzmann (PB) theory. The analytic results of the SMPB theory, for the difference, sum, and product profiles were compared with the results of the Monte-Carlo (MC) simulations [ Bhuiyan, L. B.; Outhwaite, C. W.; Henderson, D. J. Electroanal. Chem. 2007, 607, 54 ; Bhuiyan, L. B.; Henderson, D. J. Chem. Phys. 2008, 128, 117101 ], as well as of the PB theory. In general, the analytic expression of the SMPB theory gives better agreement with the MC data than the PB theory does. For the difference profile, as the electrode charge increases, the result of the PB theory departs from the MC data, but the SMPB theory still reproduces the MC data quite well, which indicates the importance of including steric effects in modeling diffuse layer properties. As for the product profile, (i) it drops to zero as the electrode charge approaches infinity; (ii) the speed of the drop increases with the ionic size, and these behaviors are in contrast with the predictions of the PB theory, where the product is identically 1.

Analytical approximations for effective relative permeability in the capillary limit

NASA Astrophysics Data System (ADS)

Rabinovich, Avinoam; Li, Boxiao; Durlofsky, Louis J.

2016-10-01

We present an analytical method for calculating two-phase effective relative permeability, krjeff, where j designates phase (here CO2 and water), under steady state and capillary-limit assumptions. These effective relative permeabilities may be applied in experimental settings and for upscaling in the context of numerical flow simulations, e.g., for CO2 storage. An exact solution for effective absolute permeability, keff, in two-dimensional log-normally distributed isotropic permeability (k) fields is the geometric mean. We show that this does not hold for krjeff since log normality is not maintained in the capillary-limit phase permeability field (Kj=k·krj) when capillary pressure, and thus the saturation field, is varied. Nevertheless, the geometric mean is still shown to be suitable for approximating krjeff when the variance of ln⁡k is low. For high-variance cases, we apply a correction to the geometric average gas effective relative permeability using a Winsorized mean, which neglects large and small Kj values symmetrically. The analytical method is extended to anisotropically correlated log-normal permeability fields using power law averaging. In these cases, the Winsorized mean treatment is applied to the gas curves for cases described by negative power law exponents (flow across incomplete layers). The accuracy of our analytical expressions for krjeff is demonstrated through extensive numerical tests, using low-variance and high-variance permeability realizations with a range of correlation structures. We also present integral expressions for geometric-mean and power law average krjeff for the systems considered, which enable derivation of closed-form series solutions for krjeff without generating permeability realizations.

Exact statistical results for binary mixing and reaction in variable density turbulence

NASA Astrophysics Data System (ADS)

Ristorcelli, J. R.

2017-02-01

We report a number of rigorous statistical results on binary active scalar mixing in variable density turbulence. The study is motivated by mixing between pure fluids with very different densities and whose density intensity is of order unity. Our primary focus is the derivation of exact mathematical results for mixing in variable density turbulence and we do point out the potential fields of application of the results. A binary one step reaction is invoked to derive a metric to asses the state of mixing. The mean reaction rate in variable density turbulent mixing can be expressed, in closed form, using the first order Favre mean variables and the Reynolds averaged density variance, ⟨ρ2⟩ . We show that the normalized density variance, ⟨ρ2⟩ , reflects the reduction of the reaction due to mixing and is a mix metric. The result is mathematically rigorous. The result is the variable density analog, the normalized mass fraction variance ⟨c2⟩ used in constant density turbulent mixing. As a consequence, we demonstrate that use of the analogous normalized Favre variance of the mass fraction, c″ ⁣2˜ , as a mix metric is not theoretically justified in variable density turbulence. We additionally derive expressions relating various second order moments of the mass fraction, specific volume, and density fields. The central role of the density specific volume covariance ⟨ρ v ⟩ is highlighted; it is a key quantity with considerable dynamical significance linking various second order statistics. For laboratory experiments, we have developed exact relations between the Reynolds scalar variance ⟨c2⟩ its Favre analog c″ ⁣2˜ , and various second moments including ⟨ρ v ⟩ . For moment closure models that evolve ⟨ρ v ⟩ and not ⟨ρ2⟩ , we provide a novel expression for ⟨ρ2⟩ in terms of a rational function of ⟨ρ v ⟩ that avoids recourse to Taylor series methods (which do not converge for large density differences). We have derived analytic results relating several other second and third order moments and see coupling between odd and even order moments demonstrating a natural and inherent skewness in the mixing in variable density turbulence. The analytic results have applications in the areas of isothermal material mixing, isobaric thermal mixing, and simple chemical reaction (in progress variable formulation).

Analytical approach for the fractional differential equations by using the extended tanh method

NASA Astrophysics Data System (ADS)

Pandir, Yusuf; Yildirim, Ayse

2018-07-01

In this study, we consider analytical solutions of space-time fractional derivative foam drainage equation, the nonlinear Korteweg-de Vries equation with time and space-fractional derivatives and time-fractional reaction-diffusion equation by using the extended tanh method. The fractional derivatives are defined in the modified Riemann-Liouville context. As a result, various exact analytical solutions consisting of trigonometric function solutions, kink-shaped soliton solutions and new exact solitary wave solutions are obtained.

Analytic reflected light curves for exoplanets

NASA Astrophysics Data System (ADS)

Haggard, Hal M.; Cowan, Nicolas B.

2018-07-01

The disc-integrated reflected brightness of an exoplanet changes as a function of time due to orbital and rotational motions coupled with an inhomogeneous albedo map. We have previously derived analytic reflected light curves for spherical harmonic albedo maps in the special case of a synchronously rotating planet on an edge-on orbit (Cowan, Fuentes & Haggard). In this paper, we present analytic reflected light curves for the general case of a planet on an inclined orbit, with arbitrary spin period and non-zero obliquity. We do so for two different albedo basis maps: bright points (δ-maps), and spherical harmonics (Y_ l^m-maps). In particular, we use Wigner D-matrices to express an harmonic light curve for an arbitrary viewing geometry as a non-linear combination of harmonic light curves for the simpler edge-on, synchronously rotating geometry. These solutions will enable future exploration of the degeneracies and information content of reflected light curves, as well as fast calculation of light curves for mapping exoplanets based on time-resolved photometry. To these ends, we make available Exoplanet Analytic Reflected Lightcurves, a simple open-source code that allows rapid computation of reflected light curves.

Quantifying entanglement in two-mode Gaussian states

NASA Astrophysics Data System (ADS)

Tserkis, Spyros; Ralph, Timothy C.

2017-12-01

Entangled two-mode Gaussian states are a key resource for quantum information technologies such as teleportation, quantum cryptography, and quantum computation, so quantification of Gaussian entanglement is an important problem. Entanglement of formation is unanimously considered a proper measure of quantum correlations, but for arbitrary two-mode Gaussian states no analytical form is currently known. In contrast, logarithmic negativity is a measure that is straightforward to calculate and so has been adopted by most researchers, even though it is a less faithful quantifier. In this work, we derive an analytical lower bound for entanglement of formation of generic two-mode Gaussian states, which becomes tight for symmetric states and for states with balanced correlations. We define simple expressions for entanglement of formation in physically relevant situations and use these to illustrate the problematic behavior of logarithmic negativity, which can lead to spurious conclusions.

Induced Eddy Currents in Simple Conductive Geometries: Mathematical Formalism Describes the Excitation of Electrical Eddy Currents in a Time-Varying Magnetic Field

DOE PAGES

Nagel, James R.

2017-12-22

In this paper, a complete mathematical formalism is introduced to describe the excitation of electrical eddy currents due to a time-varying magnetic field. The process works by applying a quasistatic approximation to Ampere's law and then segregating the magnetic field into impressed and induced terms. The result is a nonhomogeneous vector Helmholtz equation that can be analytically solved for many practical geometries. Four demonstration cases are then solved under a constant excitation field over all space—an infinite slab in one dimension, a longitudinal cylinder in two dimensions, a transverse cylinder in two dimensions, and a sphere in three dimensions. Numericalmore » simulations are also performed in parallel with analytic computations, all of which verify the accuracy of the derived expressions.« less

Mutual capture of dipolar molecules at low and very low energies. I. Approximate analytical treatment.

PubMed

Nikitin, E E; Troe, J

2010-09-16

Approximate analytical expressions are derived for the low-energy rate coefficients of capture of two identical dipolar polarizable rigid rotors in their lowest nonresonant (j(1) = 0 and j(2) = 0) and resonant (j(1) = 0,1 and j(2) = 1,0) states. The considered range extends from the quantum, ultralow energy regime, characterized by s-wave capture, to the classical regime described within fly wheel and adiabatic channel approaches, respectively. This is illustrated by the table of contents graphic (available on the Web) that shows the scaled rate coefficients for the mutual capture of rotors in the resonant state versus the reduced wave vector between the Bethe zero-energy (left arrows) and classical high-energy (right arrow) limits for different ratios δ of the dipole-dipole to dispersion interaction.

Mode selection and frequency tuning by injection in pulsed TEA-CO2 lasers

NASA Technical Reports Server (NTRS)

Flamant, P. H.; Menzies, R. T.

1983-01-01

An analytical model characterizing pulsed-TEA-CO2-laser injection locking by tunable CW-laser radiation is presented and used to explore the requirements for SLM pulse generation. Photon-density-rate equations describing the laser mechanism are analyzed in terms of the mode competition between photon densities emitted at two frequencies. The expression derived for pulsed dye lasers is extended to homogeneously broadened CO2 lasers, and locking time is defined as a function of laser parameters. The extent to which injected radiation can be detuned from the CO2 line center and continue to produce SLM pulses is investigated experimentally in terms of the analytical framework. The dependence of locking time on the detuning/pressure-broadened-halfwidth ratio is seen as important for spectroscopic applications requiring tuning within the TEA-laser line-gain bandwidth.

A thermoelastic transversely isotropic thick walled cylinder/disk application: An analytical solution and study

NASA Technical Reports Server (NTRS)

Arnold, S. M.

1989-01-01

A continuum theory is utilized to represent the thermoelastic behavior of a thick walled composite cylinder that can be idealized as transversely isotropic. A multiaxial statement of the constitutive theory employed is presented, as well as the out of the plane of isotropy, plane stress, and plane strain reductions. The derived analytical solution presented is valid for a cylindrical tube or thin disk with a concentric hole, subjected to internal and/or external pressure and a general radial temperature distribution. A specific problem examined is that of a thick walled cylinder subjected to an internal and external pressure loading and a linear radial temperature distribution. The results are expressed in nondimensional form and the effects on the response behavior are examined for various material properties, fiber orientation and types of loadings.

Nonstationary Deformation of an Elastic Layer with Mixed Boundary Conditions

NASA Astrophysics Data System (ADS)

Kubenko, V. D.

2016-11-01

The analytic solution to the plane problem for an elastic layer under a nonstationary surface load is found for mixed boundary conditions: normal stress and tangential displacement are specified on one side of the layer (fourth boundary-value problem of elasticity) and tangential stress and normal displacement are specified on the other side of the layer (second boundary-value problem of elasticity). The Laplace and Fourier integral transforms are applied. The inverse Laplace and Fourier transforms are found exactly using tabulated formulas and convolution theorems for various nonstationary loads. Explicit analytical expressions for stresses and displacements are derived. Loads applied to a constant surface area and to a surface area varying in a prescribed manner are considered. Computations demonstrate the dependence of the normal stress on time and spatial coordinates. Features of wave processes are analyzed

Evaluation of generalized degrees of freedom for sparse estimation by replica method

NASA Astrophysics Data System (ADS)

Sakata, A.

2016-12-01

We develop a method to evaluate the generalized degrees of freedom (GDF) for linear regression with sparse regularization. The GDF is a key factor in model selection, and thus its evaluation is useful in many modelling applications. An analytical expression for the GDF is derived using the replica method in the large-system-size limit with random Gaussian predictors. The resulting formula has a universal form that is independent of the type of regularization, providing us with a simple interpretation. Within the framework of replica symmetric (RS) analysis, GDF has a physical meaning as the effective fraction of non-zero components. The validity of our method in the RS phase is supported by the consistency of our results with previous mathematical results. The analytical results in the RS phase are calculated numerically using the belief propagation algorithm.

Induced Eddy Currents in Simple Conductive Geometries: Mathematical Formalism Describes the Excitation of Electrical Eddy Currents in a Time-Varying Magnetic Field

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nagel, James R.

In this paper, a complete mathematical formalism is introduced to describe the excitation of electrical eddy currents due to a time-varying magnetic field. The process works by applying a quasistatic approximation to Ampere's law and then segregating the magnetic field into impressed and induced terms. The result is a nonhomogeneous vector Helmholtz equation that can be analytically solved for many practical geometries. Four demonstration cases are then solved under a constant excitation field over all space—an infinite slab in one dimension, a longitudinal cylinder in two dimensions, a transverse cylinder in two dimensions, and a sphere in three dimensions. Numericalmore » simulations are also performed in parallel with analytic computations, all of which verify the accuracy of the derived expressions.« less

Climatic influence of background and volcanic stratosphere aerosol models

NASA Technical Reports Server (NTRS)

Deschamps, P. Y.; Herman, M.; Lenoble, J.; Tanre, D.

1982-01-01

A simple modelization of the earth atmosphere system including tropospheric and stratospheric aerosols has been derived and tested. Analytical expressions are obtained for the albedo variation due to a thin stratospheric aerosol layer. Also outlined are the physical procedures and the respective influence of the main parameters: aerosol optical thickness, single scattering albedo and asymmetry factor, and sublayer albedo. The method is applied to compute the variation of the zonal and planetary albedos due to a stratospheric layer of background H2SO4 particles and of volcanic ash.

The theory of the gravitational potential applied to orbit prediction

NASA Technical Reports Server (NTRS)

Kirkpatrick, J. C.

1976-01-01

A complete derivation of the geopotential function and its gradient is presented. Also included is the transformation of Laplace's equation from Cartesian to spherical coordinates. The analytic solution to Laplace's equation is obtained from the transformed version, in the classical manner of separating the variables. A cursory introduction to the method devised by Pines, using direction cosines to express the orientation of a point in space, is presented together with sample computer program listings for computing the geopotential function and the components of its gradient. The use of the geopotential function is illustrated.

A second-order frequency-aided digital phase-locked loop for Doppler rate tracking

NASA Astrophysics Data System (ADS)

Chie, C. M.

1980-08-01

A second-order digital phase-locked loop (DPLL) has a finite lock range which is a function of the frequency of the incoming signal to be tracked. For this reason, it is not capable of tracking an input with Doppler rate for an indefinite period of time. In this correspondence, an analytical expression for the hold-in time is derived. In addition, an all-digital scheme to alleviate this problem is proposed based on the information obtained from estimating the input signal frequency.

Imploding spherical and cylindrical shocks

NASA Astrophysics Data System (ADS)

Yousaf, M.

1986-03-01

In this paper it is shown that the value of the similarity exponent α derived analytically by Fujimoto and Mishkin [J. Fluid Mech. 89, 61 (1978); Phys. Fluids 21, 1933 (1978)] is exactly the same as that found by Stanyukovich [Unsteady Motion of Continuous Media, (Academic, New York, 1960)]. Since the result found by Stanyukovich is an approximation to α, Fujimoto and Mishkin's claim to have an exact expression of α is false. The two methods are outlined and Stanyukovich's result is simplified to show its equivalence to the work of Fujimoto and Mishkin.

Fluctuation theorem for Hamiltonian Systems: Le Chatelier's principle

NASA Astrophysics Data System (ADS)

Evans, Denis J.; Searles, Debra J.; Mittag, Emil

2001-05-01

For thermostated dissipative systems, the fluctuation theorem gives an analytical expression for the ratio of probabilities that the time-averaged entropy production in a finite system observed for a finite time takes on a specified value compared to the negative of that value. In the past, it has been generally thought that the presence of some thermostating mechanism was an essential component of any system that satisfies a fluctuation theorem. In the present paper, we point out that a fluctuation theorem can be derived for purely Hamiltonian systems, with or without applied dissipative fields.

From non-trivial geometries to power spectra and vice versa

NASA Astrophysics Data System (ADS)

Brooker, D. J.; Tsamis, N. C.; Woodard, R. P.

2018-04-01

We review a recent formalism which derives the functional forms of the primordial—tensor and scalar—power spectra of scalar potential inflationary models. The formalism incorporates the case of geometries with non-constant first slow-roll parameter. Analytic expressions for the power spectra are given that explicitly display the dependence on the geometric properties of the background. Moreover, we present the full algorithm for using our formalism, to reconstruct the model from the observed power spectra. Our techniques are applied to models possessing "features" in their potential with excellent agreement.

Leading temperature dependence of the conductance in Kondo-correlated quantum dots.

PubMed

Aligia, A A

2018-04-18

Using renormalized perturbation theory in the Coulomb repulsion, we derive an analytical expression for the leading term in the temperature dependence of the conductance through a quantum dot described by the impurity Anderson model, in terms of the renormalized parameters of the model. Taking these parameters from the literature, we compare the results with published ones calculated using the numerical renormalization group obtaining a very good agreement. The approach is superior to alternative perturbative treatments. We compare in particular to the results of a simple interpolative perturbation approach.

Quantitative analysis of random ameboid motion

NASA Astrophysics Data System (ADS)

Bödeker, H. U.; Beta, C.; Frank, T. D.; Bodenschatz, E.

2010-04-01

We quantify random migration of the social ameba Dictyostelium discoideum. We demonstrate that the statistics of cell motion can be described by an underlying Langevin-type stochastic differential equation. An analytic expression for the velocity distribution function is derived. The separation into deterministic and stochastic parts of the movement shows that the cells undergo a damped motion with multiplicative noise. Both contributions to the dynamics display a distinct response to external physiological stimuli. The deterministic component depends on the developmental state and ambient levels of signaling substances, while the stochastic part does not.

Duality and symmetry lost in solid mechanics

NASA Astrophysics Data System (ADS)

Bui, Huy Duong

2008-01-01

Some conservation laws in Solids and Fracture Mechanics present a lack of symmetry between kinematic and dynamic variables. It is shown that Duality is the right tool to re-establish the symmetry between equations and variables and to provide conservation laws of the pure divergence type which provide true path independent integrals. The loss of symmetry of some energetic expressions is exploited to derive a new method for solving some inverse problems. In particular, the earthquake inverse problem is solved analytically. To cite this article: H.D. Bui, C. R. Mecanique 336 (2008).

Elastic model for crimped collagen fibrils

NASA Technical Reports Server (NTRS)

Freed, Alan D.; Doehring, Todd C.

2005-01-01

A physiologic constitutive expression is presented in algorithmic format for the nonlinear elastic response of wavy collagen fibrils found in soft connective tissues. The model is based on the observation that crimped fibrils in a fascicle have a three-dimensional structure at the micron scale that we approximate as a helical spring. The symmetry of this wave form allows the force/displacement relationship derived from Castigliano's theorem to be solved in closed form: all integrals become analytic. Model predictions are in good agreement with experimental observations for mitral-valve chordae tendinece.

Tunneling conductance of amine-linked alkyl chains.

PubMed

Prodan, Emil; Car, Roberto

2008-06-01

The tunneling transport theory developed in ref 9 (Phys. Rev. B 2007, 76, 115102) is applied to molecular devices made of alkyl chains linked to gold electrodes via amine groups. Using the analytic expression of the tunneling conductance derived in our previous work, we identify the key physical quantities that characterize the conductance of these devices. By investigating the transport characteristics of three devices, containing four, six, and eight methyl groups, we extract the dependence of the tunneling conductance on the chain's length, which is an exponential decay law in agreement with recent experimental data.

On the dispersion characteristics of extraordinary mode in a relativistic fully degenerate electron plasma

NASA Astrophysics Data System (ADS)

Noureen, S.; Abbas, G.; Sarfraz, M.

2018-01-01

The study of relativistic degenerate plasmas is important in many astrophysical and laboratory environments. Using linearized relativistic Vlasov-Maxwell equations, a generalized expression for the plasma conductivity tensor is derived. Employing Fermi-Dirac distribution at zero temperature, the dispersion relation of the extraordinary mode in a relativistic degenerate electron plasma is investigated. The propagation characteristics are examined in different relativistic density ranges. The shifting of cutoff points due to relativistic effects is observed analytically and graphically. Non-relativistic and ultra-relativistic limiting cases are also presented.

Diversity Order Analysis of Dual-Hop Relaying with Partial Relay Selection

NASA Astrophysics Data System (ADS)

Bao, Vo Nguyen Quoc; Kong, Hyung Yun

In this paper, we study the performance of dual hop relaying in which the best relay selected by partial relay selection will help the source-destination link to overcome the channel impairment. Specifically, closed-form expressions for outage probability, symbol error probability and achievable diversity gain are derived using the statistical characteristic of the signal-to-noise ratio. Numerical investigation shows that the system achieves diversity of two regardless of relay number and also confirms the correctness of the analytical results. Furthermore, the performance loss due to partial relay selection is investigated.

Topology of large-scale structure. IV - Topology in two dimensions

NASA Technical Reports Server (NTRS)

Melott, Adrian L.; Cohen, Alexander P.; Hamilton, Andrew J. S.; Gott, J. Richard, III; Weinberg, David H.

1989-01-01

In a recent series of papers, an algorithm was developed for quantitatively measuring the topology of the large-scale structure of the universe and this algorithm was applied to numerical models and to three-dimensional observational data sets. In this paper, it is shown that topological information can be derived from a two-dimensional cross section of a density field, and analytic expressions are given for a Gaussian random field. The application of a two-dimensional numerical algorithm for measuring topology to cross sections of three-dimensional models is demonstrated.

Quantitative evaluation of cross correlation between two finite-length time series with applications to single-molecule FRET.

PubMed

Hanson, Jeffery A; Yang, Haw

2008-11-06

The statistical properties of the cross correlation between two time series has been studied. An analytical expression for the cross correlation function's variance has been derived. On the basis of these results, a statistically robust method has been proposed to detect the existence and determine the direction of cross correlation between two time series. The proposed method has been characterized by computer simulations. Applications to single-molecule fluorescence spectroscopy are discussed. The results may also find immediate applications in fluorescence correlation spectroscopy (FCS) and its variants.

Beam-width spreading of vortex beams in free space

NASA Astrophysics Data System (ADS)

Wang, Weiwei; Li, Jinhong; Duan, Meiling

2018-01-01

Based on the extended Huygens-Fresnel principle and the definition of second-order moments of the Wigner distribution function, the analytical expression for the beam-width spreading of Gaussian Schell-model (GSM) vortex beams in free space are derived, and used to study the influence of beam parameters on the beam-width spreading of GSM vortex beams. With the increment of the propagation distance, the beam-width spreading of GSM vortex beams will increase; the bigger the topological charge, spatial correlation length, wavelength and waist width are, the smaller the beam-width spreading is.

Limit regimes of ice formation in turbulent supercooled water.

PubMed

De Santi, Francesca; Olla, Piero

2017-10-01

A study of ice formation in stationary turbulent conditions is carried out in various limit regimes of crystal growth, supercooling, and ice entrainment at the water surface. Analytical expressions for the temperature, salinity, and ice concentration mean profiles are provided, and the role of fluctuations in ice production is numerically quantified. Lower bounds on the ratio of sensible heat flux to latent heat flux to the atmosphere are derived and their dependence on key parameters such as salt rejection in freezing and ice entrainment in the water column is elucidated.

Explicit robust schemes for implementation of general principal value-based constitutive models

NASA Technical Reports Server (NTRS)

Arnold, S. M.; Saleeb, A. F.; Tan, H. Q.; Zhang, Y.

1993-01-01

The issue of developing effective and robust schemes to implement general hyperelastic constitutive models is addressed. To this end, special purpose functions are used to symbolically derive, evaluate, and automatically generate the associated FORTRAN code for the explicit forms of the corresponding stress function and material tangent stiffness tensors. These explicit forms are valid for the entire deformation range. The analytical form of these explicit expressions is given here for the case in which the strain-energy potential is taken as a nonseparable polynomial function of the principle stretches.

Spectral response analysis of PVDF capacitive sensors

NASA Astrophysics Data System (ADS)

Reyes-Ramírez, B.; García-Segundo, C.; García-Valenzuela, A.

2013-06-01

We investigate the spectral response to ultrasound waves in water of low-noise capacitive sensors based on PVDF polymer piezoelectric films. First, we analyze theoretically the mechanical-to-electrical transduction as a function of the frequency of ultrasonic signals and derive an analytic expression of the sensor's transfer function. Then we present experimental results of the frequency response of a home-made PDVF in water to test signals from 1 to 20 MHz induced by a commercial hydrophone powered by a signal generator and compare with our theoretical model.

Spectral Analysis for Weighted Iterated Triangulations of Graphs

NASA Astrophysics Data System (ADS)

Chen, Yufei; Dai, Meifeng; Wang, Xiaoqian; Sun, Yu; Su, Weiyi

Much information about the structural properties and dynamical aspects of a network is measured by the eigenvalues of its normalized Laplacian matrix. In this paper, we aim to present a first study on the spectra of the normalized Laplacian of weighted iterated triangulations of graphs. We analytically obtain all the eigenvalues, as well as their multiplicities from two successive generations. As an example of application of these results, we then derive closed-form expressions for their multiplicative Kirchhoff index, Kemeny’s constant and number of weighted spanning trees.

Geometry of Thin Nematic Elastomer Sheets

NASA Astrophysics Data System (ADS)

Aharoni, Hillel; Sharon, Eran; Kupferman, Raz

A thin sheet of nematic elastomer attains 3D configurations depending on the nematic director field upon heating. In this talk we describe the intrinsic geometry of such a sheet, and derive an expression for the metric induced by general smooth nematic director fields. Furthermore, we investigate the reverse problem of constructing a director field that induces a specified 2D geometry. We provide an explicit analytical recipe for constructing any surface of revolution using this method. We demonstrate how the design of an arbitrary 2D geometry is accessible using approximate numerical methods.

Truncation of CPC solar collectors and its effect on energy collection

NASA Astrophysics Data System (ADS)

Carvalho, M. J.; Collares-Pereira, M.; Gordon, J. M.; Rabl, A.

1985-01-01

Analytic expressions are derived for the angular acceptance function of two-dimensional compound parabolic concentrator solar collectors (CPC's) of arbitrary degree of truncation. Taking into account the effect of truncation on both optical and thermal losses in real collectors, the increase in monthly and yearly collectible energy is also evaluated. Prior analyses that have ignored the correct behavior of the angular acceptance function at large angles for truncated collectors are shown to be in error by 0-2 percent in calculations of yearly collectible energy for stationary collectors.

Suppression of Beneficial Mutations in Dynamic Microbial Populations

NASA Astrophysics Data System (ADS)

Bittihn, Philip; Hasty, Jeff; Tsimring, Lev S.

2017-01-01

Quantitative predictions for the spread of mutations in bacterial populations are essential to interpret evolution experiments and to improve the stability of synthetic gene circuits. We derive analytical expressions for the suppression factor for beneficial mutations in populations that undergo periodic dilutions, covering arbitrary population sizes, dilution factors, and growth advantages in a single stochastic model. We find that the suppression factor grows with the dilution factor and depends nontrivially on the growth advantage, resulting in the preferential elimination of mutations with certain growth advantages. We confirm our results by extensive numerical simulations.

Field-emission from parabolic tips: Current distributions, the net current, and effective emission area

NASA Astrophysics Data System (ADS)

Biswas, Debabrata

2018-04-01

Field emission from nano-structured emitters primarily takes place from the tips. Using recent results on the variation of the enhancement factor around the apex [Biswas et al., Ultramicroscopy 185, 1-4 (2018)], analytical expressions for the surface distribution of net emitted electrons, as well as the total and normal energy distributions are derived in terms of the apex radius Ra and the local electric field at the apex Ea. Formulae for the net emitted current and effective emission area in terms of these quantities are also obtained.

The time delay in strong gravitational lensing with Gauss-Bonnet correction

DOE Office of Scientific and Technical Information (OSTI.GOV)

Man, Jingyun; Cheng, Hongbo, E-mail: jingyunman@mail.ecust.edu.cn, E-mail: hbcheng@ecust.edu.cn

2014-11-01

The time delay between two relativistic images in the strong gravitational lensing governed by Gauss-Bonnet gravity is studied. We make a complete analytical derivation of the expression of time delay in presence of Gauss-Bonnet coupling. With respect to Schwarzschild, the time delay decreases as a consequence of the shrinking of the photon sphere. As the coupling increases, the second term in the time delay expansion becomes more relevant. Thus time delay in strong limit encodes some new information about geometry in five-dimensional spacetime with Gauss-Bonnet correction.

Dynamics of a Two-Link Vehicle in an L-Shaped Corridor Revisited

NASA Astrophysics Data System (ADS)

Antonyuk, E. Ya.; Zabuga, A. T.

2014-03-01

The kinematics of a two-link mobile robot with three steerable wheels moving in an L-shaped corridor is analyzed. A smooth (with continuous first derivative) path is designed maintaining the optimal maneuverability of the vehicle. The motion of the vehicle along this path is planned. Analytical expressions for the reactions at the contact of the wheels with the ground are given in the general case of motion. The radius of curvature of the programmed path is shown to have a strong influence on the reactions.

Statistics of wormlike chains. II. Phase transition of polymer liquid crystals and its mixture with low molecular weight liquid crystals

NASA Astrophysics Data System (ADS)

Zhang, W. X.; Zhao, S. R.; Sun, C. P.

1997-02-01

A general self-consistent field (SCF) for the mixture of polymer and low molecular weight (LMW) molecules has been derived by variation principle. Considering a Maier-Saupe type of interaction, the analytical expressions of the SCF for polymer liquid crystals (PLCs) and the mixture of PLCs and LMW liquid crystals are obtained, from which the phase behaviors of PLCs as well as the mixture are studied. The theoretical results are in agreement with experimental results by adjusting a parameter.

A continued fraction resummation form of bath relaxation effect in the spin-boson model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gong, Zhihao; Tang, Zhoufei; Wu, Jianlan, E-mail: jianlanwu@zju.edu.cn

2015-02-28

In the spin-boson model, a continued fraction form is proposed to systematically resum high-order quantum kinetic expansion (QKE) rate kernels, accounting for the bath relaxation effect beyond the second-order perturbation. In particular, the analytical expression of the sixth-order QKE rate kernel is derived for resummation. With higher-order correction terms systematically extracted from higher-order rate kernels, the resummed quantum kinetic expansion approach in the continued fraction form extends the Pade approximation and can fully recover the exact quantum dynamics as the expansion order increases.