Visualizing second order tensor fields with hyperstreamlines
NASA Technical Reports Server (NTRS)
Delmarcelle, Thierry; Hesselink, Lambertus
1993-01-01
Hyperstreamlines are a generalization to second order tensor fields of the conventional streamlines used in vector field visualization. As opposed to point icons commonly used in visualizing tensor fields, hyperstreamlines form a continuous representation of the complete tensor information along a three-dimensional path. This technique is useful in visulaizing both symmetric and unsymmetric three-dimensional tensor data. Several examples of tensor field visualization in solid materials and fluid flows are provided.
Visualization of second order tensor fields and matrix data
NASA Technical Reports Server (NTRS)
Delmarcelle, Thierry; Hesselink, Lambertus
1992-01-01
We present a study of the visualization of 3-D second order tensor fields and matrix data. The general problem of visualizing unsymmetric real or complex Hermitian second order tensor fields can be reduced to the simultaneous visualization of a real and symmetric second order tensor field and a real vector field. As opposed to the discrete iconic techniques commonly used in multivariate data visualization, the emphasis is on exploiting the mathematical properties of tensor fields in order to facilitate their visualization and to produce a continuous representation of the data. We focus on interactively sensing and exploring real and symmetric second order tensor data by generalizing the vector notion of streamline to the tensor concept of hyperstreamline. We stress the importance of a structural analysis of the data field analogous to the techniques of vector field topology extraction in order to obtain a unique and objective representation of second order tensor fields.
Deffayet, C.; Deser, S.; Esposito-Farese, G.
2009-09-15
We extend to curved backgrounds all flat-space scalar field models that obey purely second-order equations, while maintaining their second-order dependence on both field and metric. This extension simultaneously restores to second order the, originally higher derivative, stress tensors as well. The process is transparent and uniform for all dimensions.
The most general second-order field equations of bi-scalar-tensor theory in four dimensions
NASA Astrophysics Data System (ADS)
Ohashi, Seiju; Tanahashi, Norihiro; Kobayashi, Tsutomu; Yamaguchi, Masahide
2015-07-01
The Horndeski theory is known as the most general scalar-tensor theory with second-order field equations. In this paper, we explore the bi-scalar extension of the Horndeski theory. Following Horndeski's approach, we determine all the possible terms appearing in the second-order field equations of the bi-scalar-tensor theory. We compare the field equations with those of the generalized multi-Galileons, and confirm that our theory contains new terms that are not included in the latter theory. We also discuss the construction of the Lagrangian leading to our most general field equations.
Spacetime encodings. III. Second order Killing tensors
Brink, Jeandrew
2010-01-15
This paper explores the Petrov type D, stationary axisymmetric vacuum (SAV) spacetimes that were found by Carter to have separable Hamilton-Jacobi equations, and thus admit a second-order Killing tensor. The derivation of the spacetimes presented in this paper borrows from ideas about dynamical systems, and illustrates concepts that can be generalized to higher-order Killing tensors. The relationship between the components of the Killing equations and metric functions are given explicitly. The origin of the four separable coordinate systems found by Carter is explained and classified in terms of the analytic structure associated with the Killing equations. A geometric picture of what the orbital invariants may represent is built. Requiring that a SAV spacetime admits a second-order Killing tensor is very restrictive, selecting very few candidates from the group of all possible SAV spacetimes. This restriction arises due to the fact that the consistency conditions associated with the Killing equations require that the field variables obey a second-order differential equation, as opposed to a fourth-order differential equation that imposes the weaker condition that the spacetime be SAV. This paper introduces ideas that could lead to the explicit computation of more general orbital invariants in the form of higher-order Killing tensors.
Superquadric glyphs for symmetric second-order tensors.
Schultz, Thomas; Kindlmann, Gordon L
2010-01-01
Symmetric second-order tensor fields play a central role in scientific and biomedical studies as well as in image analysis and feature-extraction methods. The utility of displaying tensor field samples has driven the development of visualization techniques that encode the tensor shape and orientation into the geometry of a tensor glyph. With some exceptions, these methods work only for positive-definite tensors (i.e. having positive eigenvalues, such as diffusion tensors). We expand the scope of tensor glyphs to all symmetric second-order tensors in two and three dimensions, gracefully and unambiguously depicting any combination of positive and negative eigenvalues. We generalize a previous method of superquadric glyphs for positive-definite tensors by drawing upon a larger portion of the superquadric shape space, supplemented with a coloring that indicates the quadratic form (including eigenvalue sign). We show that encoding arbitrary eigenvalue magnitudes requires design choices that differ fundamentally from those in previous work on traceless tensors that arise in the study of liquid crystals. Our method starts with a design of 2-D tensor glyphs guided by principles of scale-preservation and symmetry, and creates 3-D glyphs that include the 2-D glyphs in their axis-aligned cross-sections. A key ingredient of our method is a novel way of mapping from the shape space of three-dimensional symmetric second-order tensors to the unit square. We apply our new glyphs to stress tensors from mechanics, geometry tensors and Hessians from image analysis, and rate-of-deformation tensors in computational fluid dynamics. PMID:20975202
Elasto-plastic model with second order defect density tensor
NASA Astrophysics Data System (ADS)
Cleja-Ţigoiu, Sanda
2011-05-01
The paper deals with a second order finite elasto-plastic model, which involves the defect density tensor, as a measure of the extra material defects existing in the damaged microstructure. The material behaviour is described with respect to an anholonomic configuration, which is introduced through the second order plastic deformation, consisting in plastic distortion and plastic connection. The defect density tensor enters the expression of the plastic connection through its gradient and represents a measure of non-metricity. The constitutive and evolution equations are derived to be compatible with the free energy imbalance. The evolution equation for the defect density tensor is non-local and coupled with the plastic distortion.
Disformal invariance of second order scalar-tensor theories: Framing the Horndeski action
NASA Astrophysics Data System (ADS)
Bettoni, Dario; Liberati, Stefano
2013-10-01
The Horndeski action is the most general one involving a metric and a scalar field that leads to second-order field equations in four dimensions. Being the natural extension of the well-known scalar-tensor theories, its structure and properties are worth analyzing along the experience accumulated in the latter context. Here, we argue that disformal transformations play, for the Horndeski theory, a similar role to that of conformal transformations for scalar-tensor theories a là Brans-Dicke. We identify the most general transformation preserving second-order field equations and discuss the issue of viable frames for this kind of theory, in particular, the possibility to cast the action in the so-called Einstein frame. Interestingly, we find that only for a subset of the Horndeski Lagrangian such a frame exists. Finally, we investigate the transformation properties of such frames under field redefinitions and frame transformations and their reciprocal relationship.
Vector and tensor contributions to the curvature perturbation at second order
NASA Astrophysics Data System (ADS)
Carrilho, Pedro; Malik, Karim A.
2016-02-01
We derive the evolution equation for the second order curvature perturbation using standard techniques of cosmological perturbation theory. We do this for different definitions of the gauge invariant curvature perturbation, arising from different splits of the spatial metric, and compare the expressions. The results are valid at all scales and include all contributions from scalar, vector and tensor perturbations, as well as anisotropic stress, with all our results written purely in terms of gauge invariant quantities. Taking the large-scale approximation, we find that a conserved quantity exists only if, in addition to the non-adiabatic pressure, the transverse traceless part of the anisotropic stress tensor is also negligible. We also find that the version of the gauge invariant curvature perturbation which is exactly conserved is the one defined with the determinant of the spatial part of the inverse metric.
Slowly rotating scalar field wormholes: The second order approximation
Kashargin, P. E.; Sushkov, S. V.
2008-09-15
We discuss rotating wormholes in general relativity with a scalar field with negative kinetic energy. To solve the problem, we use the assumption about slow rotation. The role of a small dimensionless parameter plays the ratio of the linear velocity of rotation of the wormhole's throat and the velocity of light. We construct the rotating wormhole solution in the second-order approximation with respect to the small parameter. The analysis shows that the asymptotical mass of the rotating wormhole is greater than that of the nonrotating one, and the null energy condition violation in the rotating wormhole spacetime is weaker than that in the nonrotating one.
One-electron contributions to the g-tensor for second-order Douglas-Kroll-Hess theory
NASA Astrophysics Data System (ADS)
Sandhoefer, B.; Neese, F.
2012-09-01
The electric g-tensor is a central quantity for the interpretation of electron paramagnetic resonance spectra. In this paper, a detailed derivation of the 1-electron contributions to the g-tensor is presented in the framework of linear response theory and the second-order Douglas-Kroll-Hess (DKH) transformation. Importantly, the DKH transformation in the presence of a magnetic field is not unique. Whether or not the magnetic field is included in the required Foldy-Wouthuysen transformation, different transformation matrices and, consequently, Hamiltonians result. In this paper, a detailed comparison of both approaches is presented, paying particular attention to the mathematical properties of the resulting Hamiltonians. In contrast to previous studies that address the g-tensor in the framework of DKH theory, the resulting terms are compared to those of the conventional Pauli theory and are given a physical interpretation. Based on these mathematical and physical arguments, we establish that the proper DKH transformation for systems with constant magnetic fields is based on a gauge-invariant Foldy-Wouthuysen transformation, i.e., a Foldy-Wouthuysen transformation including the magnetic field. Calculations using density functional theory (DFT) are carried out on a set of heavy, diatomic molecules, and a set of transition-metal complexes. Based on these calculations, the performance of the relativistic calculation with and without inclusion of picture-change effects is compared. Additionally, the g-tensor is calculated for the Lanthanide dihydrides. Together with the results from the other two molecular test sets, these calculations serve to quantify the magnitude of picture-change effects and elucidate trends across the periodic table.
One-electron contributions to the g-tensor for second-order Douglas-Kroll-Hess theory.
Sandhoefer, B; Neese, F
2012-09-01
The electric g-tensor is a central quantity for the interpretation of electron paramagnetic resonance spectra. In this paper, a detailed derivation of the 1-electron contributions to the g-tensor is presented in the framework of linear response theory and the second-order Douglas-Kroll-Hess (DKH) transformation. Importantly, the DKH transformation in the presence of a magnetic field is not unique. Whether or not the magnetic field is included in the required Foldy-Wouthuysen transformation, different transformation matrices and, consequently, Hamiltonians result. In this paper, a detailed comparison of both approaches is presented, paying particular attention to the mathematical properties of the resulting Hamiltonians. In contrast to previous studies that address the g-tensor in the framework of DKH theory, the resulting terms are compared to those of the conventional Pauli theory and are given a physical interpretation. Based on these mathematical and physical arguments, we establish that the proper DKH transformation for systems with constant magnetic fields is based on a gauge-invariant Foldy-Wouthuysen transformation, i.e., a Foldy-Wouthuysen transformation including the magnetic field. Calculations using density functional theory (DFT) are carried out on a set of heavy, diatomic molecules, and a set of transition-metal complexes. Based on these calculations, the performance of the relativistic calculation with and without inclusion of picture-change effects is compared. Additionally, the g-tensor is calculated for the Lanthanide dihydrides. Together with the results from the other two molecular test sets, these calculations serve to quantify the magnitude of picture-change effects and elucidate trends across the periodic table. PMID:22957550
NASA Astrophysics Data System (ADS)
Renaux-Petel, Sébastien
2012-02-01
In this short note we explain how to use the linear equations of motion to simplify the third-order action for the cosmological fluctuations. No field redefinition is needed in this exact procedure which considerably limits the range of independent cubic operators, and hence of possible shapes of the primordial bispectrum. We demonstrate this in the context of the most general single-field scalar-tensor theory with second-order equations of motion, whose third-order action has been calculated recently in arXiv:1107.2642 and 1107.3917. In particular, we show that the three cubic operators initially pointed out in these works as new compared to k-inflation can actually be expressed in terms of standard k-inflationary operators.
Perturbations of matter fields in the second-order gauge-invariant cosmological perturbation theory
NASA Astrophysics Data System (ADS)
Nakamura, Kouji
2009-12-01
To show that the general framework of the second-order gauge-invariant perturbation theory developed by K. Nakamura [Prog. Theor. Phys. 110, 723 (2003)PTPKAV0033-068X10.1143/PTP.110.723; Prog. Theor. Phys. 113, 481 (2005)PTPKAV0033-068X10.1143/PTP.113.481] is applicable to a wide class of cosmological situations, some formulas for the perturbations of the matter fields are summarized within the framework of the second-order gauge-invariant cosmological perturbation theory in a four-dimensional homogeneous isotropic universe, which is developed in Prog. Theor. Phys. 117, 17 (2007)PTPKAV0033-068X10.1143/PTP.117.17. We derive the formulas for the perturbations of the energy-momentum tensors and equations of motion for a perfect fluid, an imperfect fluid, and a single scalar field, and show that all equations are derived in terms of gauge-invariant variables without any gauge fixing. Through these formulas, we may say that the decomposition formulas for the perturbations of any tensor field into gauge-invariant and gauge-variant parts, which are proposed in the above papers, are universal.
NASA Astrophysics Data System (ADS)
Delgado Acosta, E. G.; Banda Guzmán, V. M.; Kirchbach, M.
2015-03-01
We propose a general method for the description of arbitrary single spin- j states transforming according to ( j, 0) ⊕ (0, j) carrier spaces of the Lorentz algebra in terms of Lorentz tensors for bosons, and tensor-spinors for fermions, and by means of second-order Lagrangians. The method allows to avoid the cumbersome matrix calculus and higher ∂2 j order wave equations inherent to the Weinberg-Joos approach. We start with reducible Lorentz tensor (tensor-spinor) representation spaces hosting one sole ( j, 0) ⊕ (0, j) irreducible sector and design there a representation reduction algorithm based on one of the Casimir invariants of the Lorentz algebra. This algorithm allows us to separate neatly the pure spin- j sector of interest from the rest, while preserving the separate Lorentz and Dirac indexes. However, the Lorentz invariants are momentum independent and do not provide wave equations. Genuine wave equations are obtained by conditioning the Lorentz tensors under consideration to satisfy the Klein-Gordon equation. In so doing, one always ends up with wave equations and associated Lagrangians that are of second order in the momenta. Specifically, a spin-3/2 particle transforming as (3/2, 0) ⊕ (0, 3/2) is comfortably described by a second-order Lagrangian in the basis of the totally anti-symmetric Lorentz tensor-spinor of second rank, Ψ [ μν]. Moreover, the particle is shown to propagate causally within an electromagnetic background. In our study of (3/2, 0) ⊕ (0, 3/2) as part of Ψ [ μν] we reproduce the electromagnetic multipole moments known from the Weinberg-Joos theory. We also find a Compton differential cross-section that satisfies unitarity in forward direction. The suggested tensor calculus presents itself very computer friendly with respect to the symbolic software FeynCalc.
NASA Astrophysics Data System (ADS)
Hohenstein, Edward G.; Kokkila, Sara I. L.; Parrish, Robert M.; Martínez, Todd J.
2013-03-01
The second-order approximate coupled cluster singles and doubles method (CC2) is a valuable tool in electronic structure theory. Although the density fitting approximation has been successful in extending CC2 to larger molecules, it cannot address the steep O(N^5) scaling with the number of basis functions, N. Here, we introduce the tensor hypercontraction (THC) approximation to CC2 (THC-CC2), which reduces the scaling to O(N^4) and the storage requirements to O(N^2). We present an algorithm to efficiently evaluate the THC-CC2 correlation energy and demonstrate its quartic scaling. This implementation of THC-CC2 uses a grid-based least-squares THC (LS-THC) approximation to the density-fitted electron repulsion integrals. The accuracy of the CC2 correlation energy under these approximations is shown to be suitable for most practical applications.
NASA Astrophysics Data System (ADS)
Sei, Masaki; Nagayama, Kohei; Kajikawa, Kotaro; Ishii, Hisao; Seki, Kazuhiko; Kondo, Katsumi; Matsumoto, Yoshiyasu; Ouchi, Yukio
1998-04-01
We demonstrated full determination of second-order nonlinear susceptibility of a 4‧-n-octyl-4-cyanobiphenyl (8CB) liquid crystal (LC) monolayer adsorbed on a second-harmonic (SH) active polyimide (PI) substrate. In order to separate the SH signal of the LC film from that of the PI film, we adopted an interferometry technique of second-harmonic generation (SHG) using an ultra-thin film local oscillator. We have found a variety of phases in the components of susceptibility: those of χzii and χizi are almost the same but the phase of χzzz differs by 80° from the other two. The phases of the components of the surface susceptibility tensor are not always identical. This fact indicates that the surface SH response is more complicated than what we expected.
The Topology of Symmetric Tensor Fields
NASA Technical Reports Server (NTRS)
Levin, Yingmei; Batra, Rajesh; Hesselink, Lambertus; Levy, Yuval
1997-01-01
Combinatorial topology, also known as "rubber sheet geometry", has extensive applications in geometry and analysis, many of which result from connections with the theory of differential equations. A link between topology and differential equations is vector fields. Recent developments in scientific visualization have shown that vector fields also play an important role in the analysis of second-order tensor fields. A second-order tensor field can be transformed into its eigensystem, namely, eigenvalues and their associated eigenvectors without loss of information content. Eigenvectors behave in a similar fashion to ordinary vectors with even simpler topological structures due to their sign indeterminacy. Incorporating information about eigenvectors and eigenvalues in a display technique known as hyperstreamlines reveals the structure of a tensor field. The simplify and often complex tensor field and to capture its important features, the tensor is decomposed into an isotopic tensor and a deviator. A tensor field and its deviator share the same set of eigenvectors, and therefore they have a similar topological structure. A a deviator determines the properties of a tensor field, while the isotopic part provides a uniform bias. Degenerate points are basic constituents of tensor fields. In 2-D tensor fields, there are only two types of degenerate points; while in 3-D, the degenerate points can be characterized in a Q'-R' plane. Compressible and incompressible flows share similar topological feature due to the similarity of their deviators. In the case of the deformation tensor, the singularities of its deviator represent the area of vortex core in the field. In turbulent flows, the similarities and differences of the topology of the deformation and the Reynolds stress tensors reveal that the basic addie-viscosity assuptions have their validity in turbulence modeling under certain conditions.
Spatial Variances of Wind Fields and Their Relation to Second-Order Structure Functions and Spectra
NASA Astrophysics Data System (ADS)
King, G. P.; Vogelzang, J.; Stoffelen, A.; Portabella, M.
2014-12-01
Kinetic energy variance as a function of spatial scale for wind fields is commonly estimated either using second-order structure functions (in the spatial domain) or by spectral analysis (in the frequency domain). It will be demonstrated that neither spectra nor second-order structure functions offer a good representation of the variance as a function of scale. These difficulties can be circumvented by using a statistical quantity called spatial variance. It combines the advantages of spectral analysis and spatial statistics. In particular, when applied to observations, spatial variances have a clear interpretation and are tolerant for missing data. They can be related to second-order structure functions, both for discrete and continuous data. For data sets without missing points the relation is statistically exact. Spatial variances can also be Fourier transformed to yield a relation with spectra. The flexibility of spatial variances is used to study various sampling strategies, and to compare them with second-order structure functions and spectral variances. It is shown that the spectral sampling strategy is not seriously biased to calm conditions for scatterometer ocean surface vector winds, and that one-fifth of the second-order structure function value is a good proxy for the cumulative variance.
Second-order magnetic critical points at finite magnetic fields: Revisiting Arrott plots
NASA Astrophysics Data System (ADS)
Bustingorry, S.; Pomiro, F.; Aurelio, G.; Curiale, J.
2016-06-01
The so-called Arrott plot, which consists in plotting H /M against M2, with H the applied magnetic field and M the magnetization, is used to extract valuable information in second-order magnetic phase transitions. Besides, it is widely accepted that a negative slope in the Arrott plot is indicative of a first-order magnetic transition. This is known as the Banerjee criterion. In consequence, the zero-field transition temperature T* is reported as the characteristic first-order transition temperature. By carefully analyzing the mean-field Landau model used for studying first-order magnetic transitions, we show in this work that T* corresponds in fact to a triple point where three first-order lines meet. More importantly, this analysis reveals the existence of two symmetrical second-order critical points at finite magnetic field (Tc,±Hc) . We then show that a modified Arrott plot can be used to obtain information about these second-order critical points. To support this idea we analyze experimental data on La2 /3Ca1 /3MnO3 and discuss an estimate for the location of the triple point and the second-order critical points.
Second-order theory for nonlinear dielectric composites incorporating field fluctuations
NASA Astrophysics Data System (ADS)
Ponte Castañeda, P.
2001-12-01
This paper deals with the development of an improved second-order theory for estimating the effective behavior of nonlinear composite dielectrics. The theory makes use of the field fluctuations in the phases of the relevant ``linear comparison composite'' to generate improved Maxwell-Garnett (MGA) and effective-medium (EMA) types of approximations for nonlinear media. Similar to the earlier version of the theory, the resulting MGA and EMA predictions are exact to second-order in the contrast, but-unlike the earlier version-the estimates satisfy all known bounds. In particular, the EMA estimates exhibit a nonlinearity-independent percolation threshold, and critical exponents that are consistent with recently developed bounds on these exponents. In addition, the MGA and EMA estimates are shown to yield reasonable predictions for strongly nonlinear composites with ``threshold-type'' nonlinearities, which are extreme cases where earlier methods have been known to sometimes fail.
Second-Order Perturbation Theory for Generalized Active Space Self-Consistent-Field Wave Functions.
Ma, Dongxia; Li Manni, Giovanni; Olsen, Jeppe; Gagliardi, Laura
2016-07-12
A multireference second-order perturbation theory approach based on the generalized active space self-consistent-field (GASSCF) wave function is presented. Compared with the complete active space (CAS) and restricted active space (RAS) wave functions, GAS wave functions are more flexible and can employ larger active spaces and/or different truncations of the configuration interaction expansion. With GASSCF, one can explore chemical systems that are not affordable with either CASSCF or RASSCF. Perturbation theory to second order on top of GAS wave functions (GASPT2) has been implemented to recover the remaining electron correlation. The method has been benchmarked by computing the chromium dimer ground-state potential energy curve. These calculations show that GASPT2 gives results similar to CASPT2 even with a configuration interaction expansion much smaller than the corresponding CAS expansion. PMID:27276688
Modeling of finite-amplitude sound beams: second order fields generated by a parametric loudspeaker.
Yang, Jun; Sha, Kan; Gan, Woon-Seng; Tian, Jing
2005-04-01
The nonlinear interaction of sound waves in air has been applied to sound reproduction for audio applications. A directional audible sound can be generated by amplitude-modulating the ultrasound carrier with an audio signal, then transmitting it from a parametric loudspeaker. This brings the need of a computationally efficient model to describe the propagation of finite-amplitude sound beams for the system design and optimization. A quasilinear analytical solution capable of fast numerical evaluation is presented for the second-order fields of the sum-, difference-frequency and second harmonic components. It is based on a virtual-complex-source approach, wherein the source field is treated as an aggregation of a set of complex virtual sources located in complex distance, then the corresponding fundamental sound field is reduced to the computation of sums of simple functions by exploiting the integrability of Gaussian functions. By this result, the five-dimensional integral expressions for the second-order sound fields are simplified to one-dimensional integrals. Furthermore, a substantial analytical reduction to sums of single integrals also is derived for an arbitrary source distribution when the basis functions are expressible as a sum of products of trigonometric functions. The validity of the proposed method is confirmed by a comparison of numerical results with experimental data previously published for the rectangular ultrasonic transducer. PMID:16060510
NASA Technical Reports Server (NTRS)
Gary, S. P.; Tokar, R. L.
1985-01-01
The present investigation is concerned with the application of a second-order theory for electromagnetic instabilities in a collisionless plasma to two modes which resonate with hot ion beams. The application of the theory is strictly limited to the linear growth phase. However, the application of the theory may be extended to obtain a description of the beam at postsaturation if the wave-beam resonance is sufficiently broad in velocity space. Under the considered limitations, it is shown that, as in the cold beam case, the fluctuating fields do not gain appreciable momentum and that the primary exchange of momentum is between the beam and main component.
Visualization of 3-D tensor fields
NASA Technical Reports Server (NTRS)
Hesselink, L.
1996-01-01
Second-order tensor fields have applications in many different areas of physics, such as general relativity and fluid mechanics. The wealth of multivariate information in tensor fields makes them more complex and abstract than scalar and vector fields. Visualization is a good technique for scientists to gain new insights from them. Visualizing a 3-D continuous tensor field is equivalent to simultaneously visualizing its three eigenvector fields. In the past, research has been conducted in the area of two-dimensional tensor fields. It was shown that degenerate points, defined as points where eigenvalues are equal to each other, are the basic singularities underlying the topology of tensor fields. Moreover, it was shown that eigenvectors never cross each other except at degenerate points. Since we live in a three-dimensional world, it is important for us to understand the underlying physics of this world. In this report, we describe a new method for locating degenerate points along with the conditions for classifying them in three-dimensional space. Finally, we discuss some topological features of three-dimensional tensor fields, and interpret topological patterns in terms of physical properties.
First and second order operator splitting methods for the phase field crystal equation
Lee, Hyun Geun; Shin, Jaemin; Lee, June-Yub
2015-10-15
In this paper, we present operator splitting methods for solving the phase field crystal equation which is a model for the microstructural evolution of two-phase systems on atomic length and diffusive time scales. A core idea of the methods is to decompose the original equation into linear and nonlinear subequations, in which the linear subequation has a closed-form solution in the Fourier space. We apply a nonlinear Newton-type iterative method to solve the nonlinear subequation at the implicit time level and thus a considerably large time step can be used. By combining these subequations, we achieve the first- and second-order accuracy in time. We present numerical experiments to show the accuracy and efficiency of the proposed methods.
NASA Astrophysics Data System (ADS)
Ito, Shin-Ichi; Nagao, Hiromichi; Yamanaka, Akinori; Tsukada, Yuhki; Koyama, Toshiyuki; Inoue, Junya
Phase field (PF) method, which phenomenologically describes dynamics of microstructure evolutions during solidification and phase transformation, has progressed in the fields of hydromechanics and materials engineering. How to determine, based on observation data, an initial state and model parameters involved in a PF model is one of important issues since previous estimation methods require too much computational cost. We propose data assimilation (DA), which enables us to estimate the parameters and states by integrating the PF model and observation data on the basis of the Bayesian statistics. The adjoint method implemented on DA not only finds an optimum solution by maximizing a posterior distribution but also evaluates the uncertainty in the estimations by utilizing the second order information of the posterior distribution. We carried out an estimation test using synthetic data generated by the two-dimensional Kobayashi's PF model. The proposed method is confirmed to reproduce the true initial state and model parameters we assume in advance, and simultaneously estimate their uncertainties due to quality and quantity of the data. This result indicates that the proposed method is capable of suggesting the experimental design to achieve the required accuracy.
NASA Astrophysics Data System (ADS)
Amore, Paolo; Boyd, John P.; Fernández, Francisco M.; Rösler, Boris
2016-05-01
We apply second order finite differences to calculate the lowest eigenvalues of the Helmholtz equation, for complicated non-tensor domains in the plane, using different grids which sample exactly the border of the domain. We show that the results obtained applying Richardson and Padé-Richardson extrapolations to a set of finite difference eigenvalues corresponding to different grids allow us to obtain extremely precise values. When possible we have assessed the precision of our extrapolations comparing them with the highly precise results obtained using the method of particular solutions. Our empirical findings suggest an asymptotic nature of the FD series. In all the cases studied, we are able to report numerical results which are more precise than those available in the literature.
Second order tensor finite element
NASA Technical Reports Server (NTRS)
Oden, J. Tinsley; Fly, J.; Berry, C.; Tworzydlo, W.; Vadaketh, S.; Bass, J.
1990-01-01
The results of a research and software development effort are presented for the finite element modeling of the static and dynamic behavior of anisotropic materials, with emphasis on single crystal alloys. Various versions of two dimensional and three dimensional hybrid finite elements were implemented and compared with displacement-based elements. Both static and dynamic cases are considered. The hybrid elements developed in the project were incorporated into the SPAR finite element code. In an extension of the first phase of the project, optimization of experimental tests for anisotropic materials was addressed. In particular, the problem of calculating material properties from tensile tests and of calculating stresses from strain measurements were considered. For both cases, numerical procedures and software for the optimization of strain gauge and material axes orientation were developed.
Relativistic Lagrangian displacement field and tensor perturbations
NASA Astrophysics Data System (ADS)
Rampf, Cornelius; Wiegand, Alexander
2014-12-01
We investigate the purely spatial Lagrangian coordinate transformation from the Lagrangian to the basic Eulerian frame. We demonstrate three techniques for extracting the relativistic displacement field from a given solution in the Lagrangian frame. These techniques are (a) from defining a local set of Eulerian coordinates embedded into the Lagrangian frame; (b) from performing a specific gauge transformation; and (c) from a fully nonperturbative approach based on the Arnowitt-Deser-Misner (ADM) split. The latter approach shows that this decomposition is not tied to a specific perturbative formulation for the solution of the Einstein equations. Rather, it can be defined at the level of the nonperturbative coordinate change from the Lagrangian to the Eulerian description. Studying such different techniques is useful because it allows us to compare and develop further the various approximation techniques available in the Lagrangian formulation. We find that one has to solve the gravitational wave equation in the relativistic analysis, otherwise the corresponding Newtonian limit will necessarily contain spurious nonpropagating tensor artifacts at second order in the Eulerian frame. We also derive the magnetic part of the Weyl tensor in the Lagrangian frame, and find that it is not only excited by gravitational waves but also by tensor perturbations which are induced through the nonlinear frame dragging. We apply our findings to calculate for the first time the relativistic displacement field, up to second order, for a Λ CDM Universe in the presence of a local primordial non-Gaussian component. Finally, we also comment on recent claims about whether mass conservation in the Lagrangian frame is violated.
Covariant second-order perturbations in generalized two-field inflation
Tzavara, Eleftheria; Tent, Bartjan van; Mizuno, Shuntaro E-mail: Shuntaro.Mizuno@apc.univ-paris7.fr
2014-07-01
We examine the covariant properties of generalized models of two-field inflation, with non-canonical kinetic terms and a possibly non-trivial field metric. We demonstrate that kinetic-term derivatives and covariant field derivatives do commute in a proper covariant framework, which was not realized before in the literature. We also define a set of generalized slow-roll parameters, using a unified notation. Within this framework, we study the most general class of models that allows for well-defined adiabatic and entropic sound speeds, which we identify as the models with parallel momentum and field velocity vectors. For these models we write the exact cubic action in terms of the adiabatic and isocurvature perturbations. We thus provide the tool to calculate the exact non-Gaussianity beyond slow-roll and at any scale for these generalized models. We illustrate our general results by considering their long-wavelength limit, as well as with the example of two-field DBI inflation.
On p -form theories with gauge invariant second order field equations
NASA Astrophysics Data System (ADS)
Deffayet, Cédric; Mukohyama, Shinji; Sivanesan, Vishagan
2016-04-01
We explore field theories of a single p -form with equations of motions of order strictly equal to 2 and gauge invariance. We give a general method for the classification of such theories which are extensions to the p -forms of the Galileon models for scalars. Our classification scheme allows us to compute an upper bound on the number of different such theories depending on p and on the space-time dimension. We are also able to build a nontrivial Galileon-like theory for a 3-form with gauge invariance and an action which is polynomial into the derivatives of the form. This theory has gauge invariant field equations but an action which is not, like a Chern-Simons theory. Hence the recently discovered no-go theorem stating that there are no nontrivial gauge invariant vector Galileons (which we are also able here to confirm with our method) does not extend to other odd-p cases.
Long-wavelength properties of phase-field-crystal models with second-order dynamics
NASA Astrophysics Data System (ADS)
Heinonen, V.; Achim, C. V.; Ala-Nissila, T.
2016-05-01
The phase-field-crystal (PFC) approach extends the notion of phase-field models by describing the topology of the microscopic structure of a crystalline material. One of the consequences is that local variation of the interatomic distance creates an elastic excitation. The dynamics of these excitations poses a challenge: pure diffusive dynamics cannot describe relaxation of elastic stresses that happen through phonon emission. To this end, several different models with fast dynamics have been proposed. In this article we use the amplitude expansion of the PFC model to compare the recently proposed hydrodynamic PFC amplitude model with two simpler models with fast dynamics. We compare these different models analytically and numerically. The results suggest that in order to have proper relaxation of elastic excitations, the full hydrodynamical description of the PFC amplitudes is required.
NASA Astrophysics Data System (ADS)
Dias, F.; Assumpcao, M.
2012-12-01
The knowledge of stress field is fundamental not only to understand driving forces and plate deformation as also it helps in the study of intraplate seismicity. In Brazil, we find reverse, strike-slip and normal mechanisms that indicates a variable stress field. The stress field has been mainly obtained using focal mechanism results and a few breakout data and in-situ measurements. However the stress field is still poorly known in Brazil. Recent earthquake focal mechanisms were determinate using P-wave modeling of seismogram stacks of several teleseismic stations ( > 30°) grouped according to distance and azimuth and first motion polarities. Every record was visually inspected and those with a good signal/noise ratio (SNR) were grouped in latitude-longitude windows of ten degrees and stacked. We usually consider groups with at least two stations, but, in sometimes a good record of single station with different azimuth was also used to constrain the focal depth. The P, pP, sP wavetrains of the stacked signals were modelled using the hudson96 program of Herrman seismology package (Herrman, 2002). We also determinate moment tensor of same events in the central region. The major difficulty is to determinate focal mechanism of low magnitudes events (< 4.0 mb) using distants seismograph stations. The central region shows a purely compressional pattern which are predicted by regional theoretical models (Richardson & Coblentz, 1996 and the TD0 model of Lithgow& Bertelloni, 2004). Meanwhile in the Amazonic region we find a SHmax from E-W to SE-NW probably caused by Caribbean and South American plates interaction (Meijer, 1995). In NE region, the compression rotates following the coast line which indicates an important component regional present in stress field spreading effects due to the continental/oceanic crustal (Assumpção, 1998) and cases of stress caused by sedimentary load in Amazon Fan in agreement local theoretical models (Watts et al., 2009). We determinate the
NASA Technical Reports Server (NTRS)
Liu, Ansheng; Chuang, S.-L.; Ning, C. Z.; Woo, Alex (Technical Monitor)
1999-01-01
Second-order nonlinear optical processes including second-harmonic generation, optical rectification, and difference-frequency generation associated with intersubband transitions in wurtzite GaN/AlGaN quantum well (QW) are investigated theoretically. Taking into account the strain-induced piezoelectric (PZ) effects, we solve the electronic structure of the QW from coupled effective-mass Schrodinger equation and Poisson equation including the exchange-correlation effect under the local-density approximation. We show that the large PZ field in the QW breaks the symmetry of the confinement potential profile and leads to large second-order susceptibilities. We also show that the interband optical pump-induced electron-hole plasma results in an enhancement in the maximum value of the nonlinear coefficients and a redshift of the peak position in the nonlinear optical spectrum. By use of the difference-frequency generation, THz radiation can be generated from a GaN/Al(0.75)Ga(0.25)N with a pump laser of 1.55 micron.
NASA Astrophysics Data System (ADS)
Abdul Hakeem, A. K.; Vishnu Ganesh, N.; Ganga, B.
2015-05-01
The magnetic field effect on a steady two dimensional laminar radiative flow of an incompressible viscous water based nanofluid over a stretching/shrinking sheet with second order slip boundary condition is investigated both analytically and numerically. The governing partial differential equations are reduced to nonlinear ordinary differential equations by means of Lie symmetry group transformations. The dimensionless governing equations for this investigation are solved analytically using hyper-geometric function and numerically by the fourth order Runge-Kutta method with the shooting technique. A unique exact solution exists for momentum equation in stretching sheet case and dual solutions are obtained for shrinking sheet case which has upper and lower branches. It is found that the lower branch solution vanishes in the presence of higher magnetic field. The velocity and temperature profiles, the local skin friction coefficient and the reduced Nusselt number are examined and discussed for different spherical nanoparticles such as Au, Ag, Cu, Al, Al2 O3 and TiO2. A comparative study between the previously published results and the present analytical and numerical results for a special case is found to be in good agreement.
NASA Astrophysics Data System (ADS)
Ouyang, Wei; Mao, Weijian; Li, Xuelei; Li, Wuqun
2014-08-01
Sound velocity inversion problem based on scattering theory is formulated in terms of a nonlinear integral equation associated with scattered field. Because of its nonlinearity, in practice, linearization algorisms (Born/single scattering approximation) are widely used to obtain an approximate inversion solution. However, the linearized strategy is not congruent with seismic wave propagation mechanics in strong perturbation (heterogeneous) medium. In order to partially dispense with the weak perturbation assumption of the Born approximation, we present a new approach from the following two steps: firstly, to handle the forward scattering by taking into account the second-order Born approximation, which is related to generalized Radon transform (GRT) about quadratic scattering potential; then to derive a nonlinear quadratic inversion formula by resorting to inverse GRT. In our formulation, there is a significant quadratic term regarding scattering potential, and it can provide an amplitude correction for inversion results beyond standard linear inversion. The numerical experiments demonstrate that the linear single scattering inversion is only good in amplitude for relative velocity perturbation () of background media up to 10 %, and its inversion errors are unacceptable for the perturbation beyond 10 %. In contrast, the quadratic inversion can give more accurate amplitude-preserved recovery for the perturbation up to 40 %. Our inversion scheme is able to manage double scattering effects by estimating a transmission factor from an integral over a small area, and therefore, only a small portion of computational time is added to the original linear migration/inversion process.
Curvature tensors unified field equations on SEXn
NASA Astrophysics Data System (ADS)
Chung, Kyung Tae; Lee, Il Young
1988-09-01
We study the curvature tensors and field equations in the n-dimensional SE manifold SEXn. We obtain several basic properties of the vectors S λ and U λ and then of the SE curvature tensor and its contractions, such as a generalized Ricci identity, a generalized Bianchi identity, and two variations of the Bianchi identity satisfied by the SE Einstein tensor. Finally, a system of field equations is discussed in SEXn and one of its particular solutions is constructed and displayed.
Liu, Gang; Jayathilake, Pahala Gedara; Khoo, Boo Cheong
2014-02-01
Two nonlinear models are proposed to investigate the focused acoustic waves that the nonlinear effects will be important inside the liquid around the scatterer. Firstly, the one dimensional solutions for the widely used Westervelt equation with different coordinates are obtained based on the perturbation method with the second order nonlinear terms. Then, by introducing the small parameter (Mach number), a dimensionless formulation and asymptotic perturbation expansion via the compressible potential flow theory is applied. This model permits the decoupling between the velocity potential and enthalpy to second order, with the first potential solutions satisfying the linear wave equation (Helmholtz equation), whereas the second order solutions are associated with the linear non-homogeneous equation. Based on the model, the local nonlinear effects of focused acoustic waves on certain volume are studied in which the findings may have important implications for bubble cavitation/initiation via focused ultrasound called HIFU (High Intensity Focused Ultrasound). The calculated results show that for the domain encompassing less than ten times the radius away from the center of the scatterer, the non-linear effect exerts a significant influence on the focused high intensity acoustic wave. Moreover, at the comparatively higher frequencies, for the model of spherical wave, a lower Mach number may result in stronger nonlinear effects. PMID:24070825
Fast stray field computation on tensor grids
Exl, L.; Auzinger, W.; Bance, S.; Gusenbauer, M.; Reichel, F.; Schrefl, T.
2012-01-01
A direct integration algorithm is described to compute the magnetostatic field and energy for given magnetization distributions on not necessarily uniform tensor grids. We use an analytically-based tensor approximation approach for function-related tensors, which reduces calculations to multilinear algebra operations. The algorithm scales with N4/3 for N computational cells used and with N2/3 (sublinear) when magnetization is given in canonical tensor format. In the final section we confirm our theoretical results concerning computing times and accuracy by means of numerical examples. PMID:24910469
Lazur, V. Yu.; Myhalyna, S. I.; Reity, O. K.
2010-06-15
The problem of interaction of two quasimolecular electrons located at an arbitrary distance from each other and near different atoms (nuclei) is solved. The interaction is considered as a second-order effect of quantum electrodynamics in the coordinate representation. It is shown that consistent account of the natural condition of the interaction symmetry with respect to both electrons leads to the additional contribution to the relativistic interaction of the two quasimolecular electrons compared with both the standard Breit operator and the generalized Breit operator obtained previously [O. N. Gadomskii, Usp. Fiz. Nauk 170, 1145 (2000) [Phys. Usp. 43, 1071 (2000)
NASA Astrophysics Data System (ADS)
Espin, Johnny; Krasnov, Kirill
2015-06-01
It is known, though not commonly, that one can describe fermions using a second order in derivatives Lagrangian instead of the first order Dirac one. In this description the propagator is scalar, and the complexity is shifted to the vertex, which contains a derivative operator. In this paper we rewrite the Lagrangian of the fermionic sector of the Standard Model in such second order form. The new Lagrangian is extremely compact, and is obtained from the usual first order Lagrangian by integrating out all primed (or dotted) 2-component spinors. It thus contains just half of the 2-component spinors that appear in the usual Lagrangian, which suggests a new perspective on unification. We sketch a natural in this framework SU (2) × SU (4) ⊂ SO (9) unified theory.
Zharkov, G. F.
2001-06-01
Based on self-consistent solution of nonlinear GL equations, the phase boundary is found, which divides the regions of first- and second-order phase transitions to normal state of a superconducting cylinder of radius R, placed in magnetic field and remaining in the state of fixed vorticity m. This boundary is a complicated function of the parameters (m,R,{kappa}) ({kappa} is the GL parameter), which does not coincide with the simple phase boundary {kappa}=1/{radical}2, dividing the regions of first- and second-order phase transitions in infinite (open) superconducting systems.
Invariant Crease Lines for Topological and Structural Analysis of Tensor Fields
Tricoche, Xavier; Kindlmann, Gordon; Westin, Carl-Fredrik
2009-01-01
We introduce a versatile framework for characterizing and extracting salient structures in three-dimensional symmetric second-order tensor fields. The key insight is that degenerate lines in tensor fields, as defined by the standard topological approach, are exactly crease (ridge and valley) lines of a particular tensor invariant called mode. This reformulation allows us to apply well-studied approaches from scientific visualization or computer vision to the extraction of topological lines in tensor fields. More generally, this main result suggests that other tensor invariants, such as anisotropy measures like fractional anisotropy (FA), can be used in the same framework in lieu of mode to identify important structural properties in tensor fields. Our implementation addresses the specific challenge posed by the non-linearity of the considered scalar measures and by the smoothness requirement of the crease manifold computation. We use a combination of smooth reconstruction kernels and adaptive refinement strategy that automatically adjust the resolution of the analysis to the spatial variation of the considered quantities. Together, these improvements allow for the robust application of existing ridge line extraction algorithms in the tensor context of our problem. Results are proposed for a diffusion tensor MRI dataset, and for a benchmark stress tensor field used in engineering research. PMID:18989019
NASA Technical Reports Server (NTRS)
Landahl, M.; Loefgren, P.
1973-01-01
A second-order theory for supersonic flow past slender bodies is presented. Through the introduction of characteristic coordinates as independent variables and the expansion procedure proposed by Lin and Oswatitsch, a uniformly valid solution is obtained for the whole flow field in the axisymmetric case and for far field in the general three-dimensional case. For distances far from the body the theory is an extension of Whitham's first-order solution and for the domain close to the body it is a modification of Van Dyke's second-order solution in the axisymmetric case. From the theory useful formulas relating flow deflections to the Whitham F-function are derived, which permits one to determine the sonic boom strength from wind tunnel measurements fairly close to the body.
Visualization of tensor fields using superquadric glyphs.
Ennis, Daniel B; Kindlman, Gordon; Rodriguez, Ignacio; Helm, Patrick A; McVeigh, Elliot R
2005-01-01
The spatially varying tensor fields that arise in magnetic resonance imaging are difficult to visualize due to the multivariate nature of the data. To improve the understanding of myocardial structure and function a family of objects called glyphs, derived from superquadric parametric functions, are used to create informative and intuitive visualizations of the tensor fields. The superquadric glyphs are used to visualize both diffusion and strain tensors obtained in canine myocardium. The eigensystem of each tensor defines the glyph shape and orientation. Superquadric functions provide a continuum of shapes across four distinct eigensystems (lambda(i), sorted eigenvalues), lambda(1) = lambda(2) = lambda(3) (spherical), lambda(1) < lambda(2) = lambda(3) (oblate), lambda(1) > lambda(2) = lambda(3) (prolate), and lambda(1) > lambda(2) > lambda(3) (cuboid). The superquadric glyphs are especially useful for identifying regions of anisotropic structure and function. Diffusion tensor renderings exhibit fiber angle trends and orthotropy (three distinct eigenvalues). Visualization of strain tensors with superquadric glyphs compactly exhibits radial thickening gradients, circumferential and longitudinal shortening, and torsion combined. The orthotropic nature of many biologic tissues and their DTMRI and strain data require visualization strategies that clearly exhibit the anisotropy of the data if it is to be interpreted properly. Superquadric glyphs improve the ability to distinguish fiber orientation and tissue orthotropy compared to ellipsoids. PMID:15690516
Bates, Jefferson E.; Shiozaki, Toru
2015-01-28
We develop an efficient algorithm for four-component complete active space self-consistent field (CASSCF) methods on the basis of the Dirac equation that takes into account spin–orbit and other relativistic effects self-consistently. Orbitals are optimized using a trust-region quasi-Newton method with Hessian updates so that energies are minimized with respect to rotations among electronic orbitals and maximized with respect to rotations between electronic and positronic orbitals. Utilizing density fitting and parallel computation, we demonstrate that Dirac–Coulomb CASSCF calculations can be routinely performed on systems with 100 atoms and a few heavy-elements. The convergence behavior and wall times for octachloridodirhenate(III) and a tungsten methylidene complex are presented. In addition, the excitation energies of octachloridodirhenate(III) are reported using a state-averaged variant.
NASA Technical Reports Server (NTRS)
Landahl, M.; Soerensen, H.; Hilding, L.
1973-01-01
An experimental investigation has been carried out in a wind tunnel to test some of the results of Landahl's second order theory. The slender models consisted of a parabolic spindle, tested at M = 3, and a wing body configuration, suggested by Ferri, and tested at M = 2.7. The theory indicates that shock position and strength at an arbitrary distance can be calculated by means of near field measurements. The results show that this method is an appropriate one for simple bodies and for bodies with complicated geometries as well.
The second-order gravitational red shift
NASA Technical Reports Server (NTRS)
Jaffe, J.
1973-01-01
The direct measurement of the nonlinear term of the gravitational field equations by using very stable clocks is discussed along with measuring the perhelion advance of a planet or satellite. These are considered measurements of the second-order gravitational red shift. The exact expression for the frequency shift of light in a gravitational field is derived. Other topics discussed include: The Doppler-cancelling technique; the second-order red shift in a spherically symmetric gravitational field; finite signal transit time; and the reality and interpretation of coordinates in the second-order red shift experiment.
Numazaki, Kazuya; Imai, Hiromitsu; Morinaga, Atsuo
2010-03-15
The second-order Zeeman effect of the sodium clock transition in a weak magnetic field of less than 50 {mu}T was measured as the scalar Aharonov-Bohm phase by two-photon stimulated Raman atom interferometry. The ac Stark effect of the Raman pulse was canceled out by adopting an appropriate intensity ratio of two photons in the Raman pulse. The Ramsey fringes for the pulse separation of 7 ms were obtained with a phase uncertainty of {pi}/200 rad. The nondispersive feature of the scalar Aharonov-Bohm phase was clearly demonstrated through 18 fringes with constant amplitude. The Breit-Rabi formula of the sodium clock transition was verified to be {Delta}{nu}=(0.222{+-}0.003)x10{sup 12}xB{sup 1.998{+-}0.004} in a magnetic field of less than 50 {mu}T.
A finite field method for calculating molecular polarizability tensors for arbitrary multipole rank.
Elking, Dennis M; Perera, Lalith; Duke, Robert; Darden, Thomas; Pedersen, Lee G
2011-11-30
A finite field method for calculating spherical tensor molecular polarizability tensors α(lm;l'm') = ∂Δ(lm)/∂ϕ(l'm')* by numerical derivatives of induced molecular multipole Δ(lm) with respect to gradients of electrostatic potential ϕ(l'm')* is described for arbitrary multipole ranks l and l'. Interconversion formulae for transforming multipole moments and polarizability tensors between spherical and traceless Cartesian tensor conventions are derived. As an example, molecular polarizability tensors up to the hexadecapole-hexadecapole level are calculated for water using the following ab initio methods: Hartree-Fock (HF), Becke three-parameter Lee-Yang-Parr exchange-correlation functional (B3LYP), Møller-Plesset perturbation theory up to second order (MP2), and Coupled Cluster theory with single and double excitations (CCSD). In addition, intermolecular electrostatic and polarization energies calculated by molecular multipoles and polarizability tensors are compared with ab initio reference values calculated by the Reduced Variation Space method for several randomly oriented small molecule dimers separated by a large distance. It is discussed how higher order molecular polarizability tensors can be used as a tool for testing and developing new polarization models for future force fields. PMID:21915883
The Topology of Three-Dimensional Symmetric Tensor Fields
NASA Technical Reports Server (NTRS)
Lavin, Yingmei; Levy, Yuval; Hesselink, Lambertus
1994-01-01
We study the topology of 3-D symmetric tensor fields. The goal is to represent their complex structure by a simple set of carefully chosen points and lines analogous to vector field topology. The basic constituents of tensor topology are the degenerate points, or points where eigenvalues are equal to each other. First, we introduce a new method for locating 3-D degenerate points. We then extract the topological skeletons of the eigenvector fields and use them for a compact, comprehensive description of the tensor field. Finally, we demonstrate the use of tensor field topology for the interpretation of the two-force Boussinesq problem.
Second-Order Footsteps Illusions
Anstis, Stuart
2015-01-01
In the “footsteps illusion”, light and dark squares travel at constant speed across black and white stripes. The squares appear to move faster and slower as their contrast against the stripes varies. We now demonstrate some second-order footsteps illusions, in which all edges are defined by colors or textures—even though luminance-based neural motion detectors are blind to such edges. PMID:27551366
Second-Order Footsteps Illusions.
Kitaoka, Akiyoshi; Anstis, Stuart
2015-12-01
In the "footsteps illusion", light and dark squares travel at constant speed across black and white stripes. The squares appear to move faster and slower as their contrast against the stripes varies. We now demonstrate some second-order footsteps illusions, in which all edges are defined by colors or textures-even though luminance-based neural motion detectors are blind to such edges. PMID:27551366
Drift kinetic equation exact through second order in gyroradius expansion
Simakov, Andrei N.; Catto, Peter J.
2005-01-01
The drift kinetic equation of Hazeltine [R. D. Hazeltine, Plasma Phys. 15, 77 (1973)] for a magnetized plasma of arbitrary collisionality is widely believed to be exact through the second order in the gyroradius expansion. It is demonstrated that this equation is only exact through the first order. The reason is that when evaluating the second-order gyrophase dependent distribution function, Hazeltine neglected contributions from the first-order gyrophase dependent distribution function, and then used this incomplete expression to derive the equation for the gyrophase independent distribution function. Consequently, the second-order distribution function and the stress tensor derived by this approach are incomplete. By relaxing slightly Hazeltine's orderings one is able to obtain a drift kinetic equation accurate through the second order in the gyroradius expansion. In addition, the gyroviscous stress tensor for plasmas of arbitrary collisionality is obtained.
Feature Surfaces in Symmetric Tensor Fields Based on Eigenvalue Manifold.
Palacios, Jonathan; Yeh, Harry; Wang, Wenping; Zhang, Yue; Laramee, Robert S; Sharma, Ritesh; Schultz, Thomas; Zhang, Eugene
2016-03-01
Three-dimensional symmetric tensor fields have a wide range of applications in solid and fluid mechanics. Recent advances in the (topological) analysis of 3D symmetric tensor fields focus on degenerate tensors which form curves. In this paper, we introduce a number of feature surfaces, such as neutral surfaces and traceless surfaces, into tensor field analysis, based on the notion of eigenvalue manifold. Neutral surfaces are the boundary between linear tensors and planar tensors, and the traceless surfaces are the boundary between tensors of positive traces and those of negative traces. Degenerate curves, neutral surfaces, and traceless surfaces together form a partition of the eigenvalue manifold, which provides a more complete tensor field analysis than degenerate curves alone. We also extract and visualize the isosurfaces of tensor modes, tensor isotropy, and tensor magnitude, which we have found useful for domain applications in fluid and solid mechanics. Extracting neutral and traceless surfaces using the Marching Tetrahedra method can cause the loss of geometric and topological details, which can lead to false physical interpretation. To robustly extract neutral surfaces and traceless surfaces, we develop a polynomial description of them which enables us to borrow techniques from algebraic surface extraction, a topic well-researched by the computer-aided design (CAD) community as well as the algebraic geometry community. In addition, we adapt the surface extraction technique, called A-patches, to improve the speed of finding degenerate curves. Finally, we apply our analysis to data from solid and fluid mechanics as well as scalar field analysis. PMID:26441450
Properties of second-order geometrical aberrations
NASA Astrophysics Data System (ADS)
Grammatin, A. P.
1994-08-01
This paper analyzes the properties of second-order aberrations that arise in centered optical systems that contain an aspherical surface whose sagittal equation contains a term proportional to the cube of the distance from a surface point to the optical axis. It is shown that the second-order spherical aberration decreases from the center of the field to its edge. No astigmatism appears in wide, oblique beams in the central part of the field. Coma increases linearly from zero at the center of the field to a value equal to the spherical aberration, and then remains constant over the field. A proof is given of the possibility of correcting the image curvature by using an aspherical surface of the type described above.
K-inflationary power spectra at second order
Martin, Jérôme; Vennin, Vincent; Ringeval, Christophe E-mail: christophe.ringeval@uclouvain.be
2013-06-01
Within the class of inflationary models, k-inflation represents the most general single field framework that can be associated with an effective quadratic action for the curvature perturbations and a varying speed of sound. The incoming flow of high-precision cosmological data, such as those from the Planck satellite and small scale Cosmic Microwave Background (CMB) experiments, calls for greater accuracy in the inflationary predictions. In this work, we calculate for the first time the next-to-next-to-leading order scalar and tensor primordial power spectra in k-inflation needed in order to obtain robust constraints on the inflationary theory. The method used is the uniform approximation together with a second order expansion in the Hubble and sound flow functions. Our result is checked in various limits in which it reduces to already known situations.
NASA Astrophysics Data System (ADS)
Nakamura, K.
2007-01-01
Following the general framework of the gauge invariant perturbation theory developed in the papers [K. Nakamura, Prog. Theor. Phys. 110 (2003), 723; ibid. 113 (2005), 481], we formulate second-order gauge invariant cosmological perturbation theory in a four-dimensional homogeneous isotropic universe. We consider perturbations both in the universe dominated by a single perfect fluid and in that dominated by a single scalar field. We derive all the components of the Einstein equations in the case that the first-order vector and tensor modes are negligible. All equations are derived in terms of gauge invariant variables without any gauge fixing. These equations imply that second-order vector and tensor modes may be generated due to the mode-mode coupling of the linear-order scalar perturbations. We also briefly discuss the main progress of this work through comparison with previous works.
Higher rank antisymmetric tensor fields in Klebanov-Strassler geometry
NASA Astrophysics Data System (ADS)
Das, Ashmita; SenGupta, Soumitra
2016-05-01
In string theory, higher rank antisymmetric tensor fields appear as massless excitations of closed strings. To date, there is no experimental support in favor of their existence. In a stringy framework, starting from a warped throatlike Klebanov-Strassler geometry, we show that all the massless higher rank antisymmetric tensor fields are heavily suppressed due to the background fluxes leading to their invisibility in our Universe.
Tensor classification of structure in smoothed particle hydrodynamics density fields
NASA Astrophysics Data System (ADS)
Forgan, Duncan; Bonnell, Ian; Lucas, William; Rice, Ken
2016-04-01
As hydrodynamic simulations increase in scale and resolution, identifying structures with non-trivial geometries or regions of general interest becomes increasingly challenging. There is a growing need for algorithms that identify a variety of different features in a simulation without requiring a `by eye' search. We present tensor classification as such a technique for smoothed particle hydrodynamics (SPH). These methods have already been used to great effect in N-Body cosmological simulations, which require smoothing defined as an input free parameter. We show that tensor classification successfully identifies a wide range of structures in SPH density fields using its native smoothing, removing a free parameter from the analysis and preventing the need for tessellation of the density field, as required by some classification algorithms. As examples, we show that tensor classification using the tidal tensor and the velocity shear tensor successfully identifies filaments, shells and sheet structures in giant molecular cloud simulations, as well as spiral arms in discs. The relationship between structures identified using different tensors illustrates how different forces compete and co-operate to produce the observed density field. We therefore advocate the use of multiple tensors to classify structure in SPH simulations, to shed light on the interplay of multiple physical processes.
NASA Astrophysics Data System (ADS)
Beig, Robert; Krammer, Werner
2004-02-01
For a conformally flat 3-space, we derive a family of linear second-order partial differential operators which sends vectors into trace-free, symmetric 2-tensors. These maps, which are parametrized by conformal Killing vectors on the 3-space, are such that the divergence of the resulting tensor field depends only on the divergence of the original vector field. In particular, these maps send source-free electric fields into TT tensors. Moreover, if the original vector field is the Coulomb field on {\\bb R}^3\\backslash \\lbrace0\\rbrace , the resulting tensor fields on {\\bb R}^3\\backslash \\lbrace0\\rbrace are nothing but the family of TT tensors originally written by Bowen and York.
Calculating Second-Order Effects in MOSFET's
NASA Technical Reports Server (NTRS)
Benumof, Reuben; Zoutendyk, John A.; Coss, James R.
1990-01-01
Collection of mathematical models includes second-order effects in n-channel, enhancement-mode, metal-oxide-semiconductor field-effect transistors (MOSFET's). When dimensions of circuit elements relatively large, effects neglected safely. However, as very-large-scale integration of microelectronic circuits leads to MOSFET's shorter or narrower than 2 micrometer, effects become significant in design and operation. Such computer programs as widely-used "Simulation Program With Integrated Circuit Emphasis, Version 2" (SPICE 2) include many of these effects. In second-order models of n-channel, enhancement-mode MOSFET, first-order gate-depletion region diminished by triangular-cross-section deletions on end and augmented by circular-wedge-cross-section bulges on sides.
Kubo Formulas for Second-Order Hydrodynamic Coefficients
Moore, Guy D.; Sohrabi, Kiyoumars A.
2011-03-25
At second order in gradients, conformal relativistic hydrodynamics depends on the viscosity {eta} and on five additional ''second-order'' hydrodynamical coefficients {tau}{sub {Pi}}, {kappa}, {lambda}{sub 1}, {lambda}{sub 2}, and {lambda}{sub 3}. We derive Kubo relations for these coefficients, relating them to equilibrium, fully retarded three-point correlation functions of the stress tensor. We show that the coefficient {lambda}{sub 3} can be evaluated directly by Euclidean means and does not in general vanish.
Antisymmetric tensor field and spontaneous magnetization in holographic duality
NASA Astrophysics Data System (ADS)
Cai, Rong-Gen; Yang, Run-Qiu
2015-08-01
A real antisymmetric tensor field was introduced to realize a holographic magnetic ordered phase in our previous papers. However, a more careful analysis shows there is a vector ghost in the model. In this paper we present a modified Lagrangian density for the antisymmetric tensor, which is ghost free and causality is well defined, and keeps all the significant results in the original model qualitatively. We show this modified Lagrangian density could come from the dimensional compactification of p -form field in string/M theory. For static curved space-time, we also prove that this modified model is ghost free and does not violate causality. This new model offers a solid foundation for the application of antisymmetric tensor field in holographic duality, especially for the spontaneous magnetization.
Bulk antisymmetric tensor fields in a Randall-Sundrum model
Mukhopadhyaya, Biswarup; Sen, Somasri; SenGupta, Soumitra
2007-12-15
We consider bulk antisymmetric tensor fields of various ranks in a Randall-Sundrum scenario. We show that, rank 2 onwards, the zero-modes of the projections of these fields on the (3+1)-dimensional visible brane become increasingly weaker as the rank of the tensor increases. All such tensor fields of rank 4 or more are absent from the dynamics in four dimensions. This leaves only the zero-mode graviton to have coupling {approx}1/M{sub P} with matter, thus explaining why the large-scale behavior of the universe is governed by gravity only. We have also computed the masses of the heavier modes up to rank 3, and have shown that they are relatively less likely to have detectable accelerator signals.
NASA Astrophysics Data System (ADS)
Dong, Hui; Qiu, Longqing; Shi, Wen; Chang, Baolin; Qiu, Yang; Xu, Lu; Liu, Chao; Zhang, Yi; Krause, Hans-Joachim; Offenhäusser, Andreas; Xie, Xiaoming
2013-03-01
An ultra-low field (ULF) magnetic resonance imaging (MRI) system was set up in an urban laboratory without magnetic shielding. The measured environmental gradient fields of 1 ˜ 5 μT/m caused image distortion. We designed a gradient detection and compensation system to effectively balance the gradient tensor components. The free induction decay signal duration of tap water was thus extended from 0.3 s to 2.5 s, providing the possibility for high-resolution imaging. Two-dimensional MRI images were then obtained at 130 μT with a helium-cooled second-order superconducting quantum interference device gradiometer. This result allows us to develop an inexpensive ULF MRI system for biological studies.
Relativistic second-order dissipative hydrodynamics at finite chemical potential
NASA Astrophysics Data System (ADS)
Jaiswal, Amaresh; Friman, Bengt; Redlich, Krzysztof
2015-12-01
Starting from the Boltzmann equation in the relaxation time approximation and employing a Chapman-Enskog like expansion for the distribution function close to equilibrium, we derive second-order evolution equations for the shear stress tensor and the dissipative charge current for a system of massless quarks and gluons. The transport coefficients are obtained exactly using quantum statistics for the phase space distribution functions at non-zero chemical potential. We show that, within the relaxation time approximation, the second-order evolution equations for the shear stress tensor and the dissipative charge current can be decoupled. We find that, for large values of the ratio of chemical potential to temperature, the charge conductivity is small compared to the coefficient of shear viscosity. Moreover, we show that in the relaxation-time approximation, the limiting behaviour of the ratio of heat conductivity to shear viscosity is qualitatively similar to that obtained for a strongly coupled conformal plasma.
Homological equations for tensor fields and periodic averaging
NASA Astrophysics Data System (ADS)
Avendaño Camacho, M.; Vorobiev, Y. M.
2011-09-01
Homological equations of tensor type associated to periodic flows on a manifold are studied. The Cushman intrinsic formula [4] is generalized to the case of multivector fields and differential forms. Some applications to normal forms and the averaging method for perturbed Hamiltonian systems on slow-fast phase spaces are given.
Source of second order chromaticity in RHIC
Luo, Y.; Gu, X.; Fischer, W.; Trbojevic, D.
2011-01-01
In this note we will answer the following questions: (1) what is the source of second order chromaticities in RHIC? (2) what is the dependence of second order chromaticity on the on-momentum {beta}-beat? (3) what is the dependence of second order chromaticity on {beta}* at IP6 and IP8? To answer these questions, we use the perturbation theory to numerically calculate the contributions of each quadrupole and sextupole to the first, second, and third order chromaticities.
Binocular Combination of Second-Order Stimuli
Zhou, Jiawei; Liu, Rong; Zhou, Yifeng; Hess, Robert F.
2014-01-01
Phase information is a fundamental aspect of visual stimuli. However, the nature of the binocular combination of stimuli defined by modulations in contrast, so-called second-order stimuli, is presently not clear. To address this issue, we measured binocular combination for first- (luminance modulated) and second-order (contrast modulated) stimuli using a binocular phase combination paradigm in seven normal adults. We found that the binocular perceived phase of second-order gratings depends on the interocular signal ratio as has been previously shown for their first order counterparts; the interocular signal ratios when the two eyes were balanced was close to 1 in both first- and second-order phase combinations. However, second-order combination is more linear than previously found for first-order combination. Furthermore, binocular combination of second-order stimuli was similar regardless of whether the carriers in the two eyes were correlated, anti-correlated, or uncorrelated. This suggests that, in normal adults, the binocular phase combination of second-order stimuli occurs after the monocular extracting of the second-order modulations. The sensory balance associated with this second-order combination can be obtained from binocular phase combination measurements. PMID:24404180
Binocular combination of second-order stimuli.
Zhou, Jiawei; Liu, Rong; Zhou, Yifeng; Hess, Robert F
2014-01-01
Phase information is a fundamental aspect of visual stimuli. However, the nature of the binocular combination of stimuli defined by modulations in contrast, so-called second-order stimuli, is presently not clear. To address this issue, we measured binocular combination for first- (luminance modulated) and second-order (contrast modulated) stimuli using a binocular phase combination paradigm in seven normal adults. We found that the binocular perceived phase of second-order gratings depends on the interocular signal ratio as has been previously shown for their first order counterparts; the interocular signal ratios when the two eyes were balanced was close to 1 in both first- and second-order phase combinations. However, second-order combination is more linear than previously found for first-order combination. Furthermore, binocular combination of second-order stimuli was similar regardless of whether the carriers in the two eyes were correlated, anti-correlated, or uncorrelated. This suggests that, in normal adults, the binocular phase combination of second-order stimuli occurs after the monocular extracting of the second-order modulations. The sensory balance associated with this second-order combination can be obtained from binocular phase combination measurements. PMID:24404180
NASA Astrophysics Data System (ADS)
Milton, Kimball A.; Fulling, Stephen A.; Parashar, Prachi; Kalauni, Pushpa; Murphy, Taylor
2016-04-01
Motivated by a desire to understand quantum fluctuation energy densities and stress within a spatially varying dielectric medium, we examine the vacuum expectation value for the stress tensor of a scalar field with arbitrary conformal parameter, in the background of a given potential that depends on only one spatial coordinate. We regulate the expressions by incorporating a temporal-spatial cutoff in the (imaginary) time and transverse-spatial directions. The divergences are captured by the zeroth- and second-order WKB approximations. Then the stress tensor is "renormalized" by omitting the terms that depend on the cutoff. The ambiguities that inevitably arise in this procedure are both duly noted and restricted by imposing certain physical conditions; one result is that the renormalized stress tensor exhibits the expected trace anomaly. The renormalized stress tensor exhibits no pressure anomaly, in that the principle of virtual work is satisfied for motions in a transverse direction. We then consider a potential that defines a wall, a one-dimensional potential that vanishes for z <0 and rises like zα, α >0 , for z >0 . Previously, the stress tensor had been computed outside of the wall, whereas now we compute all components of the stress tensor in the interior of the wall. The full finite stress tensor is computed numerically for the two cases where explicit solutions to the differential equation are available, α =1 and 2. The energy density exhibits an inverse linear divergence as the boundary is approached from the inside for a linear potential, and a logarithmic divergence for a quadratic potential. Finally, the interaction between two such walls is computed, and it is shown that the attractive Casimir pressure between the two walls also satisfies the principle of virtual work (i.e., the pressure equals the negative derivative of the energy with respect to the distance between the walls).
Symmetries of second-order PDEs and conformal Killing vectors
NASA Astrophysics Data System (ADS)
Tsamparlis, Michael; Paliathanasis, Andronikos
2015-06-01
We study the Lie point symmetries of a general class of partial differential equations (PDE) of second order. An equation from this class naturally defines a second-order symmetric tensor (metric). In the case the PDE is linear on the first derivatives we show that the Lie point symmetries are given by the conformal algebra of the metric modulo a constraint involving the linear part of the PDE. Important elements in this class are the Klein-Gordon equation and the Laplace equation. We apply the general results and determine the Lie point symmetries of these equations in various general classes of Riemannian spaces. Finally we study the type II hidden symmetries of the wave equation in a Riemannian space with a Lorenzian metric.
Second-order closure models for supersonic turbulent flows
NASA Technical Reports Server (NTRS)
Speziale, Charles G.; Sarkar, Sutanu
1991-01-01
Recent work on the development of a second-order closure model for high-speed compressible flows is reviewed. This turbulent closure is based on the solution of modeled transport equations for the Favre-averaged Reynolds stress tensor and the solenoidal part of the turbulent dissipation rate. A new model for the compressible dissipation is used along with traditional gradient transport models for the Reynolds heat flux and mass flux terms. Consistent with simple asymptotic analyses, the deviatoric part of the remaining higher-order correlations in the Reynolds stress transport equations are modeled by a variable density extension of the newest incompressible models. The resulting second-order closure model is tested in a variety of compressible turbulent flows which include the decay of isotropic turbulence, homogeneous shear flow, the supersonic mixing layer, and the supersonic flat-plate turbulent boundary layer. Comparisons between the model predictions and the results of physical and numerical experiments are quite encouraging.
Chen, Zhenhua; Chen, Xun; Ying, Fuming; Gu, Junjing; Zhang, Huaiyu; Wu, Wei
2014-10-01
Using the formulas and techniques developed in Papers I and II of this series, the recently developed second-order perturbation theory based on a valence bond self-consistent field reference function (VBPT2) has been extended by using the internally contracted correction wave function. This ansatz strongly reduces the size of the interaction space compared to the uncontracted wave function and thus improves the capability of the VBPT2 method dramatically. Test calculations show that internally contracted VBPT2 using only a small number of reference valence bond functions, can give results as accuracy as the VBPT2 method and other more sophisticated methods such as full configuration interaction and multireference configuration interaction. PMID:25296795
Cosmology in generalized Horndeski theories with second-order equations of motion
NASA Astrophysics Data System (ADS)
Kase, Ryotaro; Tsujikawa, Shinji
2014-08-01
We study the cosmology of an extended version of Horndeski theories with second-order equations of motion on the flat Friedmann-Lemaître-Robertson-Walker background. In addition to a dark energy field χ associated with the gravitational sector, we take into account multiple scalar fields ϕI (I =1,2,…,N-1) characterized by the Lagrangians P(I)(XI) with XI=∂μϕI∂μϕI. These additional scalar fields can model the perfect fluids of radiation and nonrelativistic matter. We derive propagation speeds of scalar and tensor perturbations as well as conditions for the absence of ghosts. The theories beyond Horndeski induce nontrivial modifications to all the propagation speeds of N scalar fields, but the modifications to those for the matter fields ϕI are generally suppressed relative to that for the dark energy field χ. We apply our results to the covariantized Galileon with an Einstein-Hilbert term in which partial derivatives of the Minkowski Galileon are replaced by covariant derivatives. Unlike the covariant Galileon with second-order equations of motion in general space-time, the scalar propagation speed square cs12 associated with the field χ becomes negative during the matter era for late-time tracking solutions, so the two Galileon theories can be clearly distinguished at the level of linear cosmological perturbations.
Static second-order polarizabilities of aminobenzophenones and nitrobenzophenones
NASA Technical Reports Server (NTRS)
Moore, Craig E.; Cardelino, Beatriz H.
1991-01-01
Static-field theoretical studies on molecular second-order polarizabilities (beta) of benzophenone derivatives were performed. Calculations were based on the use of shaped electric fields and semiempirical Hamiltonians. Either an electron-donating (amine) or an electron-withdrawing (nitro) substituent was incorporated into a phenyl ring of benzophenone; the phenyl rings of benzophenone were oriented either coplanar or perpendicular to the carbonyl. The change in charge transfer with respect to the electrophilic character of the carbonyl group was monitored to determine its effect on the molecular second-order polarizability. Calculations were performed for all constitutional isomers of the two benzophenone derivatives.
Second-order reconstruction of the inflationary potential
NASA Technical Reports Server (NTRS)
Liddle, Andrew R.; Turner, Michael S.
1994-01-01
To first order in the deviation from scale invariance the inflationary potential and its first two derivatives can be expressed in terms of the spectral indices of the scalar and tensor perturbations, n and n(sub T), and their contributions to the variance of the quadrupole CBR temperature anisotropy, S and T. In addition, there is a 'consistency relation' between these quantities: n(sub T) = (-1/ 7)(T/S). We derive the second-order expressions for the inflationary potential and its first two derivatives and the first-order expression for its third derivative, in terms, of n, n(sub T), S, T, and dn/d ln gamma. We also obtain the second-order consistency relation, n(sub T) = (-1/7)(T/S)(1 + 0.11(T/S) + 0.15(n-1)). As an example we consider the exponential potential, the only known case where exact analytic solutions for the perturbation spectra exist. We reconstruct the potential via Taylor expansion (with coefficients calculated at both first and second order), and introduce the Pade approximate as a greatly improved alternative.
Second-order (2 +1 ) -dimensional anisotropic hydrodynamics
NASA Astrophysics Data System (ADS)
Bazow, Dennis; Heinz, Ulrich; Strickland, Michael
2014-11-01
We present a complete formulation of second-order (2 +1 ) -dimensional anisotropic hydrodynamics. The resulting framework generalizes leading-order anisotropic hydrodynamics by allowing for deviations of the one-particle distribution function from the spheroidal form assumed at leading order. We derive complete second-order equations of motion for the additional terms in the macroscopic currents generated by these deviations from their kinetic definition using a Grad-Israel-Stewart 14-moment ansatz. The result is a set of coupled partial differential equations for the momentum-space anisotropy parameter, effective temperature, the transverse components of the fluid four-velocity, and the viscous tensor components generated by deviations of the distribution from spheroidal form. We then perform a quantitative test of our approach by applying it to the case of one-dimensional boost-invariant expansion in the relaxation time approximation (RTA) in which case it is possible to numerically solve the Boltzmann equation exactly. We demonstrate that the second-order anisotropic hydrodynamics approach provides an excellent approximation to the exact (0+1)-dimensional RTA solution for both small and large values of the shear viscosity.
Nine Practices of Second Order Schools
ERIC Educational Resources Information Center
Brown, Bill; Tucker, Patrick; Williams, Thomas L.
2012-01-01
Many schools are in some stage of implementing differentiated instruction, with some already in what Carol Tomlinson describes in "The Differentiated School" as "second order change," where the entire school practices differentiation. In high-performing schools, differentiation has proved to be an effective instructional strategy; in classroom…
Second-Order Conditioning in "Drosophila"
ERIC Educational Resources Information Center
Tabone, Christopher J.; de Belle, J. Steven
2011-01-01
Associative conditioning in "Drosophila melanogaster" has been well documented for several decades. However, most studies report only simple associations of conditioned stimuli (CS, e.g., odor) with unconditioned stimuli (US, e.g., electric shock) to measure learning or establish memory. Here we describe a straightforward second-order conditioning…
Correction of second order chromaticity at Tevatron
Valishev, A.; Annala, G.; Lebedev, V.; Moore, R.S.; /Fermilab
2007-06-01
Correction of the second order betatron tune chromaticity is essential for operation at the working point near half integer resonance which is proposed as one of the ways to improve performance of the Tevatron. In this report the new chromaticity correction scheme with split sextupole families is described. Details of implementation and commissioning at the present working point are discussed.
Second order perturbations during inflation beyond slow-roll
Huston, Ian; Malik, Karim A. E-mail: k.malik@qmul.ac.uk
2011-10-01
We numerically calculate the evolution of second order cosmological perturbations for an inflationary scalar field without resorting to the slow-roll approximation or assuming large scales. In contrast to previous approaches we therefore use the full non-slow-roll source term for the second order Klein-Gordon equation which is valid on all scales. The numerical results are consistent with the ones obtained previously where slow-roll is a good approximation. We investigate the effect of localised features in the scalar field potential which break slow-roll for some portion of the evolution. The numerical package solving the second order Klein-Gordon equation has been released under an open source license and is available for download.
NASA Astrophysics Data System (ADS)
Du, J.; Chen, C.; Lesur, V.; Wang, L.
2015-07-01
General expressions of magnetic vector (MV) and magnetic gradient tensor (MGT) in terms of the first- and second-order derivatives of spherical harmonics at different degrees/orders are relatively complicated and singular at the poles. In this paper, we derived alternative non-singular expressions for the MV, the MGT and also the third-order partial derivatives of the magnetic potential field in the local north-oriented reference frame. Using our newly derived formulae, the magnetic potential, vector and gradient tensor fields and also the third-order partial derivatives of the magnetic potential field at an altitude of 300 km are calculated based on a global lithospheric magnetic field model GRIMM_L120 (GFZ Reference Internal Magnetic Model, version 0.0) with spherical harmonic degrees 16-90. The corresponding results at the poles are discussed and the validity of the derived formulas is verified using the Laplace equation of the magnetic potential field.
NASA Astrophysics Data System (ADS)
Du, J.; Chen, C.; Lesur, V.; Wang, L.
2014-12-01
General expressions of magnetic vector (MV) and magnetic gradient tensor (MGT) in terms of the first- and second-order derivatives of spherical harmonics at different degrees and orders, are relatively complicated and singular at the poles. In this paper, we derived alternative non-singular expressions for the MV, the MGT and also the higher-order partial derivatives of the magnetic field in local north-oriented reference frame. Using our newly derived formulae, the magnetic potential, vector and gradient tensor fields at an altitude of 300 km are calculated based on a global lithospheric magnetic field model GRIMM_L120 (version 0.0) and the main magnetic field model of IGRF11. The corresponding results at the poles are discussed and the validity of the derived formulas is verified using the Laplace equation of the potential field.
NASA Astrophysics Data System (ADS)
Kim, Inkoo; Lee, Yoon Sup
2014-10-01
We report the formulation and implementation of KRCASPT2, a two-component multi-configurational second-order perturbation theory based on Kramers restricted complete active space self-consistent field (KRCASSCF) reference function, in the framework of the spin-orbit relativistic effective core potential. The zeroth-order Hamiltonian is defined as the sum of nondiagonal one-electron operators with generalized two-component Fock matrix elements as scalar factors. The Kramers symmetry within the zeroth-order Hamiltonian is maintained via the use of a state-averaged density, allowing a consistent treatment of degenerate states. The explicit expressions are derived for the matrix elements of the zeroth-order Hamiltonian as well as for the perturbation vector. The use of a fully variational reference function and nondiagonal operators in relativistic multi-configurational perturbation theory is reported for the first time. A series of initial calculations are performed on the ionization potential and excitation energies of the atoms of the 6p-block; the results display a significant improvement over those from KRCASSCF, showing a closer agreement with experimental results. Accurate atomic properties of the superheavy elements of the 7p-block are also presented, and the electronic structures of the low-lying excited states are compared with those of their lighter homologues.
Second-Order Invariants and Holography
NASA Astrophysics Data System (ADS)
Luongo, Orlando; Bonanno, Luca; Iannone, Gerardo
2012-12-01
Motivated by recent works on the role of the holographic principle in cosmology, we relate a class of second-order Ricci invariants to the IR cutoff characterizing the holographic dark energy density. The choice of second-order invariants provides an invariant way to account the problem of causality for the correct cosmological cutoff, since the presence of event horizons is not an a priori assumption. We find that these models work fairly well, by fitting the observational data, through a combined cosmological test with the use of SNeIa, BAO and CMB. This class of models is also able to overcome the fine-tuning and coincidence problems. Finally, to make a comparison with other recent models, we adopt the statistical tests AIC and BIC.
The Poisson equation at second order in relativistic cosmology
Hidalgo, J.C.; Christopherson, Adam J.; Malik, Karim A. E-mail: Adam.Christopherson@nottingham.ac.uk
2013-08-01
We calculate the relativistic constraint equation which relates the curvature perturbation to the matter density contrast at second order in cosmological perturbation theory. This relativistic ''second order Poisson equation'' is presented in a gauge where the hydrodynamical inhomogeneities coincide with their Newtonian counterparts exactly for a perfect fluid with constant equation of state. We use this constraint to introduce primordial non-Gaussianity in the density contrast in the framework of General Relativity. We then derive expressions that can be used as the initial conditions of N-body codes for structure formation which probe the observable signature of primordial non-Gaussianity in the statistics of the evolved matter density field.
Remarks on the second-order Seiberg-Witten maps
Trampetic, Josip; Wohlgenannt, Michael
2007-12-15
In this brief report, we discuss the Seiberg-Witten maps up to the second order in the noncommutative parameter {theta}. They add to the recently published solutions in [A. Alboteanu, T. Ohl, and R. Rueckl, Phys. Rev. D 76, 105018 (2007).]. Expressions for the vector, fermion, and Higgs fields are given explicitly.
Comments on the present state of second-order closure models for incompressible flows
NASA Technical Reports Server (NTRS)
Speziale, Charles G.
1992-01-01
Second-order closure models account for history and nonlocal effects of the mean velocity gradients on the Reynolds stress tensor. Turbulent flows involving body forces or curvature, Reynolds stress relaxational effects, and counter-gradient transport are usually better described. The topics are presented in viewgraph form and include: (1) the Reynolds stress transport equation; (2) issues in second-order closure modeling; and (3) near wall models.
Second order Kerr-Newman time delay
NASA Astrophysics Data System (ADS)
He, G.; Lin, W.
2016-01-01
The explicit form for the post-Newtonian gravitational time delay of light signals propagating on the equatorial plane of a Kerr-Newman black hole is derived. Based on the null geodesic in Kerr-Newman spacetime, we adopt the iterative method to calculate the time delay. Our result reduces to the previous formulation for the Kerr black hole if we drop the contribution from the electrical charge. Our time-delay formula for the Reissner-Nordström geometry is different from the previous publication [Phys. Rev. D 69, 023002 (2004)], in which the largest second order contribution to the time delay is missing.
Analysis of second-order gratings
Hardy, A.; Welch, D.F.; Streifer, W. )
1989-10-01
The authors report the results of a second-order grating analysis. The gratings are used as distributed Bragg reflectors in surface-emitting lasers, which are currently being fabricated in several laboratories. The gratings provide reflection, output coupling, and power transmission to other gain segments for purposes of injection locking. The analysis determines these quantities for arbitrary-shaped grating teeth and includes the presence of a substrate reflector to reduce the radiated power in that direction. The reflector is shown to be effective, but only if it can be precisely positioned. Examples illustrating variations in dimensions, tooth shapes and heights, waveguide loss, and detuning are included.
Second-order coherence of supercontinuum light.
Genty, Goëry; Surakka, Minna; Turunen, Jari; Friberg, Ari T
2010-09-15
We analyze the coherence properties of supercontinuum generated in photonic crystal fibers by applying the second-order coherence theory of nonstationary light. Using an ensemble of simulated realizations, we construct two-frequency cross-spectral density and two-time mutual coherence functions. This allows us to introduce measures of temporal and spectral coherence. We show that, in the long-pulse regime, supercontinuum light can be decomposed into a sum of coherent and quasi-stationary contributions. Our approach and findings are also applicable in the short-pulse regime. PMID:20847777
Relativistic quantum transport coefficients for second-order viscous hydrodynamics
NASA Astrophysics Data System (ADS)
Florkowski, Wojciech; Jaiswal, Amaresh; Maksymiuk, Ewa; Ryblewski, Radoslaw; Strickland, Michael
2015-05-01
We express the transport coefficients appearing in the second-order evolution equations for bulk viscous pressure and shear stress tensor using Bose-Einstein, Boltzmann, and Fermi-Dirac statistics for the equilibrium distribution function and Grad's 14-moment approximation as well as the method of Chapman-Enskog expansion for the nonequilibrium part. Focusing on the case of transversally homogeneous and boost-invariant longitudinal expansion of the viscous medium, we compare the results obtained using the above methods with those obtained from the exact solution of the massive 0+1 d relativistic Boltzmann equation in the relaxation-time approximation. We show that compared to the 14-moment approximation, the hydrodynamic transport coefficients obtained by employing the Chapman-Enskog method lead to better agreement with the exact solution of the relativistic Boltzmann equation.
Second order optical nonlinearity in silicon by symmetry breaking
NASA Astrophysics Data System (ADS)
Cazzanelli, Massimo; Schilling, Joerg
2016-03-01
Although silicon does not possess a dipolar bulk second order nonlinear susceptibility due to its centro-symmetric crystal structure, in recent years several attempts were undertaken to create such a property in silicon. This review presents the different sources of a second order susceptibility (χ(2)) in silicon and the connected second order nonlinear effects which were investigated up to now. After an introduction, a theoretical overview discusses the second order nonlinearity in general and distinguishes between the dipolar contribution—which is usually dominating in non-centrosymmetric structures—and the quadrupolar contribution, which even exists in centro-symmetric materials. Afterwards, the classic work on second harmonic generation from silicon surfaces in reflection measurements is reviewed. Due to the abrupt symmetry breaking at surfaces and interfaces locally a dipolar second order susceptibility appears, resulting in, e.g., second harmonic generation. Since the bulk contribution is usually small, the study of this second harmonic signal allows a sensitive observation of the surface/interface conditions. The impact of covering films, strain, electric fields, and defect states at the interfaces was already investigated in this way. With the advent of silicon photonics and the search for ever faster electrooptic modulators, the interest turned to the creation of a dipolar bulk χ(2) in silicon. These efforts have been focussing on several experiments applying an inhomogeneous strain to the silicon lattice to break its centro-symmetry. Recent results suggesting the impact of electric fields which are exerted from fixed charges in adjacent covering layers are also included. After a subsequent summary on "competing" concepts using not Si but Si-related materials, the paper will end with some final conclusions, suggesting possible future research direction in this dynamically developing field.
Robust stability of second-order systems
NASA Technical Reports Server (NTRS)
Chuang, C.-H.
1995-01-01
It has been shown recently how virtual passive controllers can be designed for second-order dynamic systems to achieve robust stability. The virtual controllers were visualized as systems made up of spring, mass and damping elements. In this paper, a new approach emphasizing on the notion of positive realness to the same second-order dynamic systems is used. Necessary and sufficient conditions for positive realness are presented for scalar spring-mass-dashpot systems. For multi-input multi-output systems, we show how a mass-spring-dashpot system can be made positive real by properly choosing its output variables. In particular, sufficient conditions are shown for the system without output velocity. Furthermore, if velocity cannot be measured then the system parameters must be precise to keep the system positive real. In practice, system parameters are not always constant and cannot be measured precisely. Therefore, in order to be useful positive real systems must be robust to some degrees. This can be achieved with the design presented in this paper.
First- and second-order Poisson spots
NASA Astrophysics Data System (ADS)
Kelly, William R.; Shirley, Eric L.; Migdall, Alan L.; Polyakov, Sergey V.; Hendrix, Kurt
2009-08-01
Although Thomas Young is generally given credit for being the first to provide evidence against Newton's corpuscular theory of light, it was Augustin Fresnel who first stated the modern theory of diffraction. We review the history surrounding Fresnel's 1818 paper and the role of the Poisson spot in the associated controversy. We next discuss the boundary-diffraction-wave approach to calculating diffraction effects and show how it can reduce the complexity of calculating diffraction patterns. We briefly discuss a generalization of this approach that reduces the dimensionality of integrals needed to calculate the complete diffraction pattern of any order diffraction effect. We repeat earlier demonstrations of the conventional Poisson spot and discuss an experimental setup for demonstrating an analogous phenomenon that we call a "second-order Poisson spot." Several features of the diffraction pattern can be explained simply by considering the path lengths of singly and doubly bent paths and distinguishing between first- and second-order diffraction effects related to such paths, respectively.
Synchronization from Second Order Network Connectivity Statistics
Zhao, Liqiong; Beverlin, Bryce; Netoff, Theoden; Nykamp, Duane Q.
2011-01-01
We investigate how network structure can influence the tendency for a neuronal network to synchronize, or its synchronizability, independent of the dynamical model for each neuron. The synchrony analysis takes advantage of the framework of second order networks, which defines four second order connectivity statistics based on the relative frequency of two-connection network motifs. The analysis identifies two of these statistics, convergent connections, and chain connections, as highly influencing the synchrony. Simulations verify that synchrony decreases with the frequency of convergent connections and increases with the frequency of chain connections. These trends persist with simulations of multiple models for the neuron dynamics and for different types of networks. Surprisingly, divergent connections, which determine the fraction of shared inputs, do not strongly influence the synchrony. The critical role of chains, rather than divergent connections, in influencing synchrony can be explained by their increasing the effective coupling strength. The decrease of synchrony with convergent connections is primarily due to the resulting heterogeneity in firing rates. PMID:21779239
Flavour fields in steady state: stress tensor and free energy
NASA Astrophysics Data System (ADS)
Banerjee, Avik; Kundu, Arnab; Kundu, Sandipan
2016-02-01
The dynamics of a probe brane in a given gravitational background is governed by the Dirac-Born-Infeld action. The corresponding open string metric arises naturally in studying the fluctuations on the probe. In Gauge-String duality, it is known that in the presence of a constant electric field on the worldvolume of the probe, the open string metric acquires an event horizon and therefore the fluctuation modes on the probe experience an effective temperature. In this article, we bring together various properties of such a system to a formal definition and a subsequent narration of the effective thermodynamics and the stress tensor of the corresponding flavour fields, also including a non-vanishing chemical potential. In doing so, we point out a potentially infinitely-degenerate scheme-dependence of regularizing the free energy, which nevertheless yields a universal contribution in certain cases. This universal piece appears as the coefficient of a log-divergence in free energy when a space-filling probe brane is embedded in AdS d+1-background, for d = 2, 4, and is related to conformal anomaly. For the special case of d = 2, the universal factor has a striking resemblance to the well-known heat current formula in (1 + 1)-dimensional conformal field theory in steady-state, which endows a plausible physical interpretation to it. Interestingly, we observe a vanishing conformal anomaly in d = 6.
NASA Technical Reports Server (NTRS)
Kiehn, R. M.
1976-01-01
With respect to irreversible, non-homeomorphic maps, contravariant and covariant tensor fields have distinctly natural covariance and transformational behavior. For thermodynamic processes which are non-adiabatic, the fact that the process cannot be represented by a homeomorphic map emphasizes the logical arrow of time, an idea which encompasses a principle of retrodictive determinism for covariant tensor fields.
The Formula of Grangeat for Tensor Fields of Arbitrary Order in n Dimensions
Schuster, T.
2007-01-01
The cone beam transform of a tensor field of order m in n ≥ 2 dimensions is considered. We prove that the image of a tensor field under this transform is related to a derivative of the n-dimensional Radon transform applied to a projection of the tensor field. Actually the relation we show reduces for m = 0 and n = 3 to the well-known formula of Grangeat. In that sense, the paper contains a generalization of Grangeat's formula to arbitrary tensor fields in any dimension. We further briefly explain the importance of that formula for the problem of tensor field tomography. Unfortunately, for m > 0, an inversion method cannot be derived immediately. Thus, we point out the possibility to calculate reconstruction kernels for the cone beam transform using Grangeat's formula. PMID:17713588
Feature-based interpolation of diffusion tensor fields and application to human cardiac DT-MRI.
Yang, Feng; Zhu, Yue-Min; Magnin, Isabelle E; Luo, Jian-Hua; Croisille, Pierre; Kingsley, Peter B
2012-02-01
Diffusion tensor interpolation is an important issue in the application of diffusion tensor magnetic resonance imaging (DT-MRI) to the human heart, all the more as the points representing the myocardium of the heart are often sparse. We propose a feature-based interpolation framework for the tensor fields from cardiac DT-MRI, by taking into account inherent relationships between tensor components. In this framework, the interpolation consists in representing a diffusion tensor in terms of two tensor features, eigenvalues and orientation, interpolating the Euler angles or the quaternion relative to tensor orientation and the logarithmically transformed eigenvalues, and reconstructing the tensor to be interpolated from the interpolated eigenvalues and tensor orientations. The results obtained with the aid of both synthetic and real cardiac DT-MRI data demonstrate that the feature-based schemes based on Euler angles or quaternions not only maintain the advantages of Log-Euclidean and Riemannian interpolation as for preserving the tensor's symmetric positive-definiteness and the monotonic determinant variation, but also preserve, at the same time, the monotonicity of fractional anisotropy (FA) and mean diffusivity (MD) values, which is not the case with Euclidean, Cholesky and Log-Euclidean methods. As a result, both interpolation schemes remove the phenomenon of FA collapse, and consequently avoid introducing artificial fiber crossing, with the difference that the quaternion is independent of coordinate system while Euler angles have the property of being more suitable for sophisticated interpolations. PMID:22154961
A second order parameter for 3SAT
Sandholm, T.W.
1996-12-31
The 3-satisfiability problem (3SAT) has had a central role in the study of complexity. It was recently found that 3SAT instances transition sharply from satisfiable to nonsatisfiable as the ratio of clauses to variables increases. Because this phase transition is so sharp, the ratio - an order parameter - can be used to predict satisfiability. This paper describes a second order parameter for 3SAT. Like the classical order parameter, it can be computed in linear time, but it analyzes the structure of the problem instance more deeply. We present an analytical method for using this new order parameter in conjunction with the classical one to enhance satisfiability prediction accuracy. The assumptions of the method are verified by rigorous statistical testing. The method significantly increases the satisfiability prediction accuracy over using the classical order parameter alone. Hardness - i.e. how long it takes to determine satisfiability - results for one complete and one incomplete algorithm from the literature are also presented as a function of the two order parameters. The importance of new order parameters lies in the fact that they refine the locating of satisfiable vs. nonsatisfiable and hard vs. easy formulas in the space of all problem instances by adding a new dimension in the analysis.
First- and second-order charged particle optics
Brown, K.L.; Servranckx, R.V.
1984-07-01
Since the invention of the alternating gradient principle there has been a rapid evolution of the mathematics and physics techniques applicable to charged particle optics. In this publication we derive a differential equation and a matrix algebra formalism valid to second-order to present the basic principles governing the design of charged particle beam transport systems. A notation first introduced by John Streib is used to convey the essential principles dictating the design of such beam transport systems. For example the momentum dispersion, the momentum resolution, and all second-order aberrations are expressed as simple integrals of the first-order trajectories (matrix elements) and of the magnetic field parameters (multipole components) characterizing the system. 16 references, 30 figures.
Deformable Registration of Feature-Endowed Point Sets Based on Tensor Fields.
Wassermann, Demian; Ross, James; Washko, George; Wells, William M; San Jose-Estepar, Raul
2014-06-01
The main contribution of this work is a framework to register anatomical structures characterized as a point set where each point has an associated symmetric matrix. These matrices can represent problem-dependent characteristics of the registered structure. For example, in airways, matrices can represent the orientation and thickness of the structure. Our framework relies on a dense tensor field representation which we implement sparsely as a kernel mixture of tensor fields. We equip the space of tensor fields with a norm that serves as a similarity measure. To calculate the optimal transformation between two structures we minimize this measure using an analytical gradient for the similarity measure and the deformation field, which we restrict to be a diffeomorphism. We illustrate the value of our tensor field model by comparing our results with scalar and vector field based models. Finally, we evaluate our registration algorithm on synthetic data sets and validate our approach on manually annotated airway trees. PMID:25473253
Deformable Registration of Feature-Endowed Point Sets Based on Tensor Fields
Wassermann, Demian; Ross, James; Washko, George; Wells, William M.; San Jose-Estepar, Raul
2014-01-01
The main contribution of this work is a framework to register anatomical structures characterized as a point set where each point has an associated symmetric matrix. These matrices can represent problem-dependent characteristics of the registered structure. For example, in airways, matrices can represent the orientation and thickness of the structure. Our framework relies on a dense tensor field representation which we implement sparsely as a kernel mixture of tensor fields. We equip the space of tensor fields with a norm that serves as a similarity measure. To calculate the optimal transformation between two structures we minimize this measure using an analytical gradient for the similarity measure and the deformation field, which we restrict to be a diffeomorphism. We illustrate the value of our tensor field model by comparing our results with scalar and vector field based models. Finally, we evaluate our registration algorithm on synthetic data sets and validate our approach on manually annotated airway trees. PMID:25473253
Non-Gaussianity from the second-order cosmological perturbation
NASA Astrophysics Data System (ADS)
Lyth, David H.; Rodríguez, Yeinzon
2005-06-01
Several conserved and/or gauge-invariant quantities described as the second-order curvature perturbation have been given in the literature. We revisit various scenarios for the generation of second-order non-Gaussianity in the primordial curvature perturbation ζ, employing for the first time a unified notation and focusing on the normalization fNL of the bispectrum. When ζ first appears a few Hubble times after horizon exit, |fNL| is much less than 1 and is, therefore, negligible. Thereafter ζ (and hence fNL) is conserved as long as the pressure is a unique function of energy density (adiabatic pressure). Nonadiabatic pressure comes presumably only from the effect of fields, other than the one pointing along the inflationary trajectory, which are light during inflation (“light noninflaton fields”). During single-component inflation fNL is constant, but multicomponent inflation might generate |fNL|˜1 or bigger. Preheating can affect fNL only in atypical scenarios where it involves light noninflaton fields. The simplest curvaton scenario typically gives fNL≪-1 or fNL=+5/4. The inhomogeneous reheating scenario can give a wide range of values for fNL. Unless there is a detection, observation can eventually provide a limit |fNL|≲1, at which level it will be crucial to calculate the precise observational limit using second-order theory.
Unifying the PST and the auxiliary tensor field formulations of 4D self-duality
NASA Astrophysics Data System (ADS)
Ivanov, E. A.; Nurmagambetov, A. J.; Zupnik, B. M.
2014-04-01
We unify the Lorentz- and O(2) duality-covariant approach to 4D self-dual theories by Pasti, Sorokin and Tonin (PST) with the formulation involving an auxiliary tensor field. We present the basic features of the new hybrid approach, including symmetries of the relevant generalized PST action. Its salient peculiarity is the unique form of the realization of the PST gauge symmetries. The corresponding transformations do not affect the auxiliary tensor field, which guarantees the self-duality of the nonlinear actions in which the O(2) invariant interactions are constructed out of the tensor field.
Overconnections and the energy-tensors of gauge and gravitational fields
NASA Astrophysics Data System (ADS)
Canarutto, Daniel
2016-08-01
A geometric construction for obtaining a prolongation of a connection to a connection of a bundle of connections is presented. This determines a natural extension of the notion of canonical energy-tensor which suits gauge and gravitational fields, and shares the main properties of the energy-tensor of a matter field in the jet space formulation of Lagrangian field theory, in particular with regards to symmetries of the Poincaré-Cartan form. Accordingly, the joint energy-tensor for interacting matter and gauge fields turns out to be a natural geometric object, whose definition needs no auxiliary structures. Various topics related to energy-tensors, symmetries and the Einstein equations in a theory with interacting matter, gauge and gravitational fields can be viewed under a clarifying light. Finally, the symmetry determined by the "Komar superpotential" is expressed as a symmetry of the gravitational Poincaré-Cartan form.
Spectral expansions of homogeneous and isotropic tensor-valued random fields
NASA Astrophysics Data System (ADS)
Malyarenko, Anatoliy; Ostoja-Starzewski, Martin
2016-06-01
We establish spectral expansions of tensor-valued homogeneous and isotropic random fields in terms of stochastic integrals with respect to orthogonal scattered random measures previously known only for the case of tensor rank 0. The fields under consideration take values in the 3-dimensional Euclidean space {E^3} and in the space {S^2(E^3)} of symmetric rank 2 tensors over {E^3}. We find a link between the theory of random fields and the theory of finite-dimensional convex compact sets. These random fields furnish stepping-stone for models of rank 1 and rank 2 tensor-valued fields in continuum physics, such as displacement, velocity, stress, strain, providing appropriate conditions (such as the governing equation or positive-definiteness) are imposed.
Dressed four-wave mixing second-order Talbot effect
NASA Astrophysics Data System (ADS)
Chen, Haixia; Zhang, Xun; Zhu, Dayu; Yang, Chang; Jiang, Tao; Zheng, Huaibin; Zhang, Yanpeng
2014-10-01
We theoretically demonstrate second-order Talbot effect (SOTE) based on entangled photon pairs. The photon pairs are generated from the spontaneous parametric four-wave mixing (SPFWM) process in a cold atomic medium and can be taken as the imaging light in order to realize coincidence recording. A strong standing wave is used to create the electromagnetically induced grating in the entangled photon pairs channels. By changing the frequency detuning of the standing wave or the other optical fields participating in the process, we can manipulate the contrast of the second-order Talbot image. We use the second-order correlation function and the dressed-state picture to explain the SOTE occurring in the SPFWM process. Moreover, we demonstrate the scheme for SOTE based on the spatially correlated twin beams generated from the SPFWM process with injection. This scheme provides a convenient detection proposal for the SOTE at the cost of the image contrast. Compared to the previous self-imaging schemes, the present schemes have the characteristic of controllable image contrast and of nonlocal imaging, and thus, they might broaden their applications in imaging techniques and find applications in quantum lithography.
Conservation laws and stress-energy-momentum tensors for systems with background fields
Gratus, Jonathan; Obukhov, Yuri N.; Tucker, Robin W.
2012-10-15
This article attempts to delineate the roles played by non-dynamical background structures and Killing symmetries in the construction of stress-energy-momentum tensors generated from a diffeomorphism invariant action density. An intrinsic coordinate independent approach puts into perspective a number of spurious arguments that have historically lead to the main contenders, viz the Belinfante-Rosenfeld stress-energy-momentum tensor derived from a Noether current and the Einstein-Hilbert stress-energy-momentum tensor derived in the context of Einstein's theory of general relativity. Emphasis is placed on the role played by non-dynamical background (phenomenological) structures that discriminate between properties of these tensors particularly in the context of electrodynamics in media. These tensors are used to construct conservation laws in the presence of Killing Lie-symmetric background fields. - Highlights: Black-Right-Pointing-Pointer The role of background fields in diffeomorphism invariant actions is demonstrated. Black-Right-Pointing-Pointer Interrelations between different stress-energy-momentum tensors are emphasised. Black-Right-Pointing-Pointer The Abraham and Minkowski electromagnetic tensors are discussed in this context. Black-Right-Pointing-Pointer Conservation laws in the presence of nondynamic background fields are formulated. Black-Right-Pointing-Pointer The discussion is facilitated by the development of a new variational calculus.
Slowly rotating neutron stars in scalar-tensor theories with a massive scalar field
NASA Astrophysics Data System (ADS)
Yazadjiev, Stoytcho S.; Doneva, Daniela D.; Popchev, Dimitar
2016-04-01
In the scalar-tensor theories with a massive scalar field, the coupling constants, and the coupling functions in general, which are observationally allowed, can differ significantly from those in the massless case. This fact naturally implies that the scalar-tensor neutron stars with a massive scalar field can have rather different structure and properties in comparison with their counterparts in the massless case and in general relativity. In the present paper, we study slowly rotating neutron stars in scalar-tensor theories with a massive gravitational scalar. Two examples of scalar-tensor theories are examined—the first example is the massive Brans-Dicke theory and the second one is a massive scalar-tensor theory indistinguishable from general relativity in the weak-field limit. In the latter case, we study the effect of the scalar field mass on the spontaneous scalarization of neutron stars. Our numerical results show that the inclusion of a mass term for the scalar field indeed changes the picture drastically compared to the massless case. It turns out that mass, radius, and moment of inertia for neutron stars in massive scalar-tensor theories can differ drastically from the pure general relativistic solutions if sufficiently large masses of the scalar field are considered.
Entanglement in a second-order quantum phase transition
Vidal, Julien; Palacios, Guillaume; Mosseri, Remy
2004-02-01
We consider a system of mutually interacting spins 1/2 embedded in a transverse magnetic field which undergoes a second-order quantum phase transition. We analyze the entanglement properties and the spin squeezing of the ground state and show that, contrarily to the one-dimensional case, a cusplike singularity appears at the critical point {lambda}{sub c} in the thermodynamical limit. We also show that there exists a value {lambda}{sub 0}{>=}{lambda}{sub c} above which the ground state is not spin squeezed despite a nonvanishing concurrence.
A Finite Field Method for Calculating Molecular Polarizability Tensors for Arbitrary Multipole Rank
Elking, Dennis M.; Perera, Lalith; Duke, Robert; Darden, Thomas; Pedersen, Lee G.
2011-01-01
A finite field method for calculating spherical tensor molecular polarizability tensors αlm;l′m′ = ∂Δlm/∂ϕl′m′* by numerical derivatives of induced molecular multipole Δlm with respect to gradients of electrostatic potential ϕl′m′* is described for arbitrary multipole ranks l and l′. Inter-conversion formulae for transforming multipole moments and polarizability tensors between spherical and traceless Cartesian tensor conventions are derived. As an example, molecular polarizability tensors up to the hexadecapole-hexadecapole level are calculated for water at the HF, B3LYP, MP2, and CCSD levels. In addition, inter-molecular electrostatic and polarization energies calculated by molecular multipoles and polarizability tensors are compared to ab initio reference values calculated by the Reduced Variation Space (RVS) method for several randomly oriented small molecule dimers separated by a large distance. It is discussed how higher order molecular polarizability tensors can be used as a tool for testing and developing new polarization models for future force fields. PMID:21915883
Characterization of Tricoordinate Boron Chemical Shift Tensors: Definitive High-Field
Bryce, David L.; Wasylishen, Roderick E.; Gee, Myrlene
2001-01-01
Despite the large known chemical shift (CS) range for boron and the large number of 11B NMR studies of glasses, no boron CS tensors have been characterized to date. We report the application of solid-state NMR techniques at moderate (9.4 T) and high (17.63 T) applied magnetic field strengths to the characterization of the boron CS tensors in trimesitylborane (BMes3) and triphenyl borate (B(OPh)3). The boron CS tensor of the former compound exhibits a remarkably large span,? 121 1 ppm, which encompasses the known range of isotropic chemical shifts for tricoordinate boron compounds. Conversely, the effect of the boron CS tensor on the 11B NMR spectra of B(OPh)3 is difficult to observe and quantify even at field strengths as high
Vacuum stress-energy tensor of a massive scalar field in a wormhole spacetime
Bezerra, V. B.; Bezerra de Mello, E. R.; Khusnutdinov, N. R.; Sushkov, S. V.
2010-04-15
The vacuum average value of the stress-energy tensor of a massive scalar field with nonminimal coupling {xi} to the curvature on the short-throat flat-space wormhole background is calculated. The final analysis is made numerically. It was shown that the energy-momentum tensor does not violate the null energy condition near the throat. Therefore, the vacuum polarization cannot self-consistently support the wormhole.
Relativistic second-order dissipative fluid dynamics at finite chemical potential
NASA Astrophysics Data System (ADS)
Jaiswal, Amaresh; Friman, Bengt; Redlich, Krzysztof
2016-07-01
We employ a Chapman-Enskog like expansion for the distribution function close to equilibrium to solve the Boltzmann equation in the relaxation time approximation and subsequently derive second-order evolution equations for dissipative charge currentand shear stress tensor for a system of massless quarks and gluons. We use quantum statistics for the phase space distribution functions to calculate the transport coefficients. We show that, the second-order evolution equations for the dissipative charge current and the shear stress tensor can be decoupled. We find that, for large chemical potential, the charge conductivity is small compared to the shear viscosity. Moreover, we demonstrate that the limiting behaviour of the ratio of heat conductivity to shear viscosity is identicalto that obtained for a strongly coupled conformal plasma.
A preliminary compressible second-order closure model for high speed flows
NASA Technical Reports Server (NTRS)
Speziale, Charles G.; Sarkar, Sutanu
1989-01-01
A preliminary version of a compressible second-order closure model that was developed in connection with the National Aero-Space Plane Project is presented. The model requires the solution of transport equations for the Favre-averaged Reynolds stress tensor and dissipation rate. Gradient transport hypotheses are used for the Reynolds heat flux, mass flux, and turbulent diffusion terms. Some brief remarks are made about the direction of future research to generalize the model.
Effenberger, F.; Fichtner, H.; Scherer, K.; Barra, S.; Kleimann, J.; Strauss, R. D.
2012-05-10
The spatial diffusion of cosmic rays in turbulent magnetic fields can, in the most general case, be fully anisotropic, i.e., one has to distinguish three diffusion axes in a local, field-aligned frame. We reexamine the transformation for the diffusion tensor from this local to a global frame, in which the Parker transport equation for energetic particles is usually formulated and solved. Particularly, we generalize the transformation formulae to allow for an explicit choice of two principal local perpendicular diffusion axes. This generalization includes the 'traditional' diffusion tensor in the special case of isotropic perpendicular diffusion. For the local frame, we describe the motivation for the choice of the Frenet-Serret trihedron, which is related to the intrinsic magnetic field geometry. We directly compare the old and the new tensor elements for two heliospheric magnetic field configurations, namely the hybrid Fisk and Parker fields. Subsequently, we examine the significance of the different formulations for the diffusion tensor in a standard three-dimensional model for the modulation of galactic protons. For this, we utilize a numerical code to evaluate a system of stochastic differential equations equivalent to the Parker transport equation and present the resulting modulated spectra. The computed differential fluxes based on the new tensor formulation deviate from those obtained with the 'traditional' one (only valid for isotropic perpendicular diffusion) by up to 60% for energies below a few hundred MeV depending on heliocentric distance.
A Topologically-Informed Hyperstreamline Seeding Method for Alignment Tensor Fields.
Fu, Fred; Abukhdeir, Nasser Mohieddin
2015-03-01
A topologically-informed hyperstreamline seeding method is presented for visualization of alignment tensor fields. The method is inspired by and applied to visualization of nematic liquid crystal (LC) orientation dynamics simulations. The method distributes hyperstreamlines along domain boundaries and edges of a nearest-neighbor graph whose vertices are degenerate regions of the alignment tensor field, which correspond to orientational defects in a nematic LC domain. This is accomplished without iteration while conforming to a user-specified spacing between hyperstreamlines and avoids possible failure modes associated with hyperstreamline integration in the vicinity of degeneracies in alignment (orientational defects). It is shown that the presented seeding method enables automated hyperstreamline-based visualization of a broad range of alignment tensor fields which enhances the ability of researchers to interpret these fields and provides an alternative to using glyph-based techniques. PMID:26357072
Second-Order Fermi Acceleration and Emission in Blazar Jets
NASA Astrophysics Data System (ADS)
Asano, Katsuaki; Takahara, Fumio; Toma, Kenji; Kusunose, Masaaki; Kakuwa, Jun
The second-order Fermi acceleration (Fermi-II) driven by turbulence may be responsible for the electron acceleration in blazar jets. We test this model with time-dependent simulations, adopt it for 1ES 1101-232, and Mrk 421. The Fermi-II model with radial evolution of the electron injection rate and/or diffusion coefficient can reproduce the spectra from the radio to the gamma-ray regime. For Mrk 421, an external radio photon field with a luminosity of 4.9 begin{math} {times} 10 (38) erg s (-1) is required to agree with the observed GeV flux. The temporal variability of the diffusion coefficient or injection rate causes flare emission. The observed synchronicity of X-ray and TeV flares implies a decrease of the magnetic field in the flaring source region.
Magnetic Compensation for Second-Order Doppler Shift in LITS
NASA Technical Reports Server (NTRS)
Burt, Eric; Tjoelker, Robert
2008-01-01
The uncertainty in the frequency of a linear-ion-trap frequency standard (LITS) can be reduced substantially by use of a very small magnetic inhomogeneity tailored to compensate for the residual second-order Doppler shift. An effect associated with the relativistic time dilatation, one cause of the second-order Doppler shift, is ion motion that is attributable to the trapping radio-frequency (RF)electromagnetic field used to trap ions. The second-order Doppler shift is reduced by using a multi-pole trap; however it is still the largest source of systematic frequency shift in the latest generation of LITSs, which are among the most stable clocks in the world. The present compensation scheme reduces the frequency instability of the affected LITS to about a tenth of its previous value. The basic principles of prior generation LITSs were discussed in several prior NASA Tech Briefs articles. Below are recapitulated only those items of basic information necessary to place the present development in context. A LITS includes a microwave local oscillator, the frequency of which is stabilized by comparison with the frequency of the ground state hyperfine transition of 199Hg+ ions. The comparison involves a combination of optical and microwave excitation and interrogation of the ions in a linear ion trap in the presence of a nominally uniform magnetic field. In the current version of the LITS, there are two connected traps (see figure): (1) a quadrupole trap wherein the optical excitation and measurement take place and (2) a 12-pole trap (denoted the resonance trap), wherein the microwave interrogation takes place. The ions are initially loaded into the quadrupole trap and are thereafter shuttled between the two traps. Shuttling ions into the resonance trap allows sensitive microwave interrogation to take place well away from loading interference. The axial magnetic field for the resonance trap is generated by an electric current in a finely wound wire coil surrounded by
Scalar-tensor gravity with a non-minimally coupled Higgs field and accelerating universe
NASA Astrophysics Data System (ADS)
Sim, Jonghyun; Lee, Tae Hoon
2016-03-01
We consider general couplings, including non-minimal derivative coupling, of a Higgs boson field to scalar-tensor gravity and calculate their contributions to the energy density and pressure in Friedmann-Robertson-Walker spacetime. In a special case where the kinetic term of the Higgs field is non-minimally coupled to the Einstein tensor, we seek de Sitter solutions for the cosmic scale factor and discuss the possibility that the late-time acceleration and the inflationary era of our universe can be described by means of scalar fields with self-interactions and the Yukawa potential.
Metric energy-momentum tensor for polarizable particles in an electormagnetic field
NASA Astrophysics Data System (ADS)
Maksimenko, N. V.; Lukashevich, S. A.
2006-12-01
Within the covariant Lagrange formalism and the relativistic theory of continuous media, the metric energy-momentum tensor is obtained for spin polarizable particles interacting with an electromagnetic field. An equation of motion of the polarizable particles with a spin of 1/2 in an external electromagnetic field is derived.
Second order closure modeling of turbulent buoyant wall plumes
NASA Technical Reports Server (NTRS)
Zhu, Gang; Lai, Ming-Chia; Shih, Tsan-Hsing
1992-01-01
Non-intrusive measurements of scalar and momentum transport in turbulent wall plumes, using a combined technique of laser Doppler anemometry and laser-induced fluorescence, has shown some interesting features not present in the free jet or plumes. First, buoyancy-generation of turbulence is shown to be important throughout the flow field. Combined with low-Reynolds-number turbulence and near-wall effect, this may raise the anisotropic turbulence structure beyond the prediction of eddy-viscosity models. Second, the transverse scalar fluxes do not correspond only to the mean scalar gradients, as would be expected from gradient-diffusion modeling. Third, higher-order velocity-scalar correlations which describe turbulent transport phenomena could not be predicted using simple turbulence models. A second-order closure simulation of turbulent adiabatic wall plumes, taking into account the recent progress in scalar transport, near-wall effect and buoyancy, is reported in the current study to compare with the non-intrusive measurements. In spite of the small velocity scale of the wall plumes, the results showed that low-Reynolds-number correction is not critically important to predict the adiabatic cases tested and cannot be applied beyond the maximum velocity location. The mean and turbulent velocity profiles are very closely predicted by the second-order closure models. but the scalar field is less satisfactory, with the scalar fluctuation level underpredicted. Strong intermittency of the low-Reynolds-number flow field is suspected of these discrepancies. The trends in second- and third-order velocity-scalar correlations, which describe turbulent transport phenomena, are also predicted in general, with the cross-streamwise correlations better than the streamwise one. Buoyancy terms modeling the pressure-correlation are shown to improve the prediction slightly. The effects of equilibrium time-scale ratio and boundary condition are also discussed.
Compact Two-State-Variable Second-Order Memristor Model.
Kim, Sungho; Kim, Hee-Dong; Choi, Sung-Jin
2016-06-01
A key requirement for using memristors in functional circuits is a predictive physical model to capture the resistive switching behavior, which shall be compact enough to be implemented using a circuit simulator. Although a number of memristor models have been developed, most of these models (i.e., first-order memristor models) have utilized only a one-state-variable. However, such simplification is not adequate for accurate modeling because multiple mechanisms are involved in resistive switching. Here, a two-state-variable based second-order memristor model is presented, which considers the axial drift of the charged vacancies in an applied electric field and the radial vacancy motion caused by the thermophoresis and diffusion. In particular, this model emulates the details of the intrinsic short-term dynamics, such as decay and temporal heat summation, and therefore, it accurately predicts the resistive switching characteristics for both DC and AC input signals. PMID:27152649
Second-Order Nonlinear Optical Imaging of Chiral Crystals
Kissick, David J.; Wanapun, Debbie; Simpson, Garth J.
2012-01-01
Second-order nonlinear optical imaging of chiral crystals (SONICC) is an emerging technique for crystal imaging and characterization. We provide a brief overview of the origin of second harmonic generation signals in SONICC and discuss recent studies using SONICC for biological applications. Given that they provide near-complete suppression of any background, SONICC images can be used to determine the presence or absence of protein crystals through both manual inspection and automated analysis. Because SONICC creates high-resolution images, nucleation and growth kinetics can also be observed. SONICC can detect metastable, homochiral crystalline forms of amino acids crystallizing from racemic solutions, which confirms Ostwald’s rule of stages for crystal growth. SONICC’s selectivity, based on order, and sensitivity, based on background suppression, make it a promising technique for numerous fields concerned with chiral crystal formation. PMID:21469954
Gravitational Microlensing by Ellis Wormhole: Second Order Effects
NASA Astrophysics Data System (ADS)
Lukmanova, Regina; Kulbakova, Aliya; Izmailov, Ramil; Potapov, Alexander A.
2016-07-01
Gravitational lensing is the effect of light bending in a gravitational field. It can be used as a possible observational method to detect or exclude the existence of wormholes. In this work, we extend the work by Abe on gravitational microlensing by Ellis wormhole by including the second order deflection term. Using the lens equation and definition of Einstein radius, we find the angular locations of the physical image inside and outside Einstein ring. The work contains a comparative analysis of light curves between the Schwarzschild black hole and the Ellis wormhole that can be used to distinguish such objects though such distinctions are too minute to be observable even in the near future. We also tabulate the optical depth and event rate for lensing by bulge and Large Magellanic Cloud (LMC) stars.
Effect of carbazole as a donor moiety on the second-order nonlinearity of organic molecules
NASA Astrophysics Data System (ADS)
Meshulam, Guilia; Berkovic, Garry; Kotler, Zvi; Ben-Asuly, Amos; Mazor, Royi; Shapiro, Lev; Khodorkovsky, Vladimir
1999-10-01
The second order nonlinearity of conjugated organic molecules involving, 1,3 indandione derivatives as an acceptor moiety has been studied. Varying the donor from dialkylamino to the chemically similar substituent, N- carbazolyl resulted in a drastic reduction of electric field induced second harmonic (beta) values. For some molecules, even a small negative value of (beta) was received. Quantum chemical calculations indicate that the decrease occurs as a result of two overlapping transitions, which contribute to (beta) with opposite signs. The charge transfer band gives a positive (beta) zzz along the molecular long axis, while a transition essentially within the carbazole moiety provides a negative (beta zzz contribution to (beta EFISH. Thus, these molecules must be described with a 2D model as opposed to the 'classical' model of 1D nonlinear optical chromophores. The prediction of the 2D model was verified experimentally by using a combination of two methods, EFISH and Hyper-Rayleigh Scattering, which probe different combination of the (beta) tensor elements.
NASA Astrophysics Data System (ADS)
Reina, Borja; Vera, Raül
2015-08-01
Hartle's model describes the equilibrium configuration of a rotating isolated compact body in perturbation theory up to second order in general relativity. The interior of the body is a perfect fluid with a barotropic equation of state, no convective motions and rigid rotation. That interior is matched across its surface to an asymptotically flat vacuum exterior. Perturbations are taken to second order around a static and spherically symmetric background configuration. Apart from the explicit assumptions, the perturbed configuration is constructed upon some implicit premises, in particular the continuity of the functions describing the perturbation in terms of some background radial coordinate. In this work we revisit the model within a modern general and consistent theory of perturbative matchings to second order, which is independent of the coordinates and gauges used to describe the two regions to be joined. We explore the matching conditions up to second order in full. The main particular result we present is that the radial function m0 (in the setting of the original work) of the second order perturbation tensor, contrary to the original assumption, presents a jump at the surface of the star, which is proportional to the value of the energy density of the background configuration there. As a consequence, the change in mass δ M needed by the perturbed configuration to keep the value of the central energy density unchanged must be amended. We also discuss some subtleties that arise when studying the deformation of the star.
Design of optimal second-order state estimators
NASA Technical Reports Server (NTRS)
Joshi, Suresh M.
1991-01-01
The present consideration of the design of online computation-saving second-order state estimators for second-order vector-matrix differential systems proposes a class of such estimators which is proven to possess guaranteed convergence. A class of optimal second-order estimators is then obtained, and the conditions required for optimality are identified. The estimator proposed offers high performance in conjunction with online computation reductions sufficiently great to allow the estimation of the large number of state variables associated with control of large, flexible space structures represented by high-dimensional second-order systems.
The Morse Oscillator and Second-Order Perturbation Theory
NASA Astrophysics Data System (ADS)
Pettitt, B. A.
1998-09-01
This article shows how the energies of the Morse oscillator are obtained exactly from a second-order perturbation expansion in a harmonic oscillator basis. This exercise is recommended for its instructional value in intermediate quantum chemistry, in that the second-order term is entirely tractable, it arises within an important context (anharmonicity of vibrations), and it gives the right answer.
An analysis of the intermediate field theory of T4 tensor model
NASA Astrophysics Data System (ADS)
Nguyen, Viet Anh; Dartois, Stéphane; Eynard, Bertrand
2015-01-01
In this paper we analyze the multi-matrix model arising from the intermediate field representation of the tensor model with all quartic melonic interactions. We derive the saddle point equation and the Schwinger-Dyson constraints. We then use them to describe the leading and next-to-leading eigenvalues distribution of the matrices.
A planar second-order DC SQUID gradiometer.
Carelli, P; Chiaventi, L; Leoni, R; Pullano, M; Schirripa Spagnolo, G
1991-01-01
In this work we describe a DC SQUID gradiometer, sensitive to the second spatial derivative of the magnetic field. The sensitive area of the gradiometer is the inductive body of the DC SQUID itself. The isoflux line distribution generated by a dipolar source, obtained by performing magnetic measurements with an array of such detectors, is relatively complicated, but its localisation capability is similar to that one usually achieves with axial detector arrays. Planar gradiometers also show a better resolution for near sources and a stronger rejection of far disturbances. The final device is expected to have an inductance of a few hundreds of pH in order to obtain performances typical of a low noise DC SQUID. The pick-up coils will be the combination of four square holes of 500 microns side with a 1.05 cm baseline. Due to the magnetic field concentration (in the final device it can be a factor 10) the gradiometer will have a sensitivity of 10(-11) T m-2 Hz-1/2 and a field sensitivity of about 2 fT Hz-1/2. Some preliminary results, obtained on detectors with an intermediate area between the prototype and final device, are reported here. The process used to fabricate this second-order gradiometer is based on Nb-NbO chi-PbAuIn Josephson tunnel junctions. Some possible improvements will also be described. PMID:1807874
Edge detection of gravity field using eigenvalue analysis of gravity gradient tensor
NASA Astrophysics Data System (ADS)
Zuo, Boxin; Hu, Xiangyun
2015-03-01
In this paper, eigenvalues of the full gravity gradient tensor (GGT) are used to detect edges of geological structure. First, the solving of GGT eigenvalues is discussed; then a new edge detection method is proposed by using the eigenvalues of GGT. Comparing with the pervious edge detection method based on curvature gravity gradient tensor (CGGT), the full gravity gradient tensor contains more independent gradient components that are helpful to detect more subtle structures of the sources. The proposed method is applied to the synthetic data with and without noise to determine the locations of the edges of the mixed positive/negative contract density bodies. It has also been tested on real field data. All of the experimental results have shown that the newly proposed method is effective for edge detection.
NASA Astrophysics Data System (ADS)
Yang, Jianfei; Poot, Dirk H. J.; Arkesteijn, Georgius A. M.; Caan, Matthan W.; van Vliet, Lucas J.; Vos, Frans M.
2015-03-01
Conventionally, a single rank-2 tensor is used to assess the white matter integrity in diffusion imaging of the human brain. However, a single tensor fails to describe the diffusion in fiber crossings. Although a dual tensor model is able to do so, the low signal-to-noise ratio hampers reliable parameter estimation as the number of parameters is doubled. We present a framework for structure-adaptive tensor field filtering to enhance the statistical analysis in complex fiber structures. In our framework, a tensor model will be fitted based on an automated relevance determination method. Particularly, a single tensor model is applied to voxels in which the data seems to represent a single fiber and a dualtensor model to voxels appearing to contain crossing fibers. To improve the estimation of the model parameters we propose a structure-adaptive tensor filter that is applied to tensors belonging to the same fiber compartment only. It is demonstrated that the structure-adaptive tensor-field filter improves the continuity and regularity of the estimated tensor field. It outperforms an existing denoising approach called LMMSE, which is applied to the diffusion-weighted images. Track-based spatial statistics analysis of fiber-specific FA maps show that the method sustains the detection of more subtle changes in white matter tracts than the classical single-tensor-based analysis. Thus, the filter enhances the applicability of the dual-tensor model in diffusion imaging research. Specifically, the reliable estimation of two tensor diffusion properties facilitates fiber-specific extraction of diffusion features.
A second order accurate embedded boundary method for the wave equation with Dirichlet data
Kreiss, H O; Petersson, N A
2004-03-02
The accuracy of Cartesian embedded boundary methods for the second order wave equation in general two-dimensional domains subject to Dirichlet boundary conditions is analyzed. Based on the analysis, we develop a numerical method where both the solution and its gradient are second order accurate. We avoid the small-cell stiffness problem without sacrificing the second order accuracy by adding a small artificial term to the Dirichlet boundary condition. Long-time stability of the method is obtained by adding a small fourth order dissipative term. Several numerical examples are provided to demonstrate the accuracy and stability of the method. The method is also used to solve the two-dimensional TM{sub z} problem for Maxwell's equations posed as a second order wave equation for the electric field coupled to ordinary differential equations for the magnetic field.
Large tensor mode, field range bound and consistency in generalized G-inflation
NASA Astrophysics Data System (ADS)
Kunimitsu, Taro; Suyama, Teruaki; Watanabe, Yuki; Yokoyama, Jun'ichi
2015-08-01
We systematically show that in potential driven generalized G-inflation models, quantum corrections coming from new physics at the strong coupling scale can be avoided, while producing observable tensor modes. The effective action can be approximated by the tree level action, and as a result, these models are internally consistent, despite the fact that we introduced new mass scales below the energy scale of inflation. Although observable tensor modes are produced with sub-strong coupling scale field excursions, this is not an evasion of the Lyth bound, since the models include higher-derivative non-canonical kinetic terms, and effective rescaling of the field would result in super-Planckian field excursions. We argue that the enhanced kinetic term of the inflaton screens the interactions with other fields, keeping the system weakly coupled during inflation.
Study of second order upwind differencing in a recirculating flow
NASA Technical Reports Server (NTRS)
Vanka, S. P.
1985-01-01
The accuracy and stability of the second order upwind differencing scheme was investigated. The solution algorithm employed is based on a coupled solution of the nonlinear finite difference equations by the multigrid technique. Calculations have been made of the driven cavity flow for several Reynolds numbers and finite difference grids. In comparison with the hybrid differencing, the second order upwind differencing is somewhat more accurate but it is not monotonically accurate with mesh refinement. Also, the convergence of the solution algorithm deteriorates with the use of the second order upwind differencing.
A study of second-order supersonic flow theory
NASA Technical Reports Server (NTRS)
Van Dyke, Milton D
1952-01-01
Second-order solutions of supersonic-flow problems are sought by iteration, using the linearized solution as the first step. For plane and axially symmetric flows, particular solutions of the iteration equation are discovered which reduce the second-order problem to an equivalent linearized problem. Comparison of second-order solutions with exact and numerical results shows great improvement over linearized theory. For full three-dimensional flow, only a partial particular solution is found. The inclined cone is solved, and the possibility of treating more general problems is considered.
Method to render second order beam optics programs symplectic
Douglas, D.; Servranckx, R.V.
1984-10-01
We present evidence that second order matrix-based beam optics programs violate the symplectic condition. A simple method to avoid this difficulty, based on a generating function approach to evaluating transfer maps, is described. A simple example illustrating the non-symplectricity of second order matrix methods, and the effectiveness of our solution to the problem, is provided. We conclude that it is in fact possible to bring second order matrix optics methods to a canonical form. The procedure for doing so has been implemented in the program DIMAT, and could be implemented in programs such as TRANSPORT and TURTLE, making them useful in multiturn applications. 15 refs.
Second-order state estimation experiments using acceleration measurements
NASA Technical Reports Server (NTRS)
Belvin, W. K.
1992-01-01
The estimation of dynamic states for feedback control of structural systems using second-order differential equations and acceleration measurements is described. The formulation of the observer model, and the design of the observer gains is discussed in detail. It is shown the second-order observer is highly stable because the stability constraints on the observer gains are model independent. The limitation of the proposed observer is the need for 'nearly' collocated actuators and accelerometers. Experimental results using a control-structure interaction testbed are presented that show the second-order observer provided more stability than a Kalman filter estimator without decreasing closed-loop performance.
{open_quotes}Quadrupoled{close_quotes} materials for second-order nonlinear optics
Hubbard, S.F.; Petschek, R.G.; Singer, K.D.
1997-10-01
We describe a new approach to second-order nonlinear optical materials, namely quadrupoling. This approach is valid in the regime of Kleinman (full permutation) symmetry breaking, and thus requires a two- or three dimensional microscopic nonlinearity at wavelengths away from material resonances. This {open_quotes}quadrupolar{close_quotes} nonlinearity arises from the second rank pseudotensor of the rotationally invariant representation of the second-order nonlinear optical tensor. We have experimentally investigated candidate molecules comprised of chiral camphorquinone derivatives by measuring the scalar invariant associated with the rank two pseudotensor using hyper-Rayleigh scattering. We have found sizable scalar figures of merit for several compounds using light for which the second harmonic wavelengths are greater than 100 nm longer than the absorption peak location. At these wavelengths, the quadrupolar scalar is as large as the polar (EFISH) scalar of p-nitroaniline. Prospects for applications are discussed.
Orthogonal canonical forms for second-order systems
NASA Technical Reports Server (NTRS)
Williams, Trevor; Laub, Alan
1989-01-01
The authors prove that a linear second-order system with arbitrary damping cannot be reduced to Hessenberg-triangular form by means of orthogonal transformations, while this reduction is always possible for the modal damping commonly assumed for models of flexible structures. The type of canonical form obtainable by means of orthogonal transformations acting on a second-order system is heavily dependent on the type of damping considered. If the damping matrix is merely positive semi-definite symmetric, it is generally not possible to obtain a reduction to Hessenberg-triangular form, while this reduction is trivial for zero or Rayleigh damping. If damping is modal, however, as is commonly assumed in structural models, the reduction exists and is nontrivial. Furthermore, reduction to triangular second-order Schur form is always possible for such damping: this canonical form appears likely to have applications to second-order system theory.
Oscillation theorems for second order nonlinear forced differential equations.
Salhin, Ambarka A; Din, Ummul Khair Salma; Ahmad, Rokiah Rozita; Noorani, Mohd Salmi Md
2014-01-01
In this paper, a class of second order forced nonlinear differential equation is considered and several new oscillation theorems are obtained. Our results generalize and improve those known ones in the literature. PMID:25077054
Fast second-order consensus via predictive mechanisms
NASA Astrophysics Data System (ADS)
Wu, Jie; Zhang, Li-Yi; Bai, Yu
2015-01-01
In this paper, we discuss second-order consensus problems for multi-agent systems with dynamic agents and fixed topologies. A new second-order consensus protocol incorporating the predictive mechanism is proposed and a convergence analysis related to the real and imaginary parts of the eigenvalues of the Laplacian matrix of the corresponding network is provided. We establish a direct connection between the eigenvalues of the Laplacian matrix and parameters of the proposed second-order consensus model, and we compute a lower bound for the sum of parameters. It is proved under general conditions related to the second-order consensus protocol and the network topology that the asymptotic consensus is achieved and the convergence speed is increased via designing a state predictor. Finally, simulation results are provided to verify the effectiveness of the proposed protocol and the analytical claims.
Orthogonal canonical forms for second-order systems
NASA Technical Reports Server (NTRS)
Williams, Trevor; Laub, Alan J.
1992-01-01
It is shown that a linear second-order system with arbitrary damping cannot be reduced to Hessenberg-triangular form by means of orthogonal transformations. However, it is also shown that such an orthogonal reduction is always possible for the modal damping commonly assumed for models of flexible structures. It is shown that modally damped models can be orthogonally reduced to a new triangular second-order Schur form.
[Second-order retrospective revaluation in human contingency learning].
Numata, Keitaro; Shimazaki, Tsuneo
2009-04-01
We demonstrated second-order retrospective revaluation with three cues (T1, T2, and C) and an outcome, in human contingency learning. Experimental task, PC-controlled video game in which participants were required to observe about the relations between firing missiles and the tank destruction, consisted of three training phases and two rating phases. Groups C+ and C- consisted of same first two training phases, CT+ (cues C and T with an outcome) and T1T2+ followed by C+, or C- training for Groups C+, C-, respectively. In rating phases, it is clearly demonstrated that the judgment of predictive value for the outcome of the T2 were higher by C+ training (second-order unovershadowing) and lowered by C- training (second-order backward blocking). The results for Groups RC+ and RC-, in which the orders of the first two training phase for Groups C+ and C- were interchanged, also showed second-order unovershadowing and second-order backward blocking. These results, the robustness of second-order retrospective revaluation against the order of the first training phases, can be explained by the extended comparator hypothesis and probabilistic contrast model. However, these results cannot be explained by traditional associative learning models. PMID:19489431
The motion of ellipsoids in a second order fluid
NASA Astrophysics Data System (ADS)
Kim, S.
1985-09-01
The rigid body motion of an ellipsoid in a second order fluid (SOF) under the action of specified (time independent) external forces and torques have been obtained to first order in the Weissenberg number by inverting the resistance relations for the force an torque under specified rigid body motions. The reciprocal theorem of Lorentz was used to bypass the calculation of the O(W) velocity field. The results agree with known analytic solutions for SOF with the secondary to primary normal stress ratio of -1/2. The solution procedure was also tested by computing the torque on a translating prolate spheroid with aspect ratios ranging from slender bodies to near-spheres. One result is that for a SOF with zero secondary normal stress (Weissenberg fluid), previous asymptotic results for near-spheres were found to be accurate even at fairly large aspect ratios. New results of nondegenerate ellipsoids suggest that the orientation (as monitored by Euler angles) and trajectory of sedimenting, nonaxisymmetric particles such as ellipsoids provide useful information on the rheology of the suspending fluid.
Second order anisotropy contribution in perpendicular magnetic tunnel junctions
Timopheev, A. A.; Sousa, R.; Chshiev, M.; Nguyen, H. T.; Dieny, B.
2016-01-01
Hard-axis magnetoresistance loops were measured on perpendicular magnetic tunnel junction pillars of diameter ranging from 50 to 150 nm. By fitting these loops to an analytical model, the effective anisotropy fields in both free and reference layers were derived and their variations in temperature range between 340 K and 5 K were determined. It is found that a second-order anisotropy term of the form −K2cos4θ must be added to the conventional uniaxial –K1cos2θ term to explain the experimental data. This higher order contribution exists both in the free and reference layers. At T = 300 K, the estimated −K2/K1 ratios are 0.1 and 0.24 for the free and reference layers, respectively. The ratio is more than doubled at low temperatures changing the ground state of the reference layer from “easy-axis” to “easy-cone” regime. The easy-cone regime has clear signatures in the shape of the hard-axis magnetoresistance loops. The existence of this higher order anisotropy was also confirmed by ferromagnetic resonance experiments on FeCoB/MgO sheet films. It is of interfacial nature and is believed to be due to spatial fluctuations at the nanoscale of the first order anisotropy parameter at the FeCoB/MgO interface. PMID:27246631
Modeling Second-Order Chemical Reactions using Cellular Automata
NASA Astrophysics Data System (ADS)
Hunter, N. E.; Barton, C. C.; Seybold, P. G.; Rizki, M. M.
2012-12-01
Cellular automata (CA) are discrete, agent-based, dynamic, iterated, mathematical computational models used to describe complex physical, biological, and chemical systems. Unlike the more computationally demanding molecular dynamics and Monte Carlo approaches, which use "force fields" to model molecular interactions, CA models employ a set of local rules. The traditional approach for modeling chemical reactions is to solve a set of simultaneous differential rate equations to give deterministic outcomes. CA models yield statistical outcomes for a finite number of ingredients. The deterministic solutions appear as limiting cases for conditions such as a large number of ingredients or a finite number of ingredients and many trials. Here we present a 2-dimensional, probabilistic CA model of a second-order gas phase reaction A + B → C, using a MATLAB basis. Beginning with a random distribution of ingredients A and B, formation of C emerges as the system evolves. The reaction rate can be varied based on the probability of favorable collisions of the reagents A and B. The model permits visualization of the conversion of reagents to products, and allows one to plot concentration vs. time for A, B and C. We test hypothetical reaction conditions such as: limiting reagents, the effects of reaction probabilities, and reagent concentrations on the reaction kinetics. The deterministic solutions of the reactions emerge as statistical averages in the limit of the large number of cells in the array. Modeling results for dynamic processes in the atmosphere will be presented.
Second order anisotropy contribution in perpendicular magnetic tunnel junctions.
Timopheev, A A; Sousa, R; Chshiev, M; Nguyen, H T; Dieny, B
2016-01-01
Hard-axis magnetoresistance loops were measured on perpendicular magnetic tunnel junction pillars of diameter ranging from 50 to 150 nm. By fitting these loops to an analytical model, the effective anisotropy fields in both free and reference layers were derived and their variations in temperature range between 340 K and 5 K were determined. It is found that a second-order anisotropy term of the form -K2cos(4)θ must be added to the conventional uniaxial -K1cos(2)θ term to explain the experimental data. This higher order contribution exists both in the free and reference layers. At T = 300 K, the estimated -K2/K1 ratios are 0.1 and 0.24 for the free and reference layers, respectively. The ratio is more than doubled at low temperatures changing the ground state of the reference layer from "easy-axis" to "easy-cone" regime. The easy-cone regime has clear signatures in the shape of the hard-axis magnetoresistance loops. The existence of this higher order anisotropy was also confirmed by ferromagnetic resonance experiments on FeCoB/MgO sheet films. It is of interfacial nature and is believed to be due to spatial fluctuations at the nanoscale of the first order anisotropy parameter at the FeCoB/MgO interface. PMID:27246631
Second order anisotropy contribution in perpendicular magnetic tunnel junctions
NASA Astrophysics Data System (ADS)
Timopheev, A. A.; Sousa, R.; Chshiev, M.; Nguyen, H. T.; Dieny, B.
2016-06-01
Hard-axis magnetoresistance loops were measured on perpendicular magnetic tunnel junction pillars of diameter ranging from 50 to 150 nm. By fitting these loops to an analytical model, the effective anisotropy fields in both free and reference layers were derived and their variations in temperature range between 340 K and 5 K were determined. It is found that a second-order anisotropy term of the form ‑K2cos4θ must be added to the conventional uniaxial –K1cos2θ term to explain the experimental data. This higher order contribution exists both in the free and reference layers. At T = 300 K, the estimated ‑K2/K1 ratios are 0.1 and 0.24 for the free and reference layers, respectively. The ratio is more than doubled at low temperatures changing the ground state of the reference layer from “easy-axis” to “easy-cone” regime. The easy-cone regime has clear signatures in the shape of the hard-axis magnetoresistance loops. The existence of this higher order anisotropy was also confirmed by ferromagnetic resonance experiments on FeCoB/MgO sheet films. It is of interfacial nature and is believed to be due to spatial fluctuations at the nanoscale of the first order anisotropy parameter at the FeCoB/MgO interface.
Second order multidimensional sign-preserving remapping for ALE methods
Hill, Ryan N; Szmelter, J.
2010-12-15
A second-order conservative sign-preserving remapping scheme for Arbitrary Lagrangian-Eulerian (ALE) methods is developed utilising concepts of the Multidimensional Positive Definite Advection Transport Algorithm (MPDATA). The algorithm is inherently multidimensional, and so does not introduce splitting errors. The remapping is implemented in a two-dimensional, finite element ALE solver employing staggered quadrilateral meshes. The MPDATA remapping uses a finite volume discretization developed for volume coordinates. It is applied for the remapping of density and internal energy arranged as cell centered, and velocity as nodal, dependent variables. In the paper, the advection of scalar fields is examined first for test cases with prescribed mesh movement. A direct comparison of MPDATA with the performance of the van Leer MUSCL scheme indicates advantages of a multidimensional approach. Furthermore, distinctly different performance between basic MPDATA and the infinite gauge option is illustrated using benchmarks involving transport of a sign changing velocity field. Further development extends the application of MPDATA remapping to the full ALE solver with a staggered mesh arrangement for density, internal energy and momentum using volume coordinates. At present, two options of the algorithm - basic and infinite gauge - are implemented. To ensure a meaningful assessment, an identical Lagrangian solver and computational mesh update routines are used with either MPDATA or van Leer MUSCL remapping. The evaluation places particular focus on the abilities of both schemes to accurately model multidimensional problems. Theoretical considerations are supported with numerical examples. In addition to the prescribed mesh movement cases for advection of scalars, the demonstrations include two-dimensional Eulerian and ALE flow simulations on quadrilateral meshes with both fixed and variable timestep control. The key comparisons include the standard test cases of Sod and Noh
Gaussian Mixtures on Tensor Fields for Segmentation: Applications to Medical Imaging
de Luis-García, Rodrigo; Westin, Carl-Fredrik; Alberola-López, Carlos
2012-01-01
In this paper, we introduce a new approach for tensor field segmentation based on the definition of mixtures of Gaussians on tensors as a statistical model. Working over the well-known Geodesic Active Regions segmentation framework, this scheme presents several interesting advantages. First, it yields a more flexible model than the use of a single Gaussian distribution, which enables the method to better adapt to the complexity of the data. Second, it can work directly on tensor-valued images or, through a parallel scheme that processes independently the intensity and the local structure tensor, on scalar textured images. Two different applications have been considered to show the suitability of the proposed method for medical imaging segmentation. First, we address DT-MRI segmentation on a dataset of 32 volumes, showing a successful segmentation of the corpus callosum and favourable comparisons with related approaches in the literature. Second, the segmentation of bones from hand radiographs is studied, and a complete automatic-semiautomatic approach has been developed that makes use of anatomical prior knowledge to produce accurate segmentation results. PMID:20932717
Notes on Translational and Rotational Properties of Tensor Fields in Relativistic Quantum Mechanics
NASA Astrophysics Data System (ADS)
Dvoeglazov, V. V.
Recently, several discussions on the possible observability of 4-vector fields have been published in literature. Furthermore, several authors recently claimed existence of the helicity=0 fundamental field. We re-examine the theory of antisymmetric tensor fields and 4-vector potentials. We study the massless limits. In fact, a theoretical motivation for this venture is the old papers of Ogievetskiĭ and Polubarinov, Hayashi, and Kalb and Ramond. Ogievetskiĭ and Polubarinov proposed the concept of the notoph, whose helicity properties are complementary to those of the photon. We analyze the quantum field theory with taking into account mass dimensions of the notoph and the photon. It appears to be possible to describe both photon and notoph degrees of freedom on the basis of the modified Bargmann-Wigner formalism for the symmetric second-rank spinor. Next, we proceed to derive equations for the symmetric tensor of the second rank on the basis of the Bargmann-Wigner formalism in a straightforward way. The symmetric multispinor of the fourth rank is used. Due to serious problems with the interpretation of the results obtained on using the standard procedure we generalize it and obtain the spin-2 relativistic equations, which are consistent with the general relativity. Thus, in fact we deduced the gravitational field equations from relativistic quantum mechanics. The relations of this theory with the scalar-tensor theories of gravitation and f(R) are discussed. Particular attention has been paid to the correct definitions of the energy-momentum tensor and other Nöther currents in the electromagnetic theory, the relativistic theory of gravitation, the general relativity, and their generalizations. We estimate possible interactions, fermion-notoph, graviton-notoph, photon-notoph, and we conclude that they can probably be seen in experiments in the next few years.
Tensor hierarchy and generalized Cartan calculus in SL(3) × SL(2) exceptional field theory
NASA Astrophysics Data System (ADS)
Hohm, Olaf; Wang, Yi-Nan
2015-04-01
We construct exceptional field theory for the duality group SL(3) × SL(2). The theory is defined on a space with 8 `external' coordinates and 6 `internal' coordinates in the (3, 2) fundamental representation, leading to a 14-dimensional generalized spacetime. The bosonic theory is uniquely determined by gauge invariance under generalized external and internal diffeomorphisms. The latter invariance can be made manifest by introducing higher form gauge fields and a so-called tensor hierarchy, which we systematically develop to much higher degree than in previous studies. To this end we introduce a novel Cartan-like tensor calculus based on a covariant nil-potent differential, generalizing the exterior derivative of conventional differential geometry. The theory encodes the full D = 11 or type IIB supergravity, respectively.
Visualizing Tensor Normal Distributions at Multiple Levels of Detail.
Abbasloo, Amin; Wiens, Vitalis; Hermann, Max; Schultz, Thomas
2016-01-01
Despite the widely recognized importance of symmetric second order tensor fields in medicine and engineering, the visualization of data uncertainty in tensor fields is still in its infancy. A recently proposed tensorial normal distribution, involving a fourth order covariance tensor, provides a mathematical description of how different aspects of the tensor field, such as trace, anisotropy, or orientation, vary and covary at each point. However, this wealth of information is far too rich for a human analyst to take in at a single glance, and no suitable visualization tools are available. We propose a novel approach that facilitates visual analysis of tensor covariance at multiple levels of detail. We start with a visual abstraction that uses slice views and direct volume rendering to indicate large-scale changes in the covariance structure, and locations with high overall variance. We then provide tools for interactive exploration, making it possible to drill down into different types of variability, such as in shape or orientation. Finally, we allow the analyst to focus on specific locations of the field, and provide tensor glyph animations and overlays that intuitively depict confidence intervals at those points. Our system is demonstrated by investigating the effects of measurement noise on diffusion tensor MRI, and by analyzing two ensembles of stress tensor fields from solid mechanics. PMID:26529741
Tensor of the nonlinear polarizability of anisotropic medium and ``local'' field method
NASA Astrophysics Data System (ADS)
Lavric, V. V.; Ovander, L. N.; Shunyakov, V. T.
1983-08-01
The nonlinear polarizability tensor (NPT) for a molecular crystal of arbitrary symmetry has been obtained within the framework of polariton theory. Use of the Göppert-Mayer unitary transformation for the Hamiltonian of the crystal plus quantized electromagnetic field system made it possible to represent finally the result for the NPT in a compact form and to compare with results of semiphenomenological calculation of the NPT and to go out of the framework of the Gaitler-London approximation.
Second-order quasinormal mode of the Schwarzschild black hole
Nakano, Hiroyuki; Ioka, Kunihito
2007-10-15
We formulate and calculate the second-order quasinormal modes (QNMs) of a Schwarzschild black hole (BH). Gravitational waves (GW) from a distorted BH, the so-called ringdowns, are well understood as QNMs in general relativity. Since QNMs from binary BH mergers will be detected with a high signal-to-noise ratio by GW detectors, it is also possible to detect the second perturbative order of QNMs, generated by nonlinear gravitational interaction near the BH. In the BH perturbation approach, we derive the master Zerilli equation for the metric perturbation to second order and explicitly regularize it at the horizon and spatial infinity. We numerically solve the second-order Zerilli equation by implementing the modified Leaver continued fraction method. The second-order QNM frequencies are found to be twice the first-order ones, and the GW amplitude is up to {approx}10% that of the first order for the binary BH mergers. Since the second-order QNMs always exist, we can use their detections (i) to test the nonlinearity of general relativity, in particular, the no-hair theorem, (ii) to remove fake events in the data analysis of QNM GWs, and (iii) to measure the distance to the BH.
Second-Order Olfactory-Mediated Fear-Potentiated Startle
Paschall, Gayla Y.; Davis, Michael
2002-01-01
Recently, we reported that discrete (4-sec) olfactory cues paired with footshock serve as effective conditioned stimuli (CSs) for potentiating the acoustic startle response in rats using the fear-potentiated startle paradigm. Because odors are such salient cues for the rat, and because of the robust olfactory conditioning observed previously, the current studies investigated second-order fear conditioning using olfactory and visual cues. In Experiments 1 and 2, we used a small number of first-order and second-order training trials on separate days to investigate second-order fear-potentiated startle. Significant potentiated startle was observed in animals receiving Paired/Paired training in both studies, but surprisingly, control animals in the Unpaired/Paired group (Exp. 1) also showed significant potentiated startle to a light S2 at testing. These findings are addressed in the Discussion. Overall, the results of both experiments suggest that olfactory cues serve as efficient S1 and S2 stimuli in second-order fear-potentiated startle paradigms when only a small number of first and second-order training trials are presented. PMID:12464699
Weak value amplification via second-order correlated technique
NASA Astrophysics Data System (ADS)
Ting, Cui; Jing-Zheng, Huang; Xiang, Liu; Gui-Hua, Zeng
2016-02-01
We propose a new framework combining weak measurement and second-order correlated technique. The theoretical analysis shows that weak value amplification (WVA) experiment can also be implemented by a second-order correlated system. We then build two-dimensional second-order correlated function patterns for achieving higher amplification factor and discuss the signal-to-noise ratio influence. Several advantages can be obtained by our proposal. For instance, detectors with high resolution are not necessary. Moreover, detectors with low saturation intensity are available in WVA setup. Finally, type-one technical noise can be effectively suppressed. Project supported by the Union Research Centre of Advanced Spaceflight Technology (Grant No. USCAST2013-05), the National Natural Science Foundation of China (Grant Nos. 61170228, 61332019, and 61471239), and the High-Tech Research and Development Program of China (Grant No. 2013AA122901).
Some restrictions on the existence of second order limit language
NASA Astrophysics Data System (ADS)
Ahmad, Muhammad Azrin; Sarmin, Nor Haniza; Yusof, Yuhani; Fong, Wan Heng
2015-10-01
The cut and paste phenomenon on DNA molecules with the presence of restriction enzyme and appropriate ligase has led to the formalism of mathematical modelling of splicing system. A type of splicing system named Yusof-Goode splicing system is used to present the transparent behaviour of the DNA splicing process. The limit language that is defined as the leftover molecules after the system reaches its equilibrium point has been extended to a second order limit language. The non-existence of the second order limit language biologically has lead to this study by using mathematical approach. In this paper, the factors that restrict the formation of the second order limit language are discussed and are presented as lemmas and theorem using Y-G approach. In addition, the discussion focuses on Yusof- Goode splicing system with at most two initial strings and two rules with one cutting site and palindromic crossing site and recognition sites.
Textural segmentation, second-order statistics, and textural elements.
Beck, J
1983-01-01
Beck (1972, 1973) hypothesized that textural segmentation occurs strongly on the basis of simple properties such as brightness, color, size, and the slopes of contours and lines of the elemental descriptors of a texture or textural elements. The experiment reported supports the hypothesis that specific stimulus features, rather than second-order statistics, account for textural segmentation. The results agree with Julesz (1981a,b) who has reported evidence disproving his original conjecture of the importance of second-order statistics. Julesz (1981a,b) now hypothesizes textural segmentation to be a function of local features which he called textons. Textons are features that give textural segmentation when textures have identical second-order statistics. The two hypotheses are to date in complete agreement on the stimulus features producing textural segmentation, and the experiment reported is consistent with both. PMID:6626590
Deflection of light to second order in conformal Weyl gravity
NASA Astrophysics Data System (ADS)
Sultana, Joseph
2013-04-01
We reexamine the deflection of light in conformal Weyl gravity obtained in Sultana and Kazanas (2010), by extending the calculation based on the procedure by Rindler and Ishak, for the bending angle by a centrally concentrated spherically symmetric matter distribution, to second order in M/R, where M is the mass of the source and R is the impact parameter. It has recently been reported in Bhattacharya et al. (JCAP 09 (2010) 004; JCAP 02 (2011) 028), that when this calculation is done to second order, the term γr in the Mannheim-Kazanas metric, yields again the paradoxical contribution γR (where the bending angle is proportional to the impact parameter) obtained by standard formalisms appropriate to asymptotically flat spacetimes. We show that no such contribution is obtained for a second order calculation and the effects of the term γr in the metric are again insignificant as reported in our earlier work.
Second-order wave effects on TLP tendon tension responses
Xue, H.; Mercier, R.S.
1996-12-31
This paper presents a general procedure for analyzing the second-order wave effects on the tendon tension responses of a TLP. The approach solves both first- and second-order equation of motions for a TLP system in frequency domain. Viscous effects are included in the form of statistically linearized damping coefficients. An efficient algorithm has been devised for reducing the burden of second-order wave diffraction analysis, which selects the interacting frequency pairs according to springing frequency of interest to minimize the cost of computing quadratic transfer functions (QTFs) and allow accurate interpolation of QTFs. Moment statistics of the tension process are computed through an eigenvalue analysis. The developed method is applied to analyze the tendon tension responses of a TLP design in water depth of 3,000 ft.
Second-order subsonic airfoil theory including edge effects
NASA Technical Reports Server (NTRS)
Van Dyke, Milton D
1956-01-01
Several recent advances in plane subsonic flow theory are combined into a unified second-order theory for airfoil sections of arbitrary shape. The solution is reached in three steps: the incompressible result is found by integration, it is converted into the corresponding subsonic compressible result by means of the second-order compressibility rule, and it is rendered uniformly valid near stagnation points by further rules. Solutions for a number of airfoils are given and are compared with the results of other theories and of experiment. A straight-forward computing scheme is outlined for calculating the surface velocities and pressures on any airfoil at any angle of attack
Solution of second order supersymmetrical intertwining relations in Minkowski plane
NASA Astrophysics Data System (ADS)
Ioffe, M. V.; Kolevatova, E. V.; Nishnianidze, D. N.
2016-08-01
Supersymmetrical (SUSY) intertwining relations are generalized to the case of quantum Hamiltonians in Minkowski space. For intertwining operators (supercharges) of second order in derivatives, the intertwined Hamiltonians correspond to completely integrable systems with the symmetry operators of fourth order in momenta. In terms of components, the intertwining relations correspond to the system of nonlinear differential equations which are solvable with the simplest—constant—ansatzes for the "metric" matrix in second order part of the supercharges. The corresponding potentials are built explicitly both for diagonalizable and nondiagonalizable form of "metric" matrices, and their properties are discussed.
Controlling flexible structures with second order actuator dynamics
NASA Technical Reports Server (NTRS)
Inman, Daniel J.; Umland, Jeffrey W.; Bellos, John
1989-01-01
The control of flexible structures for those systems with actuators that are modeled by second order dynamics is examined. Two modeling approaches are investigated. First a stability and performance analysis is performed using a low order finite dimensional model of the structure. Secondly, a continuum model of the flexible structure to be controlled, coupled with lumped parameter second order dynamic models of the actuators performing the control is used. This model is appropriate in the modeling of the control of a flexible panel by proof-mass actuators as well as other beam, plate and shell like structural numbers. The model is verified with experimental measurements.
Second-order accurate nonoscillatory schemes for scalar conservation laws
NASA Technical Reports Server (NTRS)
Huynh, Hung T.
1989-01-01
Explicit finite difference schemes for the computation of weak solutions of nonlinear scalar conservation laws is presented and analyzed. These schemes are uniformly second-order accurate and nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time.
Second-order sliding mode control with experimental application.
Eker, Ilyas
2010-07-01
In this article, a second-order sliding mode control (2-SMC) is proposed for second-order uncertain plants using equivalent control approach to improve the performance of control systems. A Proportional + Integral + Derivative (PID) sliding surface is used for the sliding mode. The sliding mode control law is derived using direct Lyapunov stability approach and asymptotic stability is proved theoretically. The performance of the closed-loop system is analysed through an experimental application to an electromechanical plant to show the feasibility and effectiveness of the proposed second-order sliding mode control and factors involved in the design. The second-order plant parameters are experimentally determined using input-output measured data. The results of the experimental application are presented to make a quantitative comparison with the traditional (first-order) sliding mode control (SMC) and PID control. It is demonstrated that the proposed 2-SMC system improves the performance of the closed-loop system with better tracking specifications in the case of external disturbances, better behavior of the output and faster convergence of the sliding surface while maintaining the stability. PMID:20413118
Second-order nonlinear optical metamaterials: ABC-type nanolaminates
Alloatti, L. Kieninger, C.; Lauermann, M.; Köhnle, K.; Froelich, A.; Wegener, M.; Frenzel, T.; Freude, W.; Leuthold, J.; Koos, C.
2015-09-21
We demonstrate a concept for second-order nonlinear metamaterials that can be obtained from non-metallic centrosymmetric constituents with inherently low optical absorption. The concept is based on iterative atomic-layer deposition of three different materials, A = Al{sub 2}O{sub 3}, B = TiO{sub 2}, and C = HfO{sub 2}. The centrosymmetry of the resulting ABC stack is broken since the ABC and the inverted CBA sequences are not equivalent—a necessary condition for non-zero second-order nonlinearity. In our experiments, we find that the bulk second-order nonlinear susceptibility depends on the density of interfaces, leading to a nonlinear susceptibility of 0.26 pm/V at a wavelength of 800 nm. ABC-type nanolaminates can be deposited on virtually any substrate and offer a promising route towards engineering of second-order optical nonlinearities at both infrared and visible wavelengths.
Children's Understanding of Second-Order Mental States
ERIC Educational Resources Information Center
Miller, Scott A.
2009-01-01
The most popular topic in theory-of-mind research has been first-order false belief: the realization that it is possible to hold false beliefs about events in the world. A more advanced development is second-order false belief: the realization that it is possible to hold a false belief about someone else's belief. This article reviews research…
Second-order variational equations for N-body simulations
NASA Astrophysics Data System (ADS)
Rein, Hanno; Tamayo, Daniel
2016-07-01
First-order variational equations are widely used in N-body simulations to study how nearby trajectories diverge from one another. These allow for efficient and reliable determinations of chaos indicators such as the Maximal Lyapunov characteristic Exponent (MLE) and the Mean Exponential Growth factor of Nearby Orbits (MEGNO). In this paper we lay out the theoretical framework to extend the idea of variational equations to higher order. We explicitly derive the differential equations that govern the evolution of second-order variations in the N-body problem. Going to second order opens the door to new applications, including optimization algorithms that require the first and second derivatives of the solution, like the classical Newton's method. Typically, these methods have faster convergence rates than derivative-free methods. Derivatives are also required for Riemann manifold Langevin and Hamiltonian Monte Carlo methods which provide significantly shorter correlation times than standard methods. Such improved optimization methods can be applied to anything from radial-velocity/transit-timing-variation fitting to spacecraft trajectory optimization to asteroid deflection. We provide an implementation of first- and second-order variational equations for the publicly available REBOUND integrator package. Our implementation allows the simultaneous integration of any number of first- and second-order variational equations with the high-accuracy IAS15 integrator. We also provide routines to generate consistent and accurate initial conditions without the need for finite differencing.
Second-Order Conditioning during a Compound Extinction Treatment
ERIC Educational Resources Information Center
Pineno, Oskar; Zilski, Jessica M.; Schachtman, Todd R.
2007-01-01
Two conditioned taste aversion experiments with rats were conducted to establish if a target taste that had received a prior pairing with illness could be subject to second-order conditioning during extinction treatment in compound with a flavor that also received prior conditioning. In these experiments, the occurrence of second-order…
Green's function of the second order differential operator with involution
NASA Astrophysics Data System (ADS)
Ashyralyev, Allaberen; Sarsenbi, Abdisalam A.
2016-08-01
In the present paper, the Green's function of the second order differential operator L defined by formula L u =α u″ (x ) -u″ (-x ) =λ u (x ) ,-1
PREDICTIONS OF HIGHWAY EMISSIONS BY A SECOND ORDER CLOSURE MODEL
The dispersion of sulfur hexafluoride tracer and sulfate from automobile emissions in the immediate vicinity of a highway were estimated for conditions similar to those existing during the General Motors sulfate dispersion experiment conducted at a GM test track. A second-order c...
Second-Order Conditioning of Human Causal Learning
ERIC Educational Resources Information Center
Jara, Elvia; Vila, Javier; Maldonado, Antonio
2006-01-01
This article provides the first demonstration of a reliable second-order conditioning (SOC) effect in human causal learning tasks. It demonstrates the human ability to infer relationships between a cause and an effect that were never paired together during training. Experiments 1a and 1b showed a clear and reliable SOC effect, while Experiments 2a…
Modeling Ability Differentiation in the Second-Order Factor Model
ERIC Educational Resources Information Center
Molenaar, Dylan; Dolan, Conor V.; van der Maas, Han L. J.
2011-01-01
In this article we present factor models to test for ability differentiation. Ability differentiation predicts that the size of IQ subtest correlations decreases as a function of the general intelligence factor. In the Schmid-Leiman decomposition of the second-order factor model, we model differentiation by introducing heteroscedastic residuals,…
Pecuniary Effects, Second-Order Conditions, and the LRAC Curve.
ERIC Educational Resources Information Center
Comolli, Paul M.
2000-01-01
Explores the importance of second-order conditions in the cost-minimization problem confronting the monopsonistic employer of factor inputs. Describes an alternative approach to the presence of pecuniary effects that does not depend on the assumption that firms are monopsonistic in factor markets. (CMK)
A New Factorisation of a General Second Order Differential Equation
ERIC Educational Resources Information Center
Clegg, Janet
2006-01-01
A factorisation of a general second order ordinary differential equation is introduced from which the full solution to the equation can be obtained by performing two integrations. The method is compared with traditional methods for solving these type of equations. It is shown how the Green's function can be derived directly from the factorisation…
Forward and Backward Second-Order Pavlovian Conditioning in Honeybees
ERIC Educational Resources Information Center
Hussaini, Syed Abid; Komischke, Bernhard; Menzel, Randolf; Lachnit, Harald
2007-01-01
Second-order conditioning (SOC) is the association of a neutral stimulus with another stimulus that had previously been combined with an unconditioned stimulus (US). We used classical conditioning of the proboscis extension response (PER) in honeybees ("Apis mellifera") with odors (CS) and sugar (US). Previous SOC experiments in bees were…
Generalized Second-Order Partial Derivatives of 1/r
ERIC Educational Resources Information Center
Hnizdo, V.
2011-01-01
The generalized second-order partial derivatives of 1/r, where r is the radial distance in three dimensions (3D), are obtained using a result of the potential theory of classical analysis. Some non-spherical-regularization alternatives to the standard spherical-regularization expression for the derivatives are derived. The utility of a…
Solving Second-Order Differential Equations with Variable Coefficients
ERIC Educational Resources Information Center
Wilmer, A., III; Costa, G. B.
2008-01-01
A method is developed in which an analytical solution is obtained for certain classes of second-order differential equations with variable coefficients. By the use of transformations and by repeated iterated integration, a desired solution is obtained. This alternative method represents a different way to acquire a solution from classic power…
Estimation and Application of Spatially Variable Noise Fields in Diffusion Tensor Imaging
Landman, Bennett A.; Bazin, Pierre-Louis; Prince, Jerry L.
2009-01-01
Optimal interpretation of magnetic resonance image content often requires an estimate of the underlying image noise, which is typically realized as a spatially invariant estimate of the noise distribution. This is not an ideal practice in diffusion tensor imaging because the noise distribution is usually spatially varying due to the use of fast imaging and noise suppression techniques. A new estimation approach for spatially varying noise fields is proposed in this article. The approach is based on a noise invariance property in scenarios in which more than one image, each with potentially different signal levels, are acquired on each slice, as in diffusion weighted MRI. This technique leads to improved noise field estimates in simulations, phantom experiments and in vivo studies when compared to traditional noise field estimators that use regional variability or background intensity histograms. The proposed method reduces the noise field estimation error by a factor of 100 in simulations, shows a strong linear correlation (R2 = 0.99) between theoretical and estimated noise changes in phantoms, and demonstrates consistent (<5% variability) noise field estimates in vivo. The advantages of spatially varying noise field estimation are demonstrated for power analysis, outlier detection, and tensor estimation. PMID:19250784
Consistency of Equations in the Second-Order Gauge-Invariant Cosmological Perturbation Theory
NASA Astrophysics Data System (ADS)
Nakamura, K.
2009-06-01
Along the general framework of the gauge-invariant perturbation theory developed in the papers [K.~Nakamura, Prog.~Theor.~Phys. 110 (2003), 723; Prog.~Theor.~Phys. 113 (2005), 481], we rederive the second-order Einstein equation on four-dimensional homogeneous isotropic background universe in a gauge-invariant manner without ignoring any mode of perturbations. We consider the perturbations both in the universe dominated by the single perfect fluid and in that dominated by the single scalar field. We also confirmed the consistency of all the equations of the second-order Einstein equation and the equations of motion for matter fields, which are derived in the paper [K.~Nakamura, arXiv:0804.3840]. This confirmation implies that all the derived equations of the second order are self-consistent and these equations are correct in this sense.
First estimates of the second-order ionospheric effect on radio occultation observations
NASA Astrophysics Data System (ADS)
Vergados, Panagiotis; Pagiatakis, Spiros D.
2010-07-01
This study examines the impact of the second-order ionospheric effect on radio occultation (RO) data products. We propose a new linear combination between dual frequency GPS observables, which retrieves slant total electron content free from the second-order ionospheric effect. Our STEC values differ from those obtained by independent techniques by a maximum of 3 total electron content units (TECU), depending on the geographic location and geomagnetic activity. Additionally, we suggest an alternative method of computing the second-order ionospheric delay in RO experiments, which does not require the use of geomagnetic and ionospheric models. First estimates show that the second-order ionospheric delay for the RO experiments falls within the range [-10, -8] mm, which is of the same order of magnitude with second-order ionospheric delay estimates from ground-based experiments. Finally, as a by-product of our model, we retrieve weighted mean geomagnetic field values, which we compare with theoretical estimates computed by the International Geomagnetic Reference Field-10 (IGRF-10) model. Our estimations agree with the IGRF-10 model between 0.23% and 7.0%.
Large tensor-to-scalar ratio in small-field inflation.
Kobayashi, Takeshi; Takahashi, Tomo
2013-06-01
We show that density perturbations seeded by the inflaton can be suppressed when having additional light degrees of freedom contributing to the production of perturbations. The inflaton fluctuations affect the light field dynamics by modulating the length of the inflationary period and, hence, produce additional density perturbations in the postinflationary era. Such perturbations can cancel those generated during inflation as both originate from the same inflaton fluctuations. This allows production of large gravitational waves from small-field inflation, which is normally forbidden by the Lyth bound on the inflaton field excursion. We also find that the field bound is taken over by the light scalar when the inflaton-induced perturbations are suppressed and, thus, present a generalized form of the Lyth bound that applies to the total field space. The novel mechanism allows violation of the usual consistency relation r≤-8n(T) for the tensor spectral index. PMID:25167480
Travelling-wave Green tensor and near-field Rayleigh-wave sensitivity
NASA Astrophysics Data System (ADS)
Liu, Kui; Zhou, Ying
2016-04-01
Travelling-wave Green tensor has been widely used in calculations of synthetic seismograms and finite-frequency sensitivities of surface waves. The classic travelling-wave decomposition is based on a far-field approximation and may not be valid when applied to construct sensitivity kernels in regions close to the receiver. In this paper, we calculate synthetic seismograms and finite-frequency sensitivity kernels of Rayleigh waves based on travelling-wave representation of Green tensor that fully accounts for near-field effects. We show that far-field approximation is adequate for synthetic seismograms when the source-receiver epicentral distance is greater than the dominant wavelength. Errors in Rayleigh-wave sensitivity kernels introduced by far-field approximation are in general negligible for single-station measurements except for in a small region around the station, and the errors are more significant in sensitivity kernels for interstation measurements. In addition, interstation measurements are strongly sensitive to structures outside the region between the two stations, even for two stations along the same great circle path from the seismic source.
NASA Technical Reports Server (NTRS)
Shih, Tsan-Hsing; Lumley, John L.
1991-01-01
Recently, several second order closure models have been proposed for closing the second moment equations, in which the velocity-pressure gradient (and scalar-pressure gradient) tensor and the dissipation rate tensor are two of the most important terms. In the literature, these correlation tensors are usually decomposed into a so called rapid term and a return-to-isotropy term. Models of these terms have been used in global flow calculations together with other modeled terms. However, their individual behavior in different flows have not been fully examined because they are un-measurable in the laboratory. Recently, the development of direct numerical simulation (DNS) of turbulence has given us the opportunity to do this kind of study. With the direct numerical simulation, we may use the solution to exactly calculate the values of these correlation terms and then directly compare them with the values from their modeled formulations (models). Here, we make direct comparisons of five representative rapid models and eight return-to-isotropy models using the DNS data of forty five homogeneous flows which were done by Rogers et al. (1986) and Lee et al. (1985). The purpose of these direct comparisons is to explore the performance of these models in different flows and identify the ones which give the best performance. The modeling procedure, model constraints, and the various evaluated models are described. The detailed results of the direct comparisons are discussed, and a few concluding remarks on turbulence models are given.
A second order operator splitting method for Allen-Cahn type equations with nonlinear source terms
NASA Astrophysics Data System (ADS)
Lee, Hyun Geun; Lee, June-Yub
2015-08-01
Allen-Cahn (AC) type equations with nonlinear source terms have been applied to a wide range of problems, for example, the vector-valued AC equation for phase separation and the phase-field equation for dendritic crystal growth. In contrast to the well developed first and second order methods for the AC equation, not many second order methods are suggested for the AC type equations with nonlinear source terms due to the difficulties in dealing with the nonlinear source term numerically. In this paper, we propose a simple and stable second order operator splitting method. A core idea of the method is to decompose the original equation into three subequations with the free-energy evolution term, the heat evolution term, and a nonlinear source term, respectively. It is important to combine these three subequations in proper order to achieve the second order accuracy and stability. We propose a method with a half-time free-energy evolution solver, a half-time heat evolution solver, a full-time midpoint solver for the nonlinear source term, and a half-time heat evolution solver followed by a final half-time free-energy evolution solver. We numerically demonstrate the second order accuracy of the new numerical method through the simulations of the phase separation and the dendritic crystal growth.
An application of integral inequality to second order nonlinear oscillation
NASA Astrophysics Data System (ADS)
Kwong, Man Kam; Wong, James S. W.
A simple result concerning integral inequalities enables us to give an alternative proof of Waltman's theorem: lim t → ∞ ∝ t0a( s) ds = ∞ implies oscillation of the second order nonlinear equation y″( t) + a( t) f( y( t)) = 0; to prove an analog of Wintner's theorem that relates the nonoscillation of the second order nonlinear equations to the existence of solutions of some integral equations, assuming that lim t → ∞ ∝ t0a( s) ds exists; and to give an alternative proof and to extend a result of Butler. An often used condition on the coefficient a( t) is given a more familiar equivalent form and an oscillation criterion involving this condition is established.
Robust eigensystem assignment for second-order estimators
NASA Technical Reports Server (NTRS)
Juang, Jer-Nan; Maghami, Peiman G.
1990-01-01
An approach for the robust eigensystem assignment of flexible structures using full state or output feedback is developed. Using the second-order dynamic equations, the approach can assign the eigenvalues of the system via velocity and displacement feedbacks, or acceleration and velocity feedbacks. The eigenvalues and eigenvectors of the system are assigned, via the second-order eigenvalue problem for the structural system, in two steps. First, an orthonormal basis spanning the attainable closed-loop eigenvector space corresponding to each desired closed-loop eigenvalue is generated using the Singular Value or QR decompositions. Second, a sequential procedure is used to choose a set of closed-loop eigenvectors that are as close as possible to the column space of a well-conditioned target matrix. Among the possible choices of the target matrix, the closest unitary matrix to the open-loop eigenvector matrix appears to be a suitable choice. A numerical example is given to illustrate the proposed algorithm.
On the second order spatiochromatic structure of natural images.
Provenzi, Edoardo; Delon, Julie; Gousseau, Yann; Mazin, Baptiste
2016-03-01
We provide a theoretical analysis of some empirical facts about the second order spatiochromatic structure of natural images in color. In particular, we show that two simple assumptions on the covariance matrices of color images yield eigenvectors made by the Kronecker product of Fourier features times the triad given by luminance plus color opponent channels. The first of these assumptions is second order stationarity while the second one is commutativity between color correlation matrices. The validity of these assumptions and the predicted shape of the PCA components of color images are experimentally observed on two large image databases. As a by-product of this experimental study, we also provide novel data to support an exponential decay law of the spatiochromatic covariance between pairs of pixels as a function of their spatial distance. PMID:26024561
A Second-Order Achromat Design Based on FODO Cell
Sun, Yipeng; /SLAC
2011-08-19
Two dipole doglegs are widely used to translate the beam axis horizontally or vertically. Quadrupoles are placed between the two consecutive dipoles to match first order dispersion and provide betatron focusing. Similarly a four dipole chicane is usually employed to form a bypass region, where the beam axis is transversely shifted first, then translated back to the original axis. In order to generate an isochronous section, quadrupoles are again needed to tune the first order transfer matrix element R{sub 56} equaling zero. Usually sextupoles are needed to correct second order dispersion in the bending plane, for both the dogleg optics and the chicane (with quad) optics. In this paper, an alternative optics design is introduced, which is based on a simple FODO cell and does not need sextupoles assistance to form a second-order achromat. It may provide a similar function of either a dogleg or a bypass, by using 2 or 4 of such combined supercells.
Second order upwind Lagrangian particle method for Euler equations
Samulyak, Roman; Chen, Hsin -Chiang; Yu, Kwangmin
2016-06-01
A new second order upwind Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface / multiphase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) an upwind second-order particle-based algorithm with limiter, providing accuracy and longmore » term stability, and (c) accurate resolution of states at free interfaces. In conclusion, numerical verification tests demonstrating the convergence order for fixed domain and free surface problems are presented.« less
Analysis of the blazing effect in second-order gratings
Matsumoto, Masauki )
1992-10-01
The validity of the blazing effect for improving the output-coupling efficiency of second-order gratings for use in grating-coupled surface-emitting lasers is examined. The Floquet-Bloch expansion method is used for the analysis of finite-length gratings with asymmetric tooth shapes operated in resonance condition. It is shown that for saw-tooth gratings the blazing effect is almost lost around the second-order Bragg wavelength because the reflected guide mode generated in the distributed Bragg reflector radiates preferentially into the substrate. By using a parallelogram grating with equal tooth and groove lengths, however, a high efficiency of radiation into the air is attainable even at the Bragg wavelength while the reflectivity is reduced. 32 refs.
Experimental Measurement of the Second-Order Coherence of Supercontinuum.
Närhi, Mikko; Turunen, Jari; Friberg, Ari T; Genty, Goëry
2016-06-17
We measure experimentally the second-order coherence properties of supercontinuum generated in a photonic crystal fiber. Our approach is based on measuring separately the quasicoherent and quasistationary contributions to the cross-spectral density and mutual coherence functions using a combination of interferometric and nonlinear gating techniques. This allows us to introduce two-dimensional coherence spectrograms which provide a direct characterization and convenient visualization of the spectrotemporal coherence properties. The measured second-order coherence functions are in very good agreement with numerical simulations based on the generalized nonlinear Schrödinger equation. Our results pave the way towards the full experimental characterization of supercontinuum coherence properties. More generally, they provide a generic approach for the complete experimental measurement of the coherence of broadband sources. PMID:27367389
Experimental Measurement of the Second-Order Coherence of Supercontinuum
NASA Astrophysics Data System (ADS)
Närhi, Mikko; Turunen, Jari; Friberg, Ari T.; Genty, Goëry
2016-06-01
We measure experimentally the second-order coherence properties of supercontinuum generated in a photonic crystal fiber. Our approach is based on measuring separately the quasicoherent and quasistationary contributions to the cross-spectral density and mutual coherence functions using a combination of interferometric and nonlinear gating techniques. This allows us to introduce two-dimensional coherence spectrograms which provide a direct characterization and convenient visualization of the spectrotemporal coherence properties. The measured second-order coherence functions are in very good agreement with numerical simulations based on the generalized nonlinear Schrödinger equation. Our results pave the way towards the full experimental characterization of supercontinuum coherence properties. More generally, they provide a generic approach for the complete experimental measurement of the coherence of broadband sources.
Second-order infinitesimal bendings of surfaces of revolution with flattening at the poles
NASA Astrophysics Data System (ADS)
Sabitov, I. Kh
2014-12-01
We study infinitesimal bendings of surfaces of revolution with flattening at the poles. We begin by considering the minimal possible smoothness class C^1 both for surfaces and for deformation fields. Conditions are formulated for a given harmonic of a first-order infinitesimal bending to be extendable into a second order infinitesimal bending. We finish by stating a criterion for nonrigidity of second order for closed surfaces of revolution in the analytic class. We also give the first concrete example of such a nonrigid surface. Bibliography: 15 entries.
Second-order infinitesimal bendings of surfaces of revolution with flattening at the poles
Sabitov, I Kh
2014-12-31
We study infinitesimal bendings of surfaces of revolution with flattening at the poles. We begin by considering the minimal possible smoothness class C{sup 1} both for surfaces and for deformation fields. Conditions are formulated for a given harmonic of a first-order infinitesimal bending to be extendable into a second order infinitesimal bending. We finish by stating a criterion for nonrigidity of second order for closed surfaces of revolution in the analytic class. We also give the first concrete example of such a nonrigid surface. Bibliography: 15 entries.
Extensions and applications of a second-order landsurface parameterization
NASA Technical Reports Server (NTRS)
Andreou, S. A.; Eagleson, P. S.
1983-01-01
Extensions and applications of a second order land surface parameterization, proposed by Andreou and Eagleson are developed. Procedures for evaluating the near surface storage depth used in one cell land surface parameterizations are suggested and tested by using the model. Sensitivity analysis to the key soil parameters is performed. A case study involving comparison with an "exact" numerical model and another simplified parameterization, under very dry climatic conditions and for two different soil types, is also incorporated.
Gravitational waves from global second order phase transitions
Jr, John T. Giblin; Price, Larry R.; Siemens, Xavier; Vlcek, Brian E-mail: larryp@caltech.edu E-mail: bvlcek@uwm.edu
2012-11-01
Global second-order phase transitions are expected to produce scale-invariant gravitational wave spectra. In this manuscript we explore the dynamics of a symmetry-breaking phase transition using lattice simulations. We explicitly calculate the stochastic gravitational wave background produced during the transition and subsequent self-ordering phase. We comment on this signal as it compares to the scale-invariant spectrum produced during inflation.
Second order filter response with series coupled silica microresonators
NASA Technical Reports Server (NTRS)
Savchenkov, A.; Iitchenko, V. S.; Handley, T.; Maleki, L.
2002-01-01
We have demonstrated an approach for fabricating a photonic filter with second order response function. The filter consists of two germania-doped silica microtoroidal or microspherical resonators cascaded in series. We use UV irradiation to tune the mode of one microcavity to bring it close to the mode of the second microcavity. This approach produces a filter function with much sharper rolloff than can be obtained with the individual microresonators.
athena: Tree code for second-order correlation functions
NASA Astrophysics Data System (ADS)
Kilbinger, Martin; Bonnett, Christopher; Coupon, Jean
2014-02-01
athena is a 2d-tree code that estimates second-order correlation functions from input galaxy catalogues. These include shear-shear correlations (cosmic shear), position-shear (galaxy-galaxy lensing) and position-position (spatial angular correlation). Written in C, it includes a power-spectrum estimator implemented in Python; this script also calculates the aperture-mass dispersion. A test data set is available.
Second order parametric processes in nonlinear silica microspheres.
Xu, Yong; Han, Ming; Wang, Anbo; Liu, Zhiwen; Heflin, James R
2008-04-25
We analyze second order parametric processes in a silica microsphere coated with radially aligned nonlinear optical molecules. In a high-Q nonlinear microsphere, we discover that it is possible to achieve ultralow threshold parametric oscillation that obeys the rule of angular momentum conservation. Based on symmetry considerations, one can also implement parametric processes that naturally generate quantum entangled photon pairs. Practical issues regarding implementation of the nonlinear microsphere are also discussed. PMID:18518201
Theoretical study of second-order hyperpolarizability for nitrogen radical cation
NASA Astrophysics Data System (ADS)
Tarazkar, Maryam; Romanov, Dmitri A.; Levis, Robert J.
2015-05-01
We report calculations of the static and dynamic hyperpolarizabilities of the nitrogen radical cation in doublet state. The electronic contributions were computed analytically using density functional theory and multi-configurational self-consistent field method with extended basis sets for non-resonant excitation. The open-shell electronic system of nitrogen radical cation provides negative second-order optical nonlinearity, suggesting that the hyperpolarizability coefficient, {{γ }(2)}, in the non-resonant regime is mainly composed of combinations of virtual one-photon transitions rather than two-photon transitions. The second-order optical properties of nitrogen radical cation have been calculated as a function of bond length starting with the neutral molecular geometry (S0 minimum) and stretching the N-N triple bond, reaching the ionic D0 relaxed geometry all the way toward dissociation limit, to investigate the effect of internuclear bond distance on second-order hyperpolarizability.
An efficient second-order projection method for viscous incompressible flow
NASA Astrophysics Data System (ADS)
Bell, J. B.; Howell, L. H.; Colella, P.
1991-04-01
In this paper we describe a second-order projection method for the time-dependent, incompressible Navier-Stokes equations. The method is a second-order fractional step scheme in which one first solves diffusion-convection equations to determine intermediate velocities which are then projected onto the space of divergence-free vector fields. The diffusion-convection step uses a specialized second-order Godunov method for differencing the nonlinear convective terms that is conservative and free-stream preserving and provides a robust treatment of the nonlinearities at high Reynolds number. The projection is based on cell-centered centered difference approximations to divergence and gradient operators with the resulting linear system solved using a multigrid relaxation scheme. We apply the method to vortex spindown in a box to validate the numerical convergence of the method and to measure its overall performance.
Feasibility of a second-order gravitational red-shift experiment
NASA Technical Reports Server (NTRS)
Jaffe, J.; Vessot, R. F. C.
1976-01-01
The number of gravitation experiments undertaken since the advent of Einstein's theory of gravitation is quite small, with, so far, only the famous perihelion-advance experiment and a recent lunar-laser-ranging experiment being capable of measuring a nonlinear, second-order effect. It now appears that another distinct test of the second-order term may be feasible through the use of very stable atomic clocks. This experiment, which would measure the second-order gravitational red-shift, is a bona fide test of the field equations of gravity, not just a test of the underlying principle of equivalence. The nature of such an experiment, the basic equations, model-orbit calculations, and some tracking-accuracy requirements are presented. It is concluded that current space-probe tracking capabilities cannot determine all the necessary orbital parameters with sufficient accuracy for this experiment at the present time.
First- or second-order transition in the melting of repeat sequence DNA.
Chen, Y Z; Prohofsky, E W
1994-01-01
Both theoretical analysis and observation of the continuity of the melted fraction of base pairs indicate that the melting transition in DNA is second order. Analysis of the salt dependence of the transition by polyelectrolyte limiting laws, however, has first-order dynamics imbedded in the analysis. This paper proposes that the observation taken to be a latent heat of melting in the limiting law analysis could instead be a specific heat anomaly associated with a second-order transition. The limiting laws can be reconstructed based on a second-order transition with a specific heat anomaly. The T2M dependence of this excess heat is also consistent with its being a specific heat anomaly of a system displaying classical critical behavior. Classical critical behavior indicates that theoretical mean field approaches such as MSPA should be particularly appropriate to helix melting studies. PMID:8130338
Chen, Jiale; Gao, Zhe
2013-08-15
The second-order velocity distribution function was calculated from the second-order rf kinetic theory [Jaeger et al., Phys. Plasmas 7, 641 (2000)]. However, the nonresonant ponderomotive force in the radial direction derived from the theory is inconsistent with that from the fluid theory. The inconsistency arises from that the multiple-timescale-separation assumption fails when the second-order Vlasov equation is directly integrated along unperturbed particle orbits. A slowly ramped wave field including an adiabatic turn-on process is applied in the modified kinetic theory in this paper. Since this modification leads only to additional reactive/nonresonant response relevant with the secular resonant response from the previous kinetic theory, the correct nonresonant ponderomotive force can be obtained while all the resonant moments remain unchanged.
Free-Form Region Description with Second-Order Pooling.
Carreira, João; Caseiro, Rui; Batista, Jorge; Sminchisescu, Cristian
2015-06-01
Semantic segmentation and object detection are nowadays dominated by methods operating on regions obtained as a result of a bottom-up grouping process (segmentation) but use feature extractors developed for recognition on fixed-form (e.g. rectangular) patches, with full images as a special case. This is most likely suboptimal. In this paper we focus on feature extraction and description over free-form regions and study the relationship with their fixed-form counterparts. Our main contributions are novel pooling techniques that capture the second-order statistics of local descriptors inside such free-form regions. We introduce second-order generalizations of average and max-pooling that together with appropriate non-linearities, derived from the mathematical structure of their embedding space, lead to state-of-the-art recognition performance in semantic segmentation experiments without any type of local feature coding. In contrast, we show that codebook-based local feature coding is more important when feature extraction is constrained to operate over regions that include both foreground and large portions of the background, as typical in image classification settings, whereas for high-accuracy localization setups, second-order pooling over free-form regions produces results superior to those of the winning systems in the contemporary semantic segmentation challenges, with models that are much faster in both training and testing. PMID:26357341
Second-order parametrized-post-Newtonian Lagrangian
Benacquista, M.J. )
1992-02-15
A many-body Lagrangian to second post-Newtonian order using an extension of the parametrized-post-Newtonian (PPN) formalism is introduced and the properties of new parameters are explored. A parametrized gauge transformation is developed to permit comparison with theories of gravity in a variety of different coordinate systems. A procedure to impose Lorentz invariance on a general second-order post-Newtonian Lagrangian is developed. The Lagrangian is then constrained to possess Lorentz invariance and a Lorentz-invariant'' gauge is introduced. The constrained Lagrangian is found to be described by ten new second-order PPN parameters. When the Lagrangian is further constrained to describe theories of gravity for which test particles move along geodesics, one of the ten new parameters is given entirely in terms of first-order PPN parameters, leaving only nine PPN parameters to describe the second-order gravitational interaction. A metric'' gauge is introduced which allows the metric to be easily found from the Lagrangian and is shown to reduce to the gauge associated with the canonical formalism of Arnowitt, Deser, and Misner when the general-relativity values of the PPN parameters are used.
NASA Astrophysics Data System (ADS)
Lu, Xu; Yang, Feng-Wei; Xie, Yi
2016-07-01
We analyze strong gravitational field time delay for photons coupled to the Weyl tensor in a Schwarzschild black hole. By making use of the method of strong deflection limit, we find that these time delays between relativistic images are significantly affected by polarization directions of such a coupling. A practical problem about determination of the polarization direction by observations is investigated. It is found that if the first and second relativistic images can be resolved, the measurement of time delay can more effectively improve detectability of the polarization direction.
Boyko, Vyacheslav M; Popovych, Roman O; Shapoval, Nataliya M
2013-01-01
Lie symmetries of systems of second-order linear ordinary differential equations with constant coefficients are exhaustively described over both the complex and real fields. The exact lower and upper bounds for the dimensions of the maximal Lie invariance algebras possessed by such systems are obtained using an effective algebraic approach. PMID:23564972
Sugimoto, Satoru; Toki, Hiroshi; Ikeda, Kiyomi
2008-04-29
We study the effect of the tensor force on nuclear structure with mean-field and beyond-mean-field methods. An important correlation induced by the tensor force is a two-particle-two-hole (2p2h) correlation, which cannot be treated with a standard mean-field method. To treat the 2p2h tensor correlation, we develop a new framework [charge- and parity-projected Hartree-Fock (CPPHF) method], which is a beyond-mean-field method. In the CPPHF method, we introduce single-particle states with parity and charge mixing. The parity and charge projections are performed on a total wave function before variation. We apply the CPPHF method to oxygen isotopes including neutron-rich ones. The potential energy from the tensor force has the same order of magnitude as that from the LS force and becomes smaller with neutron number, which indicates that excess neutrons do not contribute to the 2p2h tensor correlation significantly. We also study the effect of the tensor force on spin-orbit-splitting (ls-splitting) in a neutron-rich fluorine isotope {sup 23}F. The tensor force reduces the ls-splitting for the proton d-orbits by about 3 MeV. This effect is important to reproduce the experimental value.
NASA Astrophysics Data System (ADS)
Sameer, M. Ikhdair; Majid, Hamzavi
2013-09-01
Approximate analytical solutions of the Dirac equation for Tietz—Hua (TH) potential including Coulomb-like tensor (CLT) potential with arbitrary spin—orbit quantum number κ are obtained within the Pekeris approximation scheme to deal with the spin—orbit coupling terms κ(κ ± 1)r-2. Under the exact spin and pseudospin symmetric limitation, bound state energy eigenvalues and associated unnormalized two-component wave functions of the Dirac particle in the field of both attractive and repulsive TH potential with tensor potential are found using the parametric Nikiforov—Uvarov (NU) method. The cases of the Morse oscillator with tensor potential, the generalized Morse oscillator with tensor potential, and the non-relativistic limits have been investigated.
Determination of robust stability margin for second-order systems
NASA Technical Reports Server (NTRS)
Chuang, C.-H.; Kau, C.-T.; Juang, Jer-Nan
1992-01-01
Robust stabilization of uncertain systems has been extensively investigated and the stability test for the whole set of uncertain parameters has been reduced to a finite number of test points, four points for the characteristic polynomial with independent coefficients. As a result the robust stability margin can be determined using a reasonable amount of computation. It is impossible to apply the results of the test to a practical system as the coefficients of the characteristic polynomial for a physical system are usually functions of uncertain parameters. However, many physical systems may be represented by a second-order mass-spring-damper system with a special multilinear form in its characteristic polynomial. This paper investigates second-order mass-spring-damper systems and the reduction of the number of test points. It is shown that such a system with arbritrary compensators always has a multilinear characteristic polynomial. It is also shown that a line in the two-dimensional parameter space forms the boundary after the mapping of a multilinear characteristic polynomial and this interior extreme line forms a conic curve in the complex plane. The boundary of uncertain domain for a multilinear polynomial with two uncertainty parameters can be determined analytically using this curve, and the four sides image of a square of the uncertain parameter. Therefore, the stability margin may be determined by checking the intersections of the boundary with the zero point. A similar procedure can be used for second-order systems with more than two uncertainty parameters when parameter optimization is used in determining the boundary.
NASA Astrophysics Data System (ADS)
Studynka, J.; Chadima, M.; Hrouda, F.; Suza, P.
2013-12-01
Low-field magnetic susceptibility of diamagnetic and paramagnetic minerals as well as that of pure magnetite and all single-domain ferromagnetic (s.l.) minerals is field-independent. In contrast, magnetic susceptibility of multi-domain pyrrhotite, hematite and titanomagnetite may significantly depend on the field intensity. Hence, the AMS data acquired in various fields have a great potential to separate the magnetic fabric carried by the latter group of minerals from the whole-rock fabric. The determination of the field variation of AMS consist of separate measurements of each sample in several fields within the Rayleigh Law range and subsequent processing in which the field-independent and field-dependent susceptibility tensors are calculated. The disadvantage of this technique is that each sample must be measured several times in various positions, which is relatively laborious and time consuming. Recently, a new 3D rotator was developed for the MFK1 Kappabridges which rotates the sample simultaneously about two axes with different velocities. The measurement is fully automated in such a way that, once the sample is mounted into the rotator, it requires no additional positioning to measure the full AMS tensor. The important advantage of the 3D rotator is that it enables to measure AMS in a sequence of pre-set field intensities without any operator manipulation. Whole procedure is computer-controlled and, once a sequence of measurements is finished, the acquired data are immediately processed and visualized. Examples of natural rocks demonstrating various types of field dependence of AMS are given.
Detection of a diffusive cloak via second-order statistics
NASA Astrophysics Data System (ADS)
Koirala, Milan; Yamilov, Alexey
2016-08-01
We propose a scheme to detect the diffusive cloak proposed by Schittny et al [Science 345, 427 (2014)]. We exploit the fact that diffusion of light is an approximation that disregards wave interference. The long-range contribution to intensity correlation is sensitive to locations of paths crossings and the interference inside the medium, allowing one to detect the size and position, including the depth, of the diffusive cloak. Our results also suggest that it is possible to separately manipulate the first- and the second-order statistics of wave propagation in turbid media.
Detection of a diffusive cloak via second-order statistics.
Koirala, Milan; Yamilov, Alexey
2016-08-15
We propose a scheme to detect the diffusive cloak proposed by Schittny et al. [Science345, 427 (2014).SCIEAS0036-807510.1126/science.1254524]. We exploit the fact that diffusion of light is an approximation that disregards wave interference. The long-range contribution to intensity correlation is sensitive to the locations of path crossings and the interference inside the medium, allowing one to detect the size and position, including the depth, of the diffusive cloak. Our results also suggest that it is possible to separately manipulate the first- and the second-order statistics of wave propagation in turbid media. PMID:27519108
Stochastic systems with delay: Perturbation theory for second order statistics
NASA Astrophysics Data System (ADS)
Frank, T. D.
2016-03-01
Within the framework of delay Fokker-Planck equations, a perturbation theoretical method is developed to determine second-order statistical quantities such as autocorrelation functions for stochastic systems with delay. Two variants of the perturbation theoretical approach are presented. The first variant is based on a non-local Fokker-Planck operator. The second variant requires to solve a Fokker-Planck equation with source term. It is shown that the two variants yield consistent results. The perturbation theoretical approaches are applied to study negative autocorrelations that are induced by feedback delays and mediated by the strength of the fluctuating forces that act on the feedback systems.
Supersonic second order analysis and optimization program user's manual
NASA Technical Reports Server (NTRS)
Clever, W. C.
1984-01-01
Approximate nonlinear inviscid theoretical techniques for predicting aerodynamic characteristics and surface pressures for relatively slender vehicles at supersonic and moderate hypersonic speeds were developed. Emphasis was placed on approaches that would be responsive to conceptual configuration design level of effort. Second order small disturbance theory was utilized to meet this objective. Numerical codes were developed for analysis and design of relatively general three dimensional geometries. Results from the computations indicate good agreement with experimental results for a variety of wing, body, and wing-body shapes. Case computational time of one minute on a CDC 176 are typical for practical aircraft arrangement.
Second-order kinetic Kohn-Sham lattice model
NASA Astrophysics Data System (ADS)
Solórzano, S.; Mendoza, M.; Herrmann, H. J.
2016-06-01
In this work, we introduce a semi-implicit second-order correction scheme to the kinetic Kohn-Sham lattice model. This approach is validated by performing realistic exchange-correlation energy calculations of atoms and dimers of the first two rows of the Periodic Table, finding good agreement with the expected values. Additionally, we simulate the ethane molecule, where we recover the bond lengths and compare the results with standard methods. Finally, we discuss the current applicability of pseudopotentials within the lattice kinetic Kohn-Sham approach.
P-stable boundary value methods for second order IVPs
NASA Astrophysics Data System (ADS)
Aceto, Lidia; Ghelardoni, Paolo; Magherini, Cecilia
2012-09-01
We introduce a family of Linear Multistep Methods used as Boundary Value Methods for the numerical solution of initial value problems for second order ordinary differential equations of special type. The aim is to obtain P-stable methods with arbitrary order of accuracy. This result allows to overcome the order barrier established by Lambert and Watson which limited to p = 2 the maximum order of a P-stable Linear Multistep Method. In addition, an extension of the methods in the Exponential Fitting framework is also considered.
A second-order impact model for forest fire regimes.
Maggi, Stefano; Rinaldi, Sergio
2006-09-01
We present a very simple "impact" model for the description of forest fires and show that it can mimic the known characteristics of wild fire regimes in savannas, boreal forests, and Mediterranean forests. Moreover, the distribution of burned biomasses in model generated fires resemble those of burned areas in numerous large forests around the world. The model has also the merits of being the first second-order model for forest fires and the first example of the use of impact models in the study of ecosystems. PMID:16723147
A critique of some recent second-order turbulence closure models for compressible boundary layers
NASA Technical Reports Server (NTRS)
Rubesin, M. W.; Crisalli, A. J.; Horstman, C. C.; Acharya, M.; Lanfranco, M. J.
1977-01-01
Computations based on two recently developed second-order turbulence closure models are compared with a series of boundary-layer experiments and with predictions of these experiments using an algebraic mixing length model. One of the models employs an eddy viscosity, whereas the other evaluates components of the Reynolds stress tensor. For flat plates, the computations are compared with the van Driest skin-friction transformation to assess the handling of compressibility. For boundary layers in pressure gradients, four experiments at Mach 4 and one at Mach 6.7 are used as the bases for comparison. In general, both models represent mean velocities and skin friction reasonably well, but represent the turbulence shear stress less accurately.
A-posteriori error estimation for second order mechanical systems
NASA Astrophysics Data System (ADS)
Ruiner, Thomas; Fehr, Jörg; Haasdonk, Bernard; Eberhard, Peter
2012-06-01
One important issue for the simulation of flexible multibody systems is the reduction of the flexible bodies degrees of freedom. As far as safety questions are concerned knowledge about the error introduced by the reduction of the flexible degrees of freedom is helpful and very important. In this work, an a-posteriori error estimator for linear first order systems is extended for error estimation of mechanical second order systems. Due to the special second order structure of mechanical systems, an improvement of the a-posteriori error estimator is achieved. A major advantage of the a-posteriori error estimator is that the estimator is independent of the used reduction technique. Therefore, it can be used for moment-matching based, Gramian matrices based or modal based model reduction techniques. The capability of the proposed technique is demonstrated by the a-posteriori error estimation of a mechanical system, and a sensitivity analysis of the parameters involved in the error estimation process is conducted.
Second-order modeling of arsenite transport in soils
NASA Astrophysics Data System (ADS)
Zhang, Hua; Magdi Selim, H.
2011-11-01
Rate limited processes including kinetic adsorption-desorption can greatly impact the fate and behavior of toxic arsenic compounds in heterogeneous soils. In this study, miscible displacement column experiments were carried out to investigate the extent of reactivity during transport of arsenite in soils. Arsenite breakthrough curves (BTCs) of Olivier and Windsor soils exhibited strong retardation with diffusive effluent fronts followed by slow release or tailing during leaching. Such behavior is indicative of the dominance of kinetic retention reactions for arsenite transport in the soil columns. Sharp decrease or increase in arsenite concentration in response to flow interruptions (stop-flow) further verified that non-equilibrium conditions are dominant. After some 40-60 pore volumes of continued leaching, 30-70% of the applied arsenite was retained by the soil in the columns. Furthermore, continued arsenite slow release for months was evident by the high levels of residual arsenite concentrations observed during leaching. In contrast, arsenite transport in a reference sand material exhibited no retention where complete mass recovery in the effluent solution was attained. A second-order model (SOM) which accounts for equilibrium, reversible, and irreversible retention mechanisms was utilized to describe arsenite transport results from the soil columns. Based on inverse and predictive modeling results, the SOM model successfully depicted arsenite BTCs from several soil columns. Based on inverse and predictive modeling results, a second-order model which accounts for kinetic reversible and irreversible reactions is recommended for describing arsenite transport in soils.
Rapid 3-D forward modeling of gravity and gravity gradient tensor fields
NASA Astrophysics Data System (ADS)
Longwei, C.; Dai, S.; Zhang, Q.
2014-12-01
Three-dimensional inversion are the key process in gravity exploration. In the commonly used scheme of inversion, the subsurface of the earth is usually divided into many small prism blocks (or grids) with variable density values. A key task in gravity inversion is to calculate the composite fields (gravity and gravity gradient tensor) generated by all these grids, this is known as forward modeling. In general forward modeling is memory-demanding and time-consuming. One scheme to rapidly calculate the fields is to implement it in Fourier domain and use fast Fourier transform algorithm. The advantage of the Fourier domain method is, obviously, much faster. However, the intrinsic edge effect of the Fourier domain method degrades the precision of the calculated fields. We have developed an innovative scheme to directly calculate the fields in spatial domain. There are two key points in this scheme. One key point is spatial discretization. Spatial convolution formula is discretized using an approach similar to normal difference method. A key idea during discretization is to use the analytical formula of a cubic prism, and this makes the resultant discrete formula have clear physical meaning: it embodies the superposition principle of the fields and is the exact formula to calculate the fields generated by all grids. The discretization only requires the grids have the same dimension in horizontal directions, and grids in different layers may have different dimension in vertical direction, and this offers more flexibility for inversion. Another key point is discrete convolution calculation. We invoke a high efficient two-dimensional discrete convolution algorithm, and it guarantees both time-saving and memory-saving. Its memory cost has the same order as the number of grids. Numerical test result shows that for a model with a dimension of 1000x1000x201 grids, it takes about 300s to calculate the fields on 1000x1000 field points in a personal computer with 3.4-GHz CPU
Hanbury Brown-Twiss interferometry and second-order correlations of inflaton quanta
Giovannini, Massimo
2011-01-15
The quantum theory of optical coherence is applied to the scrutiny of the statistical properties of the relic inflaton quanta. After adapting the description of the quantized scalar and tensor modes of the geometry to the analysis of intensity correlations, the normalized degrees of first-order and second-order coherence are computed in the concordance paradigm and are shown to encode faithfully the statistical properties of the initial quantum state. The strongly bunched curvature phonons are not only super-Poissonian but also superchaotic. Testable inequalities are derived in the limit of large-angular scales and can be physically interpreted in the light of the tenets of Hanbury Brown-Twiss interferometry. The quantum mechanical results are compared and contrasted with different situations including the one where intensity correlations are the result of a classical stochastic process. The survival of second-order correlations (not necessarily related to the purity of the initial quantum state) is addressed by defining a generalized ensemble where super-Poissonian statistics is an intrinsic property of the density matrix and turns out to be associated with finite volume effects which are expected to vanish in the thermodynamic limit.
Sachs-Wolfe at second order: the CMB bispectrum on large angular scales
Boubekeur, Lotfi; Creminelli, Paolo; D'Amico, Guido; Noreña, Jorge; Vernizzi, Filippo E-mail: creminel@ictp.it E-mail: norena@sissa.it
2009-08-01
We calculate the Cosmic Microwave Background anisotropy bispectrum on large angular scales in the absence of primordial non-Gaussianities, assuming exact matter dominance and extending at second order the classic Sachs-Wolfe result δT/T = Φ/3. The calculation is done in Poisson gauge. Besides intrinsic contributions calculated at last scattering, one must consider integrated effects. These are associated to lensing, and to the time dependence of the potentials (Rees-Sciama) and of the vector and tensor components of the metric generated at second order. The bispectrum is explicitly computed in the flat-sky approximation. It scales as l{sup −4} in the scale invariant limit and the shape dependence of its various contributions is represented in 3d plots. Although all the contributions to the bispectrum are parametrically of the same order, the full bispectrum is dominated by lensing. In the squeezed limit it corresponds to f{sub NL}{sup local} = −1/6−cos(2θ), where θ is the angle between the short and the long modes; the angle dependent contribution comes from lensing. In the equilateral limit it corresponds to f{sub NL}{sup equil} ≅ 3.13.
Tilt effects on moment tensor inversion in the near field of active volcanoes
NASA Astrophysics Data System (ADS)
van Driel, M.; Wassermann, J.; Pelties, C.; Schiemenz, A.; Igel, H.
2015-09-01
Dynamic tilts (rotational motion around horizontal axes) change the projection of local gravity onto the horizontal components of seismometers. This causes sensitivity of these components to tilt, especially at low frequencies. We analyse the consequences of this effect onto moment tensor inversion for very long period (vlp) events in the near field of active volcanoes on the basis of synthetic examples using the station distribution of a real deployed seismic network and the topography of Mt. Merapi volcano (Java, Indonesia). The examples show that for periods in the vlp range of 10-30 s tilt can have a strong effect on the moment tensor inversion, although its effect on the horizontal seismograms is significant only for few stations. We show that tilts can be accurately computed using the spectral element method and include them in the Green's functions. The (simulated) tilts might be largely influenced by strain-tilt coupling (stc). However, due to the frequency dependence of the tilt contribution to the horizontal seismograms, only the largest tilt signals affect the source inversion in the vlp frequency range. As these are less sensitive to stc than the weaker signals, the effect of stc can likely be neglected in this application. In the converse argument, this is not necessarily true for longer periods, where the horizontal seismograms are dominated by the tilt signal and rotational sensors would be necessary to account for it. As these are not yet commercially available, this study underlines the necessity for the development of such instruments.
3D extension of Tensorial Polar Decomposition. Application to (photo-)elasticity tensors
NASA Astrophysics Data System (ADS)
Desmorat, Rodrigue; Desmorat, Boris
2016-06-01
The orthogonalized harmonic decomposition of symmetric fourth-order tensors (i.e. having major and minor indicial symmetries, such as elasticity tensors) is completed by a representation of harmonic fourth-order tensors H by means of two second-order harmonic (symmetric deviatoric) tensors only. A similar decomposition is obtained for non-symmetric tensors (i.e. having minor indicial symmetry only, such as photo-elasticity tensors or elasto-plasticity tangent operators) introducing a fourth-order major antisymmetric traceless tensor Z. The tensor Z is represented by means of one harmonic second-order tensor and one antisymmetric second-order tensor only. Representations of totally symmetric (rari-constant), symmetric and major antisymmetric fourth-order tensors are simple particular cases of the proposed general representation. Closed-form expressions for tensor decomposition are given in the monoclinic case. Practical applications to elasticity and photo-elasticity monoclinic tensors are finally presented. xml:lang="fr"
Moreno-Torres, M.; Anguiano, M.; Grasso, M.; Van Giai, N.; Liang, H.; De Donno, V.
2010-06-15
Tensor effects in shell evolution are studied within the mean-field approach. Particular attention is paid to the analysis of the magic gaps in different regions of the nuclear chart, namely, Z,N=8, 20, and 28. Hartree-Fock calculations with Skyrme and Gogny interactions are performed where the tensor term has a zero and finite range, respectively. Results obtained with and without the tensor component are compared between them and with the experimental data, when available. To complete this analysis, the tensor effect is also investigated within the relativistic Hartree-Fock model, where the exchange of rho mesons and pions is taken into account. It turns out that the tensor effect in the evolution of the magic gaps can be more easily identified in the cases Z,N=8 and 20, whereas the interpretation of the effect is more complicated for Z or N= 28. Consequently, we indicate the regions defined by the magic numbers 8 and 20 as suitable for fitting the tensor parameters in a mean-field approach: We suggest to include explicitly the data associated to these gap evolutions in the fitting procedures. In general, with the parametrizations used in this work (which have not been fitted on these data), the mean-field results obtained with the tensor contribution do not reproduce the experimental trend, that is, the reduction of the gaps at 8 and 20 that is observed when going toward the drip lines. Since some of the considered nuclei have N=Z, a discussion will be devoted to the interpretation of the experimental data concerning these nuclei and to the Wigner-energy correction.
Second order distorted born approximation for backscattering from a layer of discrete random medium
NASA Technical Reports Server (NTRS)
Lang, Roger H.; Saatchi, Sasan S.
1993-01-01
In recent years there has been increasing interest in scattering and depolarization characteristics of the vegetation canopies. Scattering models applied to the microwave remote sensing of vegetation canopies showed that multiple scattering effects can be important in simulating the backscattering coefficients correctly. In particular, in most applications, the cross-polarized backscattering coefficients are often underestimated by single scattering models. Recently, there have been concerted efforts to include the second order terms in the radiative transfer models of vegetation canopies in order to account for multiple scattering within the canopy. The coherent wave theory approach is extended to include multiple scattering effects to predict the coherent and incoherent backscattering contributions from a layer of vegetation canopy. The problem is initially formulated in terms of the exact equation for the correlation function of the field, i.e., the Bethe-Salpeter equation. Using fractional volume as a small parameter, a Foldy type approximation is made to obtain a more manageable correlation equation. This equation is iterated to obtain first and second order solutions. The iteration procedure assumes the variance of the field fluctuations are small compared to the coherent intensity. This assumption proved to be particularly successful in computing backscattering coefficients. First and second order backscattering coefficients are calculated from the iterants of the correlation equation. It is shown that the first order coefficients are the same as the distorted Born results used previously by the authors. These results contained enhancement terms in the direct-reflected contributions. The important contributions to second order backscattering are examined and interpreted in terms of scattering diagrams. Examples of situations in which second order backscattering coefficients are important are given.
A novel unsplit perfectly matched layer for the second-order acoustic wave equation.
Ma, Youneng; Yu, Jinhua; Wang, Yuanyuan
2014-08-01
When solving acoustic field equations by using numerical approximation technique, absorbing boundary conditions (ABCs) are widely used to truncate the simulation to a finite space. The perfectly matched layer (PML) technique has exhibited excellent absorbing efficiency as an ABC for the acoustic wave equation formulated as a first-order system. However, as the PML was originally designed for the first-order equation system, it cannot be applied to the second-order equation system directly. In this article, we aim to extend the unsplit PML to the second-order equation system. We developed an efficient unsplit implementation of PML for the second-order acoustic wave equation based on an auxiliary-differential-equation (ADE) scheme. The proposed method can benefit to the use of PML in simulations based on second-order equations. Compared with the existing PMLs, it has simpler implementation and requires less extra storage. Numerical results from finite-difference time-domain models are provided to illustrate the validity of the approach. PMID:24794509
Kim, Wangdo; Kohles, Sean S.
2009-01-01
Tracking tissue deformation is often hampered by material inhomogeneity, so local measurements tend to be insufficient thus lending to the necessity of full-field optical measurements. This study presents a novel approach to factoring heterogeneous deformation of soft and hard tissues in a fracture callus by introducing an anisotropic metric derived from the deformation gradient tensor (F). The deformation gradient tensor contains all the information available in a Green-Lagrange strain tensor, plus the rigid-body rotational components. A recent study [Bottlang et al., J. Biomech. 41(3), 2008] produced full-field strains within ovine fracture calluses acquired through the application of electronic speckle pattern interferometery (ESPI). The technique is based on infinitesimal strain approximation (Engineering Strain) whose scheme is not independent of rigid body rotation. In this work, for rotation extraction, the stretch and rotation tensors were separately determined from F by the polar decomposition theorem. Interfragmentary motions in a fracture gap were characterized by the two distinct mechanical factors (stretch and rotation) at each material point through full-field mapping. In the composite nature of bone and soft tissue, collagen arrangements are hypothesized such that fibers locally aligned with principal directions will stretch and fibers not aligned with the principal direction will rotate and stretch. This approach has revealed the deformation gradient tensor as an appropriate quantification of strain within callus bony and fibrous tissue via optical measurements. PMID:19647826
A new approach for second-order perturbation theory.
Tomlinson, David G; Asadchev, Andrey; Gordon, Mark S
2016-05-30
A new second-order perturbation theory (MP2) approach is presented for closed shell energy evaluations. The new algorithm has a significantly lower memory footprint, a lower FLOP (floating point operations) count, and a lower time to solution compared to previously implemented parallel MP2 methods in the GAMESS software package. Additionally, this algorithm features an adaptive approach for the disk/distributed memory storage of the MP2 amplitudes. The algorithm works well on a single workstation, small cluster, and large Cray cluster, and it allows one to perform large calculations with thousands of basis functions in a matter of hours on a single workstation. The same algorithm has been adapted for graphical processing unit (GPU) architecture. The performance of the new GPU algorithm is also discussed. © 2016 Wiley Periodicals, Inc. PMID:26940648
Absorbing boundary conditions for second-order hyperbolic equations
NASA Technical Reports Server (NTRS)
Jiang, Hong; Wong, Yau Shu
1989-01-01
A uniform approach to construct absorbing artificial boundary conditions for second-order linear hyperbolic equations is proposed. The nonlocal boundary condition is given by a pseudodifferential operator that annihilates travelling waves. It is obtained through the dispersion relation of the differential equation by requiring that the initial-boundary value problem admits the wave solutions travelling in one direction only. Local approximation of this global boundary condition yields an nth-order differential operator. It is shown that the best approximations must be in the canonical forms which can be factorized into first-order operators. These boundary conditions are perfectly absorbing for wave packets propagating at certain group velocities. A hierarchy of absorbing boundary conditions is derived for transonic small perturbation equations of unsteady flows. These examples illustrate that the absorbing boundary conditions are easy to derive, and the effectiveness is demonstrated by the numerical experiments.
Absorbing boundary conditions for second-order hyperbolic equations
NASA Technical Reports Server (NTRS)
Jiang, Hong; Wong, Yau Shu
1990-01-01
A uniform approach to construct absorbing artificial boundary conditions for second-order linear hyperbolic equations is proposed. The nonlocal boundary condition is given by a pseudodifferential operator that annihilates travelling waves. It is obtained through the dispersion relation of the differential equation by requiring that the initial-boundary value problem admits the wave solutions travelling in one direction only. Local approximation of this global boundary condition yields an nth-order differential operator. It is shown that the best approximations must be in the canonical forms which can be factorized into first-order operators. These boundary conditions are perfectly absorbing for wave packets propagating at certain group velocities. A hierarchy of absorbing boundary conditions is derived for transonic small perturbation equations of unsteady flows. These examples illustrate that the absorbing boundary conditions are easy to derive, and the effectiveness is demonstrated by the numerical experiments.
Analysis of implicit second-order upwind-biased stencils
NASA Technical Reports Server (NTRS)
Roberts, Thomas W.; Warren, Gary P.
1993-01-01
Truncation error and stability properties of several implicit upwind schemes for the two-dimensional Euler equations are examined. The schemes use linear data reconstruction methods to achieve second-order flux integrations where the implicit Jacobian operators are first order. The stability properties of the schemes are examined by a Von Neumann analysis of the linearized, constant-coefficient Euler equations. The choice of the data reconstruction method used to evaluate the flux integral has a dramatic effect on the convergence properties of the implicit solution method. In particular, the typical one-dimensional data reconstruction methods used with structured grids exhibit poor convergence properties compared to the unstructured grid method considered. Of the schemes examined, the one with the superior convergence properties is well-suited for both unstructured and structured grids, which has important implications for the design of implicit methods.
New implicitly solvable potential produced by second order shape invariance
Cannata, F.; Ioffe, M.V.; Kolevatova, E.V.; Nishnianidze, D.N.
2015-05-15
The procedure proposed recently by Bougie et al. (2010) to study the general form of shape invariant potentials in one-dimensional Supersymmetric Quantum Mechanics (SUSY QM) is generalized to the case of Higher Order SUSY QM with supercharges of second order in momentum. A new shape invariant potential is constructed by this method. It is singular at the origin, it grows at infinity, and its spectrum depends on the choice of connection conditions in the singular point. The corresponding Schrödinger equation is solved explicitly: the wave functions are constructed analytically, and the energy spectrum is defined implicitly via the transcendental equation which involves Confluent Hypergeometric functions. - Highlights: • New potential with 2nd order irreducible shape invariance was constructed. • The connection conditions at the singularity of potential were obtained. • The explicit expressions for all wave functions were derived. • The implicit equation for the energy spectrum was obtained.
Pump power dependence of second order correlation in nondegenerate SPDC
NASA Astrophysics Data System (ADS)
Kim, Charles; Kanner, Gary
2009-08-01
We observed the second order correlation peak for nondegenerate spontaneous parametric down conversion (SPDC) of a pulsed pump at 532 nm into 810 nm and 1550 nm entangled beams. We used a Si avalanche photodiode (APD) to detect the 810 nm photons, and an InGaAs APD to detect those at 1550 nm. We defined both a visibility and signal-to-noise ratio (SNR) based on the data, which were obtained at various pump powers. In contrast to classical imaging systems, for which SNR increases monotonically with transmitted power, the SNR for the correlation peak in our setup exhibited a gradual decay as the pump power increased. We derived an empirical relation for the SNR, which was inversely proportional to the square root of pump power.
Second-order UV contamination in astronomical spectra
NASA Astrophysics Data System (ADS)
Gutierrez-Moreno, A.; Heathcote, S.; Moreno, H.; Hamuy, M.
1994-11-01
Observations of the planetary nebula Me 2-1 are used to discuss some effects of contamination produced by the second-order ultraviolet spectrum in the first-order red, longer than approximately 6200 A. It is found that this contamination is present with some instrumental setups, and that the contaminating lines are displaced in wavelength with respect to the predictions from diffraction theory by an amount which varies slightly for different setups. From the theoretical and empirical analysis of these displacements, it is concluded that the shift in wavelength is due to the presence of chromatic difference in magnification in spectrograph cameras including any transmission element, as in the case in the semi-solid Schmidt-Cassegrain cameras used in our observations.
Invariant classification of second-order conformally flat superintegrable systems
NASA Astrophysics Data System (ADS)
Capel, J. J.; Kress, J. M.
2014-12-01
In this paper we continue the work of Kalnins et al in classifying all second-order conformally-superintegrable (Laplace-type) systems over conformally flat spaces, using tools from algebraic geometry and classical invariant theory. The results obtained show, through Stäckel equivalence, that the list of known nondegenerate superintegrable systems over three-dimensional conformally flat spaces is complete. In particular, a seven-dimensional manifold is determined such that each point corresponds to a conformal class of superintegrable systems. This manifold is foliated by the nonlinear action of the conformal group in three dimensions. Two systems lie in the same conformal class if and only if they lie in the same leaf of the foliation. This foliation is explicitly described using algebraic varieties formed from representations of the conformal group. The proof of these results rely heavily on Gröbner basis calculations using the computer algebra software packages Maple and Singular.
WEAK GALERKIN METHODS FOR SECOND ORDER ELLIPTIC INTERFACE PROBLEMS.
Mu, Lin; Wang, Junping; Wei, Guowei; Ye, Xiu; Zhao, Shan
2013-10-01
Weak Galerkin methods refer to general finite element methods for partial differential equations (PDEs) in which differential operators are approximated by their weak forms as distributions. Such weak forms give rise to desirable flexibilities in enforcing boundary and interface conditions. A weak Galerkin finite element method (WG-FEM) is developed in this paper for solving elliptic PDEs with discontinuous coefficients and interfaces. Theoretically, it is proved that high order numerical schemes can be designed by using the WG-FEM with polynomials of high order on each element. Extensive numerical experiments have been carried to validate the WG-FEM for solving second order elliptic interface problems. High order of convergence is numerically confirmed in both L 2 and L ∞ norms for the piecewise linear WG-FEM. Special attention is paid to solve many interface problems, in which the solution possesses a certain singularity due to the nonsmoothness of the interface. A challenge in research is to design nearly second order numerical methods that work well for problems with low regularity in the solution. The best known numerical scheme in the literature is of order [Formula: see text] to [Formula: see text] for the solution itself in L ∞ norm. It is demonstrated that the WG-FEM of the lowest order, i.e., the piecewise constant WG-FEM, is capable of delivering numerical approximations that are of order [Formula: see text] to [Formula: see text] in the L ∞ norm for C (1) or Lipschitz continuous interfaces associated with a C (1) or H (2) continuous solution. PMID:24072935
Localization and mass spectra of various matter fields on scalar-tensor brane
Xie, Qun-Ying; Zhao, Zhen-Hua; Zhong, Yi; Yang, Jie; Zhou, Xiang-Nan
2015-03-10
Recently, a new scalar-tensor braneworld model was presented in [http://dx.doi.org/10.1103/PhysRevD.86.127502]. It not only solves the gauge hierarchy problem but also reproduces a correct Friedmann-like equation on the brane. In this new model, there are two different brane solutions, for which the mass spectra of gravity on the brane are the same. In this paper, we investigate localization and mass spectra of various bulk matter fields (i.e., scalar, vector, Kalb-Ramond, and fermion fields) on the brane. It is shown that the zero modes of all the matter fields can be localized on the positive tension brane under some conditions, and the mass spectra of each kind of bulk matter field for the two brane solutions are different except for some special cases, which implies that the two brane solutions are not physically equivalent. When the coupling constants between the dilaton and bulk matter fields take special values, the mass spectra for both solutions are the same, and the scalar and vector zero modes are localized on the negative tension brane, while the KR zero mode is still localized on the positive tension brane.