Approximate Model of Zone Sedimentation

NASA Astrophysics Data System (ADS)

Dzianik, František

2011-12-01

The process of zone sedimentation is affected by many factors that are not possible to express analytically. For this reason, the zone settling is evaluated in practice experimentally or by application of an empirical mathematical description of the process. The paper presents the development of approximate model of zone settling, i.e. the general function which should properly approximate the behaviour of the settling process within its entire range and at the various conditions. Furthermore, the specification of the model parameters by the regression analysis of settling test results is shown. The suitability of the model is reviewed by graphical dependencies and by statistical coefficients of correlation. The approximate model could by also useful on the simplification of process design of continual settling tanks and thickeners.

A Mathematica program for the approximate analytical solution to a nonlinear undamped Duffing equation by a new approximate approach

NASA Astrophysics Data System (ADS)

Wu, Dongmei; Wang, Zhongcheng

2006-03-01

According to Mickens [R.E. Mickens, Comments on a Generalized Galerkin's method for non-linear oscillators, J. Sound Vib. 118 (1987) 563], the general HB (harmonic balance) method is an approximation to the convergent Fourier series representation of the periodic solution of a nonlinear oscillator and not an approximation to an expansion in terms of a small parameter. Consequently, for a nonlinear undamped Duffing equation with a driving force Bcos(ωx), to find a periodic solution when the fundamental frequency is identical to ω, the corresponding Fourier series can be written as y˜(x)=∑n=1m acos[(2n-1)ωx]. How to calculate the coefficients of the Fourier series efficiently with a computer program is still an open problem. For HB method, by substituting approximation y˜(x) into force equation, expanding the resulting expression into a trigonometric series, then letting the coefficients of the resulting lowest-order harmonic be zero, one can obtain approximate coefficients of approximation y˜(x) [R.E. Mickens, Comments on a Generalized Galerkin's method for non-linear oscillators, J. Sound Vib. 118 (1987) 563]. But for nonlinear differential equations such as Duffing equation, it is very difficult to construct higher-order analytical approximations, because the HB method requires solving a set of algebraic equations for a large number of unknowns with very complex nonlinearities. To overcome the difficulty, forty years ago, Urabe derived a computational method for Duffing equation based on Galerkin procedure [M. Urabe, A. Reiter, Numerical computation of nonlinear forced oscillations by Galerkin's procedure, J. Math. Anal. Appl. 14 (1966) 107-140]. Dooren obtained an approximate solution of the Duffing oscillator with a special set of parameters by using Urabe's method [R. van Dooren, Stabilization of Cowell's classic finite difference method for numerical integration, J. Comput. Phys. 16 (1974) 186-192]. In this paper, in the frame of the general HB method

Born approximation in linear-time invariant system

NASA Astrophysics Data System (ADS)

Gumjudpai, Burin

2017-09-01

An alternative way of finding the LTI’s solution with the Born approximation, is investigated. We use Born approximation in the LTI and in the transformed LTI in form of Helmholtz equation. General solution are considered as infinite series or Feynman graph. Slow-roll approximation are explored. Transforming the LTI system into Helmholtz equation, approximated general solution can be found for any given forms of force with its initial value.

More on approximations of Poisson probabilities

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

Kao, C

1980-05-01

Calculation of Poisson probabilities frequently involves calculating high factorials, which becomes tedious and time-consuming with regular calculators. The usual way to overcome this difficulty has been to find approximations by making use of the table of the standard normal distribution. A new transformation proposed by Kao in 1978 appears to perform better for this purpose than traditional transformations. In the present paper several approximation methods are stated and compared numerically, including an approximation method that utilizes a modified version of Kao's transformation. An approximation based on a power transformation was found to outperform those based on the square-root type transformationsmore » as proposed in literature. The traditional Wilson-Hilferty approximation and Makabe-Morimura approximation are extremely poor compared with this approximation. 4 tables. (RWR)« less

Spline approximation, Part 1: Basic methodology

NASA Astrophysics Data System (ADS)

Ezhov, Nikolaj; Neitzel, Frank; Petrovic, Svetozar

2018-04-01

In engineering geodesy point clouds derived from terrestrial laser scanning or from photogrammetric approaches are almost never used as final results. For further processing and analysis a curve or surface approximation with a continuous mathematical function is required. In this paper the approximation of 2D curves by means of splines is treated. Splines offer quite flexible and elegant solutions for interpolation or approximation of "irregularly" distributed data. Depending on the problem they can be expressed as a function or as a set of equations that depend on some parameter. Many different types of splines can be used for spline approximation and all of them have certain advantages and disadvantages depending on the approximation problem. In a series of three articles spline approximation is presented from a geodetic point of view. In this paper (Part 1) the basic methodology of spline approximation is demonstrated using splines constructed from ordinary polynomials and splines constructed from truncated polynomials. In the forthcoming Part 2 the notion of B-spline will be explained in a unique way, namely by using the concept of convex combinations. The numerical stability of all spline approximation approaches as well as the utilization of splines for deformation detection will be investigated on numerical examples in Part 3.

Surface expression of the Chicxulub crater

PubMed

Pope, K O; Ocampo, A C; Kinsland, G L; Smith, R

1996-06-01

Analyses of geomorphic, soil, and topographic data from the northern Yucatan Peninsula, Mexico, confirm that the buried Chicxulub impact crater has a distinct surface expression and that carbonate sedimentation throughout the Cenozoic has been influenced by the crater. Late Tertiary sedimentation was mostly restricted to the region within the buried crater, and a semicircular moat existed until at least Pliocene time. The topographic expression of the crater is a series of features concentric with the crater. The most prominent is an approximately 83-km-radius trough or moat containing sinkholes (the Cenote ring). Early Tertiary surfaces rise abruptly outside the moat and form a stepped topography with an outer trough and ridge crest at radii of approximately 103 and approximately 129 km, respectively. Two discontinuous troughs lie within the moat at radii of approximately 41 and approximately 62 km. The low ridge between the inner troughs corresponds to the buried peak ring. The moat corresponds to the outer edge of the crater floor demarcated by a major ring fault. The outer trough and the approximately 62-km-radius inner trough also mark buried ring faults. The ridge crest corresponds to the topographic rim of the crater as modified by postimpact processes. These interpretations support previous findings that the principal impact basin has a diameter of approximately 180 km, but concentric, low-relief slumping extends well beyond this diameter and the eroded crater rim may extend to a diameter of approximately 260 km.

The Quantum Approximation Optimization Algorithm for MaxCut: A Fermionic View

NASA Technical Reports Server (NTRS)

Wang, Zhihui; Hadfield, Stuart; Jiang, Zhang; Rieffel, Eleanor G.

2017-01-01

Farhi et al. recently proposed a class of quantum algorithms, the Quantum Approximate Optimization Algorithm (QAOA), for approximately solving combinatorial optimization problems. A level-p QAOA circuit consists of steps in which a classical Hamiltonian, derived from the cost function, is applied followed by a mixing Hamiltonian. The 2p times for which these two Hamiltonians are applied are the parameters of the algorithm. As p increases, however, the parameter search space grows quickly. The success of the QAOA approach will depend, in part, on finding effective parameter-setting strategies. Here, we analytically and numerically study parameter setting for QAOA applied to MAXCUT. For level-1 QAOA, we derive an analytical expression for a general graph. In principle, expressions for higher p could be derived, but the number of terms quickly becomes prohibitive. For a special case of MAXCUT, the Ring of Disagrees, or the 1D antiferromagnetic ring, we provide an analysis for arbitrarily high level. Using a Fermionic representation, the evolution of the system under QAOA translates into quantum optimal control of an ensemble of independent spins. This treatment enables us to obtain analytical expressions for the performance of QAOA for any p. It also greatly simplifies numerical search for the optimal values of the parameters. By exploring symmetries, we identify a lower-dimensional sub-manifold of interest; the search effort can be accordingly reduced. This analysis also explains an observed symmetry in the optimal parameter values. Further, we numerically investigate the parameter landscape and show that it is a simple one in the sense of having no local optima.

Convergence and approximate calculation of average degree under different network sizes for decreasing random birth-and-death networks

NASA Astrophysics Data System (ADS)

Long, Yin; Zhang, Xiao-Jun; Wang, Kui

2018-05-01

In this paper, convergence and approximate calculation of average degree under different network sizes for decreasing random birth-and-death networks (RBDNs) are studied. First, we find and demonstrate that the average degree is convergent in the form of power law. Meanwhile, we discover that the ratios of the back items to front items of convergent reminder are independent of network link number for large network size, and we theoretically prove that the limit of the ratio is a constant. Moreover, since it is difficult to calculate the analytical solution of the average degree for large network sizes, we adopt numerical method to obtain approximate expression of the average degree to approximate its analytical solution. Finally, simulations are presented to verify our theoretical results.

Impulse approximation in nuclear pion production reactions: Absence of a one-body operator

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

Bolton, Daniel R.; Miller, Gerald A.

2011-06-15

The impulse approximation of pion production reactions is studied by developing a relativistic formalism, consistent with that used to define the nucleon-nucleon potential. For plane wave initial states we find that the usual one-body (1B) expression O{sub 1B} is replaced by O{sub 2B}=-iK(m{sub {pi}}/2)O{sub 1B}/m{sub {pi}}, where K(m{sub {pi}}/2) is the sum of all irreducible contributions to nucleon-nucleon scattering with energy transfer of m{sub {pi}}/2. We show that O{sub 2B}{approx_equal}O{sub 1B} for plane wave initial states. For distorted waves, we find that the usual operator is replaced with a sum of two-body operators that are well approximated by the operatormore » O{sub 2B}. Our new formalism solves the (previously ignored) problem of energy transfer forbidding a one-body impulse operator. Using a purely one pion exchange deuteron, the net result is that the impulse amplitude for np{yields}d{pi}{sup 0} at threshold is enhanced by a factor of approximately two. This amplitude is added to the larger ''rescattering'' amplitude and, although experimental data remain in disagreement, the theoretical prediction of the threshold cross section is brought closer to (and in agreement with) the data.« less

Finding gene regulatory network candidates using the gene expression knowledge base.

PubMed

Venkatesan, Aravind; Tripathi, Sushil; Sanz de Galdeano, Alejandro; Blondé, Ward; Lægreid, Astrid; Mironov, Vladimir; Kuiper, Martin

2014-12-10

Network-based approaches for the analysis of large-scale genomics data have become well established. Biological networks provide a knowledge scaffold against which the patterns and dynamics of 'omics' data can be interpreted. The background information required for the construction of such networks is often dispersed across a multitude of knowledge bases in a variety of formats. The seamless integration of this information is one of the main challenges in bioinformatics. The Semantic Web offers powerful technologies for the assembly of integrated knowledge bases that are computationally comprehensible, thereby providing a potentially powerful resource for constructing biological networks and network-based analysis. We have developed the Gene eXpression Knowledge Base (GeXKB), a semantic web technology based resource that contains integrated knowledge about gene expression regulation. To affirm the utility of GeXKB we demonstrate how this resource can be exploited for the identification of candidate regulatory network proteins. We present four use cases that were designed from a biological perspective in order to find candidate members relevant for the gastrin hormone signaling network model. We show how a combination of specific query definitions and additional selection criteria derived from gene expression data and prior knowledge concerning candidate proteins can be used to retrieve a set of proteins that constitute valid candidates for regulatory network extensions. Semantic web technologies provide the means for processing and integrating various heterogeneous information sources. The GeXKB offers biologists such an integrated knowledge resource, allowing them to address complex biological questions pertaining to gene expression. This work illustrates how GeXKB can be used in combination with gene expression results and literature information to identify new potential candidates that may be considered for extending a gene regulatory network.

Uniform analytic approximation of Wigner rotation matrices

NASA Astrophysics Data System (ADS)

Hoffmann, Scott E.

2018-02-01

We derive the leading asymptotic approximation, for low angle θ, of the Wigner rotation matrix elements, dm1m2 j(θ ) , uniform in j, m1, and m2. The result is in terms of a Bessel function of integer order. We numerically investigate the error for a variety of cases and find that the approximation can be useful over a significant range of angles. This approximation has application in the partial wave analysis of wavepacket scattering.

Technique for evaluation of the strong potential Born approximation for electron capture

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

Sil, N.C.; McGuire, J.H.

1985-04-01

A technique is presented for evaluating differential cross sections in the strong potential Born (SPB) approximation. Our final expression is expressed as a finite sum of one-dimensional integrals, expressible as a finite sum of derivatives of hypergeometric functions.

Asymptotic solution of the diffusion equation in slender impermeable tubes of revolution. I. The leading-term approximation

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

Traytak, Sergey D., E-mail: sergtray@mail.ru; Le STUDIUM; Semenov Institute of Chemical Physics RAS, 4 Kosygina St., 117977 Moscow

The anisotropic 3D equation describing the pointlike particles diffusion in slender impermeable tubes of revolution with cross section smoothly depending on the longitudinal coordinate is the object of our study. We use singular perturbations approach to find the rigorous asymptotic expression for the local particles concentration as an expansion in the ratio of the characteristic transversal and longitudinal diffusion relaxation times. The corresponding leading-term approximation is a generalization of well-known Fick-Jacobs approximation. This result allowed us to delineate the conditions on temporal and spatial scales under which the Fick-Jacobs approximation is valid. A striking analogy between solution of our problemmore » and the method of inner-outer expansions for low Knudsen numbers gas kinetic theory is established. With the aid of this analogy we clarify the physical and mathematical meaning of the obtained results.« less

Parameter inference in small world network disease models with approximate Bayesian Computational methods

NASA Astrophysics Data System (ADS)

Walker, David M.; Allingham, David; Lee, Heung Wing Joseph; Small, Michael

2010-02-01

Small world network models have been effective in capturing the variable behaviour of reported case data of the SARS coronavirus outbreak in Hong Kong during 2003. Simulations of these models have previously been realized using informed “guesses” of the proposed model parameters and tested for consistency with the reported data by surrogate analysis. In this paper we attempt to provide statistically rigorous parameter distributions using Approximate Bayesian Computation sampling methods. We find that such sampling schemes are a useful framework for fitting parameters of stochastic small world network models where simulation of the system is straightforward but expressing a likelihood is cumbersome.

Approximate geodesic distances reveal biologically relevant structures in microarray data.

PubMed

Nilsson, Jens; Fioretos, Thoas; Höglund, Mattias; Fontes, Magnus

2004-04-12

Genome-wide gene expression measurements, as currently determined by the microarray technology, can be represented mathematically as points in a high-dimensional gene expression space. Genes interact with each other in regulatory networks, restricting the cellular gene expression profiles to a certain manifold, or surface, in gene expression space. To obtain knowledge about this manifold, various dimensionality reduction methods and distance metrics are used. For data points distributed on curved manifolds, a sensible distance measure would be the geodesic distance along the manifold. In this work, we examine whether an approximate geodesic distance measure captures biological similarities better than the traditionally used Euclidean distance. We computed approximate geodesic distances, determined by the Isomap algorithm, for one set of lymphoma and one set of lung cancer microarray samples. Compared with the ordinary Euclidean distance metric, this distance measure produced more instructive, biologically relevant, visualizations when applying multidimensional scaling. This suggests the Isomap algorithm as a promising tool for the interpretation of microarray data. Furthermore, the results demonstrate the benefit and importance of taking nonlinearities in gene expression data into account.

Double power series method for approximating cosmological perturbations

NASA Astrophysics Data System (ADS)

Wren, Andrew J.; Malik, Karim A.

2017-04-01

We introduce a double power series method for finding approximate analytical solutions for systems of differential equations commonly found in cosmological perturbation theory. The method was set out, in a noncosmological context, by Feshchenko, Shkil' and Nikolenko (FSN) in 1966, and is applicable to cases where perturbations are on subhorizon scales. The FSN method is essentially an extension of the well known Wentzel-Kramers-Brillouin (WKB) method for finding approximate analytical solutions for ordinary differential equations. The FSN method we use is applicable well beyond perturbation theory to solve systems of ordinary differential equations, linear in the derivatives, that also depend on a small parameter, which here we take to be related to the inverse wave-number. We use the FSN method to find new approximate oscillating solutions in linear order cosmological perturbation theory for a flat radiation-matter universe. Together with this model's well-known growing and decaying Mészáros solutions, these oscillating modes provide a complete set of subhorizon approximations for the metric potential, radiation and matter perturbations. Comparison with numerical solutions of the perturbation equations shows that our approximations can be made accurate to within a typical error of 1%, or better. We also set out a heuristic method for error estimation. A Mathematica notebook which implements the double power series method is made available online.

Structural physical approximation for the realization of the optimal singlet fraction with two measurements

NASA Astrophysics Data System (ADS)

Adhikari, Satyabrata

2018-04-01

Structural physical approximation (SPA) has been exploited to approximate nonphysical operation such as partial transpose. It has already been studied in the context of detection of entanglement and found that if the minimum eigenvalue of SPA to partial transpose is less than 2/9 then the two-qubit state is entangled. We find application of SPA to partial transpose in the estimation of the optimal singlet fraction. We show that the optimal singlet fraction can be expressed in terms of the minimum eigenvalue of SPA to partial transpose. We also show that the optimal singlet fraction can be realized using Hong-Ou-Mandel interferometry with only two detectors. Further we have shown that the generated hybrid entangled state between a qubit and a binary coherent state can be used as a resource state in quantum teleportation.

Spline approximations for nonlinear hereditary control systems

NASA Technical Reports Server (NTRS)

Daniel, P. L.

1982-01-01

A sline-based approximation scheme is discussed for optimal control problems governed by nonlinear nonautonomous delay differential equations. The approximating framework reduces the original control problem to a sequence of optimization problems governed by ordinary differential equations. Convergence proofs, which appeal directly to dissipative-type estimates for the underlying nonlinear operator, are given and numerical findings are summarized.

The Kubo-Greenwood formula as a result of the random phase approximation for the electrons of the metal

NASA Astrophysics Data System (ADS)

Ivliev, S. V.

2017-12-01

For calculation of short laser pulse absorption in metal the imaginary part of permittivity, which is simply related to the conductivity, is required. Currently to find the static and dynamic conductivity the Kubo-Greenwood formula is most commonly used. It describes the electromagnetic energy absorption in the one-electron approach. In the present study, this formula is derived directly from the expression for the permittivity expression in the random phase approximation, which in fact is equivalent to the method of the mean field. The detailed analysis of the role of electron-electron interaction in the calculation of the matrix elements of the velocity operator is given. It is shown that in the one-electron random phase approximation the single-particle conductive electron wave functions in the field of fixed ions should be used. The possibility of considering the exchange and correlation effects by means of an amendment to a local function field is discussed.

Lognormal Approximations of Fault Tree Uncertainty Distributions.

PubMed

El-Shanawany, Ashraf Ben; Ardron, Keith H; Walker, Simon P

2018-01-26

Fault trees are used in reliability modeling to create logical models of fault combinations that can lead to undesirable events. The output of a fault tree analysis (the top event probability) is expressed in terms of the failure probabilities of basic events that are input to the model. Typically, the basic event probabilities are not known exactly, but are modeled as probability distributions: therefore, the top event probability is also represented as an uncertainty distribution. Monte Carlo methods are generally used for evaluating the uncertainty distribution, but such calculations are computationally intensive and do not readily reveal the dominant contributors to the uncertainty. In this article, a closed-form approximation for the fault tree top event uncertainty distribution is developed, which is applicable when the uncertainties in the basic events of the model are lognormally distributed. The results of the approximate method are compared with results from two sampling-based methods: namely, the Monte Carlo method and the Wilks method based on order statistics. It is shown that the closed-form expression can provide a reasonable approximation to results obtained by Monte Carlo sampling, without incurring the computational expense. The Wilks method is found to be a useful means of providing an upper bound for the percentiles of the uncertainty distribution while being computationally inexpensive compared with full Monte Carlo sampling. The lognormal approximation method and Wilks's method appear attractive, practical alternatives for the evaluation of uncertainty in the output of fault trees and similar multilinear models. © 2018 Society for Risk Analysis.

Exact and approximate stochastic simulation of intracellular calcium dynamics.

PubMed

Wieder, Nicolas; Fink, Rainer H A; Wegner, Frederic von

2011-01-01

In simulations of chemical systems, the main task is to find an exact or approximate solution of the chemical master equation (CME) that satisfies certain constraints with respect to computation time and accuracy. While Brownian motion simulations of single molecules are often too time consuming to represent the mesoscopic level, the classical Gillespie algorithm is a stochastically exact algorithm that provides satisfying results in the representation of calcium microdomains. Gillespie's algorithm can be approximated via the tau-leap method and the chemical Langevin equation (CLE). Both methods lead to a substantial acceleration in computation time and a relatively small decrease in accuracy. Elimination of the noise terms leads to the classical, deterministic reaction rate equations (RRE). For complex multiscale systems, hybrid simulations are increasingly proposed to combine the advantages of stochastic and deterministic algorithms. An often used exemplary cell type in this context are striated muscle cells (e.g., cardiac and skeletal muscle cells). The properties of these cells are well described and they express many common calcium-dependent signaling pathways. The purpose of the present paper is to provide an overview of the aforementioned simulation approaches and their mutual relationships in the spectrum ranging from stochastic to deterministic algorithms.

Fostering Formal Commutativity Knowledge with Approximate Arithmetic

PubMed Central

Hansen, Sonja Maria; Haider, Hilde; Eichler, Alexandra; Godau, Claudia; Frensch, Peter A.; Gaschler, Robert

2015-01-01

How can we enhance the understanding of abstract mathematical principles in elementary school? Different studies found out that nonsymbolic estimation could foster subsequent exact number processing and simple arithmetic. Taking the commutativity principle as a test case, we investigated if the approximate calculation of symbolic commutative quantities can also alter the access to procedural and conceptual knowledge of a more abstract arithmetic principle. Experiment 1 tested first graders who had not been instructed about commutativity in school yet. Approximate calculation with symbolic quantities positively influenced the use of commutativity-based shortcuts in formal arithmetic. We replicated this finding with older first graders (Experiment 2) and third graders (Experiment 3). Despite the positive effect of approximation on the spontaneous application of commutativity-based shortcuts in arithmetic problems, we found no comparable impact on the application of conceptual knowledge of the commutativity principle. Overall, our results show that the usage of a specific arithmetic principle can benefit from approximation. However, the findings also suggest that the correct use of certain procedures does not always imply conceptual understanding. Rather, the conceptual understanding of commutativity seems to lag behind procedural proficiency during elementary school. PMID:26560311

Radiative transfer in falling snow: A two-stream approximation

NASA Astrophysics Data System (ADS)

Koh, Gary

1989-04-01

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

On approximate formulas for the electrostatic force between two conducting spheres

NASA Astrophysics Data System (ADS)

Sliško, Josip; Brito-Orta, Raúl A.

1998-04-01

A series expression for the electrostatic force between two charged conducting spheres having equal radii and charges is derived using the method of electrical images. This expression is a special case of that for two spheres with arbitrary charges and radii, found by Maxwell using zonal harmonics. Keeping in mind the use of approximate formulas for the interpretation of classroom measurements of the electrostatic force between spheres, we comment on two incorrect approximate formulas and examine the contribution of the first few non-Coulomb terms of the correct formula by comparing with values obtained using a computational approach.

Rational trigonometric approximations using Fourier series partial sums

NASA Technical Reports Server (NTRS)

Geer, James F.

1993-01-01

A class of approximations (S(sub N,M)) to a periodic function f which uses the ideas of Pade, or rational function, approximations based on the Fourier series representation of f, rather than on the Taylor series representation of f, is introduced and studied. Each approximation S(sub N,M) is the quotient of a trigonometric polynomial of degree N and a trigonometric polynomial of degree M. The coefficients in these polynomials are determined by requiring that an appropriate number of the Fourier coefficients of S(sub N,M) agree with those of f. Explicit expressions are derived for these coefficients in terms of the Fourier coefficients of f. It is proven that these 'Fourier-Pade' approximations converge point-wise to (f(x(exp +))+f(x(exp -)))/2 more rapidly (in some cases by a factor of 1/k(exp 2M)) than the Fourier series partial sums on which they are based. The approximations are illustrated by several examples and an application to the solution of an initial, boundary value problem for the simple heat equation is presented.

Quantum approximate optimization algorithm for MaxCut: A fermionic view

NASA Astrophysics Data System (ADS)

Wang, Zhihui; Hadfield, Stuart; Jiang, Zhang; Rieffel, Eleanor G.

2018-02-01

Farhi et al. recently proposed a class of quantum algorithms, the quantum approximate optimization algorithm (QAOA), for approximately solving combinatorial optimization problems (E. Farhi et al., arXiv:1411.4028; arXiv:1412.6062; arXiv:1602.07674). A level-p QAOA circuit consists of p steps; in each step a classical Hamiltonian, derived from the cost function, is applied followed by a mixing Hamiltonian. The 2 p times for which these two Hamiltonians are applied are the parameters of the algorithm, which are to be optimized classically for the best performance. As p increases, parameter optimization becomes inefficient due to the curse of dimensionality. The success of the QAOA approach will depend, in part, on finding effective parameter-setting strategies. Here we analytically and numerically study parameter setting for the QAOA applied to MaxCut. For the level-1 QAOA, we derive an analytical expression for a general graph. In principle, expressions for higher p could be derived, but the number of terms quickly becomes prohibitive. For a special case of MaxCut, the "ring of disagrees," or the one-dimensional antiferromagnetic ring, we provide an analysis for an arbitrarily high level. Using a fermionic representation, the evolution of the system under the QAOA translates into quantum control of an ensemble of independent spins. This treatment enables us to obtain analytical expressions for the performance of the QAOA for any p . It also greatly simplifies the numerical search for the optimal values of the parameters. By exploring symmetries, we identify a lower-dimensional submanifold of interest; the search effort can be accordingly reduced. This analysis also explains an observed symmetry in the optimal parameter values. Further, we numerically investigate the parameter landscape and show that it is a simple one in the sense of having no local optima.

Approximation Of Multi-Valued Inverse Functions Using Clustering And Sugeno Fuzzy Inference

NASA Technical Reports Server (NTRS)

Walden, Maria A.; Bikdash, Marwan; Homaifar, Abdollah

1998-01-01

Finding the inverse of a continuous function can be challenging and computationally expensive when the inverse function is multi-valued. Difficulties may be compounded when the function itself is difficult to evaluate. We show that we can use fuzzy-logic approximators such as Sugeno inference systems to compute the inverse on-line. To do so, a fuzzy clustering algorithm can be used in conjunction with a discriminating function to split the function data into branches for the different values of the forward function. These data sets are then fed into a recursive least-squares learning algorithm that finds the proper coefficients of the Sugeno approximators; each Sugeno approximator finds one value of the inverse function. Discussions about the accuracy of the approximation will be included.

Variational Gaussian approximation for Poisson data

NASA Astrophysics Data System (ADS)

Arridge, Simon R.; Ito, Kazufumi; Jin, Bangti; Zhang, Chen

2018-02-01

The Poisson model is frequently employed to describe count data, but in a Bayesian context it leads to an analytically intractable posterior probability distribution. In this work, we analyze a variational Gaussian approximation to the posterior distribution arising from the Poisson model with a Gaussian prior. This is achieved by seeking an optimal Gaussian distribution minimizing the Kullback-Leibler divergence from the posterior distribution to the approximation, or equivalently maximizing the lower bound for the model evidence. We derive an explicit expression for the lower bound, and show the existence and uniqueness of the optimal Gaussian approximation. The lower bound functional can be viewed as a variant of classical Tikhonov regularization that penalizes also the covariance. Then we develop an efficient alternating direction maximization algorithm for solving the optimization problem, and analyze its convergence. We discuss strategies for reducing the computational complexity via low rank structure of the forward operator and the sparsity of the covariance. Further, as an application of the lower bound, we discuss hierarchical Bayesian modeling for selecting the hyperparameter in the prior distribution, and propose a monotonically convergent algorithm for determining the hyperparameter. We present extensive numerical experiments to illustrate the Gaussian approximation and the algorithms.

CMB-lensing beyond the Born approximation

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

Marozzi, Giovanni; Fanizza, Giuseppe; Durrer, Ruth

2016-09-01

We investigate the weak lensing corrections to the cosmic microwave background temperature anisotropies considering effects beyond the Born approximation. To this aim, we use the small deflection angle approximation, to connect the lensed and unlensed power spectra, via expressions for the deflection angles up to third order in the gravitational potential. While the small deflection angle approximation has the drawback to be reliable only for multipoles ℓ ∼< 2500, it allows us to consistently take into account the non-Gaussian nature of cosmological perturbation theory beyond the linear level. The contribution to the lensed temperature power spectrum coming from the non-Gaussianmore » nature of the deflection angle at higher order is a new effect which has not been taken into account in the literature so far. It turns out to be the leading contribution among the post-Born lensing corrections. On the other hand, the effect is smaller than corrections coming from non-linearities in the matter power spectrum, and its imprint on CMB lensing is too small to be seen in present experiments.« less

Local approximation of a metapopulation's equilibrium.

PubMed

Barbour, A D; McVinish, R; Pollett, P K

2018-04-18

We consider the approximation of the equilibrium of a metapopulation model, in which a finite number of patches are randomly distributed over a bounded subset [Formula: see text] of Euclidean space. The approximation is good when a large number of patches contribute to the colonization pressure on any given unoccupied patch, and when the quality of the patches varies little over the length scale determined by the colonization radius. If this is the case, the equilibrium probability of a patch at z being occupied is shown to be close to [Formula: see text], the equilibrium occupation probability in Levins's model, at any point [Formula: see text] not too close to the boundary, if the local colonization pressure and extinction rates appropriate to z are assumed. The approximation is justified by giving explicit upper and lower bounds for the occupation probabilities, expressed in terms of the model parameters. Since the patches are distributed randomly, the occupation probabilities are also random, and we complement our bounds with explicit bounds on the probability that they are satisfied at all patches simultaneously.

A Survey of Techniques for Approximate Computing

DOE PAGES

Mittal, Sparsh

2016-03-18

Approximate computing trades off computation quality with the effort expended and as rising performance demands confront with plateauing resource budgets, approximate computing has become, not merely attractive, but even imperative. Here, we present a survey of techniques for approximate computing (AC). We discuss strategies for finding approximable program portions and monitoring output quality, techniques for using AC in different processing units (e.g., CPU, GPU and FPGA), processor components, memory technologies etc., and programming frameworks for AC. Moreover, we classify these techniques based on several key characteristics to emphasize their similarities and differences. Finally, the aim of this paper is tomore » provide insights to researchers into working of AC techniques and inspire more efforts in this area to make AC the mainstream computing approach in future systems.« less

Partially Coherent Scattering in Stellar Chromospheres. Part 4; Analytic Wing Approximations

NASA Technical Reports Server (NTRS)

Gayley, K. G.

1993-01-01

Simple analytic expressions are derived to understand resonance-line wings in stellar chromospheres and similar astrophysical plasmas. The results are approximate, but compare well with accurate numerical simulations. The redistribution is modeled using an extension of the partially coherent scattering approximation (PCS) which we term the comoving-frame partially coherent scattering approximation (CPCS). The distinction is made here because Doppler diffusion is included in the coherent/noncoherent decomposition, in a form slightly improved from the earlier papers in this series.

On the convergence of local approximations to pseudodifferential operators with applications

NASA Technical Reports Server (NTRS)

Hagstrom, Thomas

1994-01-01

We consider the approximation of a class pseudodifferential operators by sequences of operators which can be expressed as compositions of differential operators and their inverses. We show that the error in such approximations can be bounded in terms of L(1) error in approximating a convolution kernel, and use this fact to develop convergence results. Our main result is a finite time convergence analysis of the Engquist-Majda Pade approximants to the square root of the d'Alembertian. We also show that no spatially local approximation to this operator can be convergent uniformly in time. We propose some temporally local but spatially nonlocal operators with better long time behavior. These are based on Laguerre and exponential series.

Approximate reasoning using terminological models

NASA Technical Reports Server (NTRS)

Yen, John; Vaidya, Nitin

1992-01-01

Term Subsumption Systems (TSS) form a knowledge-representation scheme in AI that can express the defining characteristics of concepts through a formal language that has a well-defined semantics and incorporates a reasoning mechanism that can deduce whether one concept subsumes another. However, TSS's have very limited ability to deal with the issue of uncertainty in knowledge bases. The objective of this research is to address issues in combining approximate reasoning with term subsumption systems. To do this, we have extended an existing AI architecture (CLASP) that is built on the top of a term subsumption system (LOOM). First, the assertional component of LOOM has been extended for asserting and representing uncertain propositions. Second, we have extended the pattern matcher of CLASP for plausible rule-based inferences. Third, an approximate reasoning model has been added to facilitate various kinds of approximate reasoning. And finally, the issue of inconsistency in truth values due to inheritance is addressed using justification of those values. This architecture enhances the reasoning capabilities of expert systems by providing support for reasoning under uncertainty using knowledge captured in TSS. Also, as definitional knowledge is explicit and separate from heuristic knowledge for plausible inferences, the maintainability of expert systems could be improved.

Hamilton's Principle and Approximate Solutions to Problems in Classical Mechanics

ERIC Educational Resources Information Center

Schlitt, D. W.

1977-01-01

Shows how to use the Ritz method for obtaining approximate solutions to problems expressed in variational form directly from the variational equation. Application of this method to classical mechanics is given. (MLH)

Transcriptional regulation of IGF-I expression in skeletal muscle

NASA Technical Reports Server (NTRS)

McCall, G. E.; Allen, D. L.; Haddad, F.; Baldwin, K. M.

2003-01-01

The present study investigated the role of transcription in the regulation of insulin-like growth factor (IGF)-I expression in skeletal muscle. RT-PCR was used to determine endogenous expression of IGF-I pre-mRNA and mRNA in control (Con) and functionally overloaded (FO) rat plantaris. The transcriptional activities of five different-length IGF-I promoter fragments controlling transcription of a firefly luciferase (FLuc) reporter gene were tested in vitro by transfection of myoblasts or in vivo during FO by direct gene transfer into the plantaris. Increased endogenous IGF-I gene transcription during 7 days of plantaris FO was evidenced by an approximately 140-160% increase (P < 0.0001) in IGF-I pre-mRNA (a transcriptional marker). IGF-I mRNA expression also increased by approximately 90% (P < 0.0001), and it was correlated (R = 0.93; P < 0.0001) with the pre-mRNA increases. The three longest IGF-I exon 1 promoters induced reporter gene expression in proliferating C2C12 and L6E9 myoblasts. In differentiated L6E9 myotubes, promoter activity increased approximately two- to threefold over myoblasts. Overexpression of calcineurin and MyoD increased the activity of the -852/+192 promoter in C2C12 myotubes by approximately 5- and approximately 18-fold, respectively. However, FO did not induce these exogenous promoter fragments. Nevertheless, the present findings are consistent with the hypothesis that the IGF-I gene is transcriptionally regulated during muscle hypertrophy in vivo as evidenced by the induction of the endogenous IGF-I pre-mRNA during plantaris FO. The exon 1 promoter region of the IGF-I gene is sufficient to direct inducible expression in vitro; however, an in vivo response to FO may require elements outside the -852/+346 region of the exon 1 IGF-I promoter or features inherent to the endogenous IGF-I gene.

Origin of Quantum Criticality in Yb-Al-Au Approximant Crystal and Quasicrystal

NASA Astrophysics Data System (ADS)

Watanabe, Shinji; Miyake, Kazumasa

2016-06-01

To get insight into the mechanism of emergence of unconventional quantum criticality observed in quasicrystal Yb15Al34Au51, the approximant crystal Yb14Al35Au51 is analyzed theoretically. By constructing a minimal model for the approximant crystal, the heavy quasiparticle band is shown to emerge near the Fermi level because of strong correlation of 4f electrons at Yb. We find that charge-transfer mode between 4f electron at Yb on the 3rd shell and 3p electron at Al on the 4th shell in Tsai-type cluster is considerably enhanced with almost flat momentum dependence. The mode-coupling theory shows that magnetic as well as valence susceptibility exhibits χ ˜ T-0.5 for zero-field limit and is expressed as a single scaling function of the ratio of temperature to magnetic field T/B over four decades even in the approximant crystal when some condition is satisfied by varying parameters, e.g., by applying pressure. The key origin is clarified to be due to strong locality of the critical Yb-valence fluctuation and small Brillouin zone reflecting the large unit cell, giving rise to the extremely-small characteristic energy scale. This also gives a natural explanation for the quantum criticality in the quasicrystal corresponding to the infinite limit of the unit-cell size.

The consequences of improperly describing oscillator strengths beyond the electric dipole approximation

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

Lestrange, Patrick J.; Egidi, Franco; Li, Xiaosong, E-mail: xsli@uw.edu

2015-12-21

The interaction between a quantum mechanical system and plane wave light is usually modeled within the electric dipole approximation. This assumes that the intensity of the incident field is constant over the length of the system and transition probabilities are described in terms of the electric dipole transition moment. For short wavelength spectroscopies, such as X-ray absorption, the electric dipole approximation often breaks down. Higher order multipoles are then included to describe transition probabilities. The square of the magnetic dipole and electric quadrupole are often included, but this results in an origin-dependent expression for the oscillator strength. The oscillator strengthmore » can be made origin-independent if all terms through the same order in the wave vector are retained. We will show the consequences and potential pitfalls of using either of these two expressions. It is shown that the origin-dependent expression may violate the Thomas-Reiche-Kuhn sum rule and the origin-independent expression can result in negative transition probabilities.« less

The consequences of improperly describing oscillator strengths beyond the electric dipole approximation.

PubMed

Lestrange, Patrick J; Egidi, Franco; Li, Xiaosong

2015-12-21

The interaction between a quantum mechanical system and plane wave light is usually modeled within the electric dipole approximation. This assumes that the intensity of the incident field is constant over the length of the system and transition probabilities are described in terms of the electric dipole transition moment. For short wavelength spectroscopies, such as X-ray absorption, the electric dipole approximation often breaks down. Higher order multipoles are then included to describe transition probabilities. The square of the magnetic dipole and electric quadrupole are often included, but this results in an origin-dependent expression for the oscillator strength. The oscillator strength can be made origin-independent if all terms through the same order in the wave vector are retained. We will show the consequences and potential pitfalls of using either of these two expressions. It is shown that the origin-dependent expression may violate the Thomas-Reiche-Kuhn sum rule and the origin-independent expression can result in negative transition probabilities.

The consequences of improperly describing oscillator strengths beyond the electric dipole approximation

NASA Astrophysics Data System (ADS)

Lestrange, Patrick J.; Egidi, Franco; Li, Xiaosong

2015-12-01

The interaction between a quantum mechanical system and plane wave light is usually modeled within the electric dipole approximation. This assumes that the intensity of the incident field is constant over the length of the system and transition probabilities are described in terms of the electric dipole transition moment. For short wavelength spectroscopies, such as X-ray absorption, the electric dipole approximation often breaks down. Higher order multipoles are then included to describe transition probabilities. The square of the magnetic dipole and electric quadrupole are often included, but this results in an origin-dependent expression for the oscillator strength. The oscillator strength can be made origin-independent if all terms through the same order in the wave vector are retained. We will show the consequences and potential pitfalls of using either of these two expressions. It is shown that the origin-dependent expression may violate the Thomas-Reiche-Kuhn sum rule and the origin-independent expression can result in negative transition probabilities.

An approximation technique for predicting the transient response of a second order nonlinear equation

NASA Technical Reports Server (NTRS)

Laurenson, R. M.; Baumgarten, J. R.

1975-01-01

An approximation technique has been developed for determining the transient response of a nonlinear dynamic system. The nonlinearities in the system which has been considered appear in the system's dissipation function. This function was expressed as a second order polynomial in the system's velocity. The developed approximation is an extension of the classic Kryloff-Bogoliuboff technique. Two examples of the developed approximation are presented for comparative purposes with other approximation methods.

Accuracy of expressions for the fill factor of a solar cell in terms of open-circuit voltage and ideality factor

NASA Astrophysics Data System (ADS)

Leilaeioun, Mehdi; Holman, Zachary C.

2016-09-01

An approximate expression proposed by Green predicts the maximum obtainable fill factor (FF) of a solar cell from its open-circuit voltage (Voc). The expression was originally suggested for silicon solar cells that behave according to a single-diode model and, in addition to Voc, it requires an ideality factor as input. It is now commonly applied to silicon cells by assuming a unity ideality factor—even when the cells are not in low injection—as well as to non-silicon cells. Here, we evaluate the accuracy of the expression in several cases. In particular, we calculate the recombination-limited FF and Voc of hypothetical silicon solar cells from simulated lifetime curves, and compare the exact FF to that obtained with the approximate expression using assumed ideality factors. Considering cells with a variety of recombination mechanisms, wafer doping densities, and photogenerated current densities reveals the range of conditions under which the approximate expression can safely be used. We find that the expression is unable to predict FF generally: For a typical silicon solar cell under one-sun illumination, the error is approximately 6% absolute with an assumed ideality factor of 1. Use of the expression should thus be restricted to cells under very low or very high injection.

Approximate matching of structured motifs in DNA sequences.

PubMed

El-Mabrouk, Nadia; Raffinot, Mathieu; Duchesne, Jean-Eudes; Lajoie, Mathieu; Luc, Nicolas

2005-04-01

Several methods have been developed for identifying more or less complex RNA structures in a genome. All these methods are based on the search for conserved primary and secondary sub-structures. In this paper, we present a simple formal representation of a helix, which is a combination of sequence and folding constraints, as a constrained regular expression. This representation allows us to develop a well-founded algorithm that searches for all approximate matches of a helix in a genome. The algorithm is based on an alignment graph constructed from several copies of a pushdown automaton, arranged one on top of another. This is a first attempt to take advantage of the possibilities of pushdown automata in the context of approximate matching. The worst time complexity is O(krpn), where k is the error threshold, n the size of the genome, p the size of the secondary expression, and r its number of union symbols. We then extend the algorithm to search for pseudo-knots and secondary structures containing an arbitrary number of helices.

Finding Groups in Gene Expression Data

PubMed Central

2005-01-01

The vast potential of the genomic insight offered by microarray technologies has led to their widespread use since they were introduced a decade ago. Application areas include gene function discovery, disease diagnosis, and inferring regulatory networks. Microarray experiments enable large-scale, high-throughput investigations of gene activity and have thus provided the data analyst with a distinctive, high-dimensional field of study. Many questions in this field relate to finding subgroups of data profiles which are very similar. A popular type of exploratory tool for finding subgroups is cluster analysis, and many different flavors of algorithms have been used and indeed tailored for microarray data. Cluster analysis, however, implies a partitioning of the entire data set, and this does not always match the objective. Sometimes pattern discovery or bump hunting tools are more appropriate. This paper reviews these various tools for finding interesting subgroups. PMID:16046827

Testing the Ginzburg-Landau approximation for three-flavor crystalline color superconductivity

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

Mannarelli, Massimo; Sharma, Rishi; Rajagopal, Krishna

2006-06-01

It is an open challenge to analyze the crystalline color superconducting phases that may arise in cold dense, but not asymptotically dense, three-flavor quark matter. At present the only approximation within which it seems possible to compare the free energies of the myriad possible crystal structures is the Ginzburg-Landau approximation. Here, we test this approximation on a particularly simple 'crystal' structure in which there are only two condensates {approx}{delta}exp(iq{sub 2}{center_dot}r) and {approx}{delta}exp(iq{sub 3}{center_dot}r) whose position-space dependence is that of two plane waves with wave vectors q{sub 2} and q{sub 3} at arbitrary angles. For this case, we are able tomore » solve the mean-field gap equation without making a Ginzburg-Landau approximation. We find that the Ginzburg-Landau approximation works in the {delta}{yields}0 limit as expected, find that it correctly predicts that {delta} decreases with increasing angle between q{sub 2} and q{sub 3} meaning that the phase with q{sub 2} parallel q{sub 3} has the lowest free energy, and find that the Ginzburg-Landau approximation is conservative in the sense that it underestimates {delta} at all values of the angle between q{sub 2} and q{sub 3}.« less

Approximate analytical solution for induction heating of solid cylinders

DOE PAGES

Jankowski, Todd Andrew; Pawley, Norma Helen; Gonzales, Lindsey Michal; ...

2015-10-20

An approximate solution to the mathematical model for induction heating of a solid cylinder in a cylindrical induction coil is presented here. The coupled multiphysics model includes equations describing the electromagnetic field in the heated object, a heat transfer simulation to determine temperature of the heated object, and an AC circuit simulation of the induction heating power supply. A multiple-scale perturbation method is used to solve the multiphysics model. The approximate analytical solution yields simple closed-form expressions for the electromagnetic field and heat generation rate in the solid cylinder, for the equivalent impedance of the associated tank circuit, and formore » the frequency response of a variable frequency power supply driving the tank circuit. The solution developed here is validated by comparing predicted power supply frequency to both experimental measurements and calculated values from finite element analysis for heating of graphite cylinders in an induction furnace. The simple expressions from the analytical solution clearly show the functional dependence of the power supply frequency on the material properties of the load and the geometrical characteristics of the furnace installation. In conclusion, the expressions developed here provide physical insight into observations made during load signature analysis of induction heating.« less

libFLASM: a software library for fixed-length approximate string matching.

PubMed

Ayad, Lorraine A K; Pissis, Solon P P; Retha, Ahmad

2016-11-10

Approximate string matching is the problem of finding all factors of a given text that are at a distance at most k from a given pattern. Fixed-length approximate string matching is the problem of finding all factors of a text of length n that are at a distance at most k from any factor of length ℓ of a pattern of length m. There exist bit-vector techniques to solve the fixed-length approximate string matching problem in time [Formula: see text] and space [Formula: see text] under the edit and Hamming distance models, where w is the size of the computer word; as such these techniques are independent of the distance threshold k or the alphabet size. Fixed-length approximate string matching is a generalisation of approximate string matching and, hence, has numerous direct applications in computational molecular biology and elsewhere. We present and make available libFLASM, a free open-source C++ software library for solving fixed-length approximate string matching under both the edit and the Hamming distance models. Moreover we describe how fixed-length approximate string matching is applied to solve real problems by incorporating libFLASM into established applications for multiple circular sequence alignment as well as single and structured motif extraction. Specifically, we describe how it can be used to improve the accuracy of multiple circular sequence alignment in terms of the inferred likelihood-based phylogenies; and we also describe how it is used to efficiently find motifs in molecular sequences representing regulatory or functional regions. The comparison of the performance of the library to other algorithms show how it is competitive, especially with increasing distance thresholds. Fixed-length approximate string matching is a generalisation of the classic approximate string matching problem. We present libFLASM, a free open-source C++ software library for solving fixed-length approximate string matching. The extensive experimental results presented here

Analytic approximation for random muffin-tin alloys

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

Mills, R.; Gray, L.J.; Kaplan, T.

1983-03-15

The methods introduced in a previous paper under the name of ''traveling-cluster approximation'' (TCA) are applied, in a multiple-scattering approach, to the case of a random muffin-tin substitutional alloy. This permits the iterative part of a self-consistent calculation to be carried out entirely in terms of on-the-energy-shell scattering amplitudes. Off-shell components of the mean resolvent, needed for the calculation of spectral functions, are obtained by standard methods involving single-site scattering wave functions. The single-site TCA is just the usual coherent-potential approximation, expressed in a form particularly suited for iteration. A fixed-point theorem is proved for the general t-matrix TCA, ensuringmore » convergence upon iteration to a unique self-consistent solution with the physically essential Herglotz properties.« less

Best uniform approximation to a class of rational functions

NASA Astrophysics Data System (ADS)

Zheng, Zhitong; Yong, Jun-Hai

2007-10-01

We explicitly determine the best uniform polynomial approximation to a class of rational functions of the form 1/(x-c)2+K(a,b,c,n)/(x-c) on [a,b] represented by their Chebyshev expansion, where a, b, and c are real numbers, n-1 denotes the degree of the best approximating polynomial, and K is a constant determined by a, b, c, and n. Our result is based on the explicit determination of a phase angle [eta] in the representation of the approximation error by a trigonometric function. Moreover, we formulate an ansatz which offers a heuristic strategies to determine the best approximating polynomial to a function represented by its Chebyshev expansion. Combined with the phase angle method, this ansatz can be used to find the best uniform approximation to some more functions.

Approximation for limit cycles and their isochrons.

PubMed

Demongeot, Jacques; Françoise, Jean-Pierre

2006-12-01

Local analysis of trajectories of dynamical systems near an attractive periodic orbit displays the notion of asymptotic phase and isochrons. These notions are quite useful in applications to biosciences. In this note, we give an expression for the first approximation of equations of isochrons in the setting of perturbations of polynomial Hamiltonian systems. This method can be generalized to perturbations of systems that have a polynomial integral factor (like the Lotka-Volterra equation).

Histological findings and immunohistochemical surfactant protein A (SP-A) expression in asphyxia: its application in the diagnosis of drowning.

PubMed

Pérez-Cárceles, M D; Sibón, A; Vizcaya, M A; Osuna, E; Gómez-Zapata, M; Luna, A; Martínez-Díaz, F

2008-09-01

The histopathological alterations that permit the diagnosis of death by asphyxia are very unspecific, although pulmonary alterations are of great importance in this respect. The postmortem diagnosis of drowning, particularly, continues to be one of the most difficult in forensic pathology. The aim of this study is to jointly evaluate microscopic findings and immunohistochemical surfactant protein A (SP-A) expression in the upper and lower lobes of lungs in different causes of death, and their possible application to the diagnosis of drowning. We studied 120 cadavers from subjects with a mean age of 48.73 years (SD 19.45; range 2-86 years), and with a mean post-mortem interval of 30 hours (SD 39.59; range 3-216 hours). According to the scene, cause and circumstances of death, and autopsy findings, cases were classified into groups as follows: (a) drowning (n=47); (b) other asphyxia (n=44) and (c) other causes (n=29). In the upper and lower lobes of lungs, histological studies of H&E staining and immunohistochemical surfactant protein A expression were made. The presence and severity of congestion, haemorrhage and oedema, together with immunohistochemical SP-A expression, may have a diagnostic value in differentiating asphyxia and drowning from other causes of death, and drowning from other types of asphyxia. Our findings suggest that both lobes should be investigated to establish the diagnosis, although the findings in the upper lobe might be the most important for differentiating the exact cause of death.

Finding Meaning in Written Emotional Expression by Family Caregivers of Persons With Dementia.

PubMed

Butcher, Howard K; Gordon, Jean K; Ko, Ji Woon; Perkhounkova, Yelena; Cho, Jun Young; Rinner, Andrew; Lutgendorf, Susan

2016-12-01

This study tested the effect of written emotional expression on the ability to find meaning in caregiving and the effects of finding meaning on emotional state and psychological burden in 91 dementia family caregivers. In a pretest-posttest design, participants were randomly assigned to either an experimental or a comparison group. Experimental caregivers (n = 57) wrote about their deepest thoughts and feelings about caring for a family member with dementia, whereas those in the comparison group (n = 34) wrote about nonemotional topics. Results showed enhanced meaning-making abilities in experimental participants relative to comparison participants, particularly for those who used more positive emotion words. Improved meaning-making ability was in turn associated with psychological benefits at posttest, but experimental participants did not show significantly more benefit than comparison participants. We explore the mediating roles of the meaning-making process as well as some of the background characteristics of the individual caregivers and their caregiving environments. © The Author(s) 2016.

Differential global gene expression in red and white skeletal muscle

NASA Technical Reports Server (NTRS)

Campbell, W. G.; Gordon, S. E.; Carlson, C. J.; Pattison, J. S.; Hamilton, M. T.; Booth, F. W.

2001-01-01

The differences in gene expression among the fiber types of skeletal muscle have long fascinated scientists, but for the most part, previous experiments have only reported differences of one or two genes at a time. The evolving technology of global mRNA expression analysis was employed to determine the potential differential expression of approximately 3,000 mRNAs between the white quad (white muscle) and the red soleus muscle (mixed red muscle) of female ICR mice (30-35 g). Microarray analysis identified 49 mRNA sequences that were differentially expressed between white and mixed red skeletal muscle, including newly identified differential expressions between muscle types. For example, the current findings increase the number of known, differentially expressed mRNAs for transcription factors/coregulators by nine and signaling proteins by three. The expanding knowledge of the diversity of mRNA expression between white and mixed red muscle suggests that there could be quite a complex regulation of phenotype between muscles of different fiber types.

Highly Accurate Analytical Approximate Solution to a Nonlinear Pseudo-Oscillator

NASA Astrophysics Data System (ADS)

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

2017-07-01

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

A test of the adhesion approximation for gravitational clustering

NASA Technical Reports Server (NTRS)

Melott, Adrian L.; Shandarin, Sergei; Weinberg, David H.

1993-01-01

We quantitatively compare a particle implementation of the adhesion approximation to fully non-linear, numerical 'N-body' simulations. Our primary tool, cross-correlation of N-body simulations with the adhesion approximation, indicates good agreement, better than that found by the same test performed with the Zel-dovich approximation (hereafter ZA). However, the cross-correlation is not as good as that of the truncated Zel-dovich approximation (TZA), obtained by applying the Zel'dovich approximation after smoothing the initial density field with a Gaussian filter. We confirm that the adhesion approximation produces an excessively filamentary distribution. Relative to the N-body results, we also find that: (a) the power spectrum obtained from the adhesion approximation is more accurate than that from ZA or TZA, (b) the error in the phase angle of Fourier components is worse than that from TZA, and (c) the mass distribution function is more accurate than that from ZA or TZA. It appears that adhesion performs well statistically, but that TZA is more accurate dynamically, in the sense of moving mass to the right place.

A test of the adhesion approximation for gravitational clustering

NASA Technical Reports Server (NTRS)

Melott, Adrian L.; Shandarin, Sergei F.; Weinberg, David H.

1994-01-01

We quantitatively compare a particle implementation of the adhesion approximation to fully nonlinear, numerical 'N-body' simulations. Our primary tool, cross-correlation of N-body simulations with the adhesion approximation, indicates good agreement, better than that found by the same test performed with the Zel'dovich approximation (hereafter ZA). However, the cross-correlation is not as good as that of the truncated Zel'dovich approximation (TZA), obtained by applying the Zel'dovich approximation after smoothing the initial density field with a Gaussian filter. We confirm that the adhesion approximation produces an excessively filamentary distribution. Relative to the N-body results, we also find that: (a) the power spectrum obtained from the adhesion approximation is more accurate that that from ZA to TZA, (b) the error in the phase angle of Fourier components is worse that that from TZA, and (c) the mass distribution function is more accurate than that from ZA or TZA. It appears that adhesion performs well statistically, but that TZA is more accurate dynamically, in the sense of moving mass to the right place.

A new approximate sum rule for bulk alloy properties

NASA Technical Reports Server (NTRS)

Bozzolo, Guillermo; Ferrante, John

1991-01-01

A new, approximate sum rule is introduced for determining bulk properties of multicomponent systems, in terms of the pure components properties. This expression is applied for the study of lattice parameters, cohesive energies, and bulk moduli of binary alloys. The correct experimental trends (i.e., departure from average values) are predicted in all cases.

Elastic scattering of virtual photons via a quark loop in the double-logarithmic approximation

NASA Astrophysics Data System (ADS)

Ermolaev, B. I.; Ivanov, D. Yu.; Troyan, S. I.

2018-04-01

We calculate the amplitude of elastic photon-photon scattering via a single quark loop in the double-logarithmic approximation, presuming all external photons to be off-shell and unpolarized. At the same time we account for the running coupling effects. We consider this process in the forward kinematics at arbitrary relations between t and the external photon virtualities. We obtain explicit expressions for the photon-photon scattering amplitudes in all double-logarithmic kinematic regions. Then we calculate the small-x asymptotics of the obtained amplitudes and compare them with the parent amplitudes, thereby fixing the applicability regions of the asymptotics, i.e., fixing the applicability region for the nonvacuum Reggeons. We find that these Reggeons should be used at x <10-8 only.

Revised Thomas-Fermi approximation for singular potentials

NASA Astrophysics Data System (ADS)

Dufty, James W.; Trickey, S. B.

2016-08-01

Approximations for the many-fermion free-energy density functional that include the Thomas-Fermi (TF) form for the noninteracting part lead to singular densities for singular external potentials (e.g., attractive Coulomb). This limitation of the TF approximation is addressed here by a formal map of the exact Euler equation for the density onto an equivalent TF form characterized by a modified Kohn-Sham potential. It is shown to be a "regularized" version of the Kohn-Sham potential, tempered by convolution with a finite-temperature response function. The resulting density is nonsingular, with the equilibrium properties obtained from the total free-energy functional evaluated at this density. This new representation is formally exact. Approximate expressions for the regularized potential are given to leading order in a nonlocality parameter, and the limiting behavior at high and low temperatures is described. The noninteracting part of the free energy in this approximation is the usual Thomas-Fermi functional. These results generalize and extend to finite temperatures the ground-state regularization by R. G. Parr and S. Ghosh [Proc. Natl. Acad. Sci. U.S.A. 83, 3577 (1986), 10.1073/pnas.83.11.3577] and by L. R. Pratt, G. G. Hoffman, and R. A. Harris [J. Chem. Phys. 88, 1818 (1988), 10.1063/1.454105] and formally systematize the finite-temperature regularization given by the latter authors.

Approximate effective nonlinear coefficient of second-harmonic generation in KTiOPO(4).

PubMed

Asaumi, K

1993-10-20

A simplified approximate expression for the effective nonlinear coefficient of type-II second-harmonicgeneration in KTiOPO(4) was obtained by observing that the difference between the refractive indices n(x) and n(y) is 1 order of magnitude smaller than the difference between n(z) and n(y) (or n(x)). The agreement of this approximate equation with the true definition is good, with a maximum discrepancy of 4%.

New Hardness Results for Diophantine Approximation

NASA Astrophysics Data System (ADS)

Eisenbrand, Friedrich; Rothvoß, Thomas

We revisit simultaneous Diophantine approximation, a classical problem from the geometry of numbers which has many applications in algorithms and complexity. The input to the decision version of this problem consists of a rational vector α ∈ ℚ n , an error bound ɛ and a denominator bound N ∈ ℕ + . One has to decide whether there exists an integer, called the denominator Q with 1 ≤ Q ≤ N such that the distance of each number Q ·α i to its nearest integer is bounded by ɛ. Lagarias has shown that this problem is NP-complete and optimization versions have been shown to be hard to approximate within a factor n c/ loglogn for some constant c > 0. We strengthen the existing hardness results and show that the optimization problem of finding the smallest denominator Q ∈ ℕ + such that the distances of Q·α i to the nearest integer are bounded by ɛ is hard to approximate within a factor 2 n unless {textrm{P}} = NP.

Geometrical-optics approximation of forward scattering by gradient-index spheres.

PubMed

Li, Xiangzhen; Han, Xiang'e; Li, Renxian; Jiang, Huifen

2007-08-01

By means of geometrical optics we present an approximation method for acceleration of the computation of the scattering intensity distribution within a forward angular range (0-60 degrees ) for gradient-index spheres illuminated by a plane wave. The incident angle of reflected light is determined by the scattering angle, thus improving the approximation accuracy. The scattering angle and the optical path length are numerically integrated by a general-purpose integrator. With some special index models, the scattering angle and the optical path length can be expressed by a unique function and the calculation is faster. This method is proved effective for transparent particles with size parameters greater than 50. It fails to give good approximation results at scattering angles whose refractive rays are in the backward direction. For different index models, the geometrical-optics approximation is effective only for forward angles, typically those less than 60 degrees or when the refractive-index difference of a particle is less than a certain value.

Correlation Energies from the Two-Component Random Phase Approximation.

PubMed

Kühn, Michael

2014-02-11

The correlation energy within the two-component random phase approximation accounting for spin-orbit effects is derived. The resulting plasmon equation is rewritten-analogously to the scalar relativistic case-in terms of the trace of two Hermitian matrices for (Kramers-restricted) closed-shell systems and then represented as an integral over imaginary frequency using the resolution of the identity approximation. The final expression is implemented in the TURBOMOLE program suite. The code is applied to the computation of equilibrium distances and vibrational frequencies of heavy diatomic molecules. The efficiency is demonstrated by calculation of the relative energies of the Oh-, D4h-, and C5v-symmetric isomers of Pb6. Results within the random phase approximation are obtained based on two-component Kohn-Sham reference-state calculations, using effective-core potentials. These values are finally compared to other two-component and scalar relativistic methods, as well as experimental data.

On Born approximation in black hole scattering

NASA Astrophysics Data System (ADS)

Batic, D.; Kelkar, N. G.; Nowakowski, M.

2011-12-01

A massless field propagating on spherically symmetric black hole metrics such as the Schwarzschild, Reissner-Nordström and Reissner-Nordström-de Sitter backgrounds is considered. In particular, explicit formulae in terms of transcendental functions for the scattering of massless scalar particles off black holes are derived within a Born approximation. It is shown that the conditions on the existence of the Born integral forbid a straightforward extraction of the quasi normal modes using the Born approximation for the scattering amplitude. Such a method has been used in literature. We suggest a novel, well defined method, to extract the large imaginary part of quasinormal modes via the Coulomb-like phase shift. Furthermore, we compare the numerically evaluated exact scattering amplitude with the Born one to find that the approximation is not very useful for the scattering of massless scalar, electromagnetic as well as gravitational waves from black holes.

New approximate orientation averaging of the water molecule interacting with the thermal neutron

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

Markovic, M.I.; Minic, D.M.; Rakic, A.D.

1992-02-01

This paper reports that exactly describing the time of thermal neutron collisions with water molecules, orientation averaging is performed by an exact method (EOA{sub k}) and four approximate methods (two well known and two less known). Expressions for the microscopic scattering kernel are developed. The two well-known approximate orientation averaging methods are Krieger-Nelkin (K-N) and Koppel-Young (K-Y). The results obtained by one of the two proposed approximate orientation averaging methods agree best with the corresponding results obtained by EOA{sub k}. The largest discrepancies between the EOA{sub k} results and the results of the approximate methods are obtained using the well-knowmore » K-N approximate orientation averaging method.« less

Inverse eigenproblem for R-symmetric matrices and their approximation

NASA Astrophysics Data System (ADS)

Yuan, Yongxin

2009-11-01

Let be a nontrivial involution, i.e., R=R-1[not equal to]±In. We say that is R-symmetric if RGR=G. The set of all -symmetric matrices is denoted by . In this paper, we first give the solvability condition for the following inverse eigenproblem (IEP): given a set of vectors in and a set of complex numbers , find a matrix such that and are, respectively, the eigenvalues and eigenvectors of A. We then consider the following approximation problem: Given an n×n matrix , find such that , where is the solution set of IEP and ||[dot operator]|| is the Frobenius norm. We provide an explicit formula for the best approximation solution by means of the canonical correlation decomposition.

Flexible Approximation Model Approach for Bi-Level Integrated System Synthesis

NASA Technical Reports Server (NTRS)

Sobieszczanski-Sobieski, Jaroslaw; Kim, Hongman; Ragon, Scott; Soremekun, Grant; Malone, Brett

2004-01-01

Bi-Level Integrated System Synthesis (BLISS) is an approach that allows design problems to be naturally decomposed into a set of subsystem optimizations and a single system optimization. In the BLISS approach, approximate mathematical models are used to transfer information from the subsystem optimizations to the system optimization. Accurate approximation models are therefore critical to the success of the BLISS procedure. In this paper, new capabilities that are being developed to generate accurate approximation models for BLISS procedure will be described. The benefits of using flexible approximation models such as Kriging will be demonstrated in terms of convergence characteristics and computational cost. An approach of dealing with cases where subsystem optimization cannot find a feasible design will be investigated by using the new flexible approximation models for the violated local constraints.

Approximated analytical solution to an Ebola optimal control problem

NASA Astrophysics Data System (ADS)

Hincapié-Palacio, Doracelly; Ospina, Juan; Torres, Delfim F. M.

2016-11-01

An analytical expression for the optimal control of an Ebola problem is obtained. The analytical solution is found as a first-order approximation to the Pontryagin Maximum Principle via the Euler-Lagrange equation. An implementation of the method is given using the computer algebra system Maple. Our analytical solutions confirm the results recently reported in the literature using numerical methods.

Mean-field approximation for spacing distribution functions in classical systems

NASA Astrophysics Data System (ADS)

González, Diego Luis; Pimpinelli, Alberto; Einstein, T. L.

2012-01-01

We propose a mean-field method to calculate approximately the spacing distribution functions p(n)(s) in one-dimensional classical many-particle systems. We compare our method with two other commonly used methods, the independent interval approximation and the extended Wigner surmise. In our mean-field approach, p(n)(s) is calculated from a set of Langevin equations, which are decoupled by using a mean-field approximation. We find that in spite of its simplicity, the mean-field approximation provides good results in several systems. We offer many examples illustrating that the three previously mentioned methods give a reasonable description of the statistical behavior of the system. The physical interpretation of each method is also discussed.

Interaction and charge transfer between dielectric spheres: Exact and approximate analytical solutions.

PubMed

Lindén, Fredrik; Cederquist, Henrik; Zettergren, Henning

2016-11-21

We present exact analytical solutions for charge transfer reactions between two arbitrarily charged hard dielectric spheres. These solutions, and the corresponding exact ones for sphere-sphere interaction energies, include sums that describe polarization effects to infinite orders in the inverse of the distance between the sphere centers. In addition, we show that these exact solutions may be approximated by much simpler analytical expressions that are useful for many practical applications. This is exemplified through calculations of Langevin type cross sections for forming a compound system of two colliding spheres and through calculations of electron transfer cross sections. We find that it is important to account for dielectric properties and finite sphere sizes in such calculations, which for example may be useful for describing the evolution, growth, and dynamics of nanometer sized dielectric objects such as molecular clusters or dust grains in different environments including astrophysical ones.

Approximate scaling properties of RNA free energy landscapes

NASA Technical Reports Server (NTRS)

Baskaran, S.; Stadler, P. F.; Schuster, P.

1996-01-01

RNA free energy landscapes are analysed by means of "time-series" that are obtained from random walks restricted to excursion sets. The power spectra, the scaling of the jump size distribution, and the scaling of the curve length measured with different yard stick lengths are used to describe the structure of these "time series". Although they are stationary by construction, we find that their local behavior is consistent with both AR(1) and self-affine processes. Random walks confined to excursion sets (i.e., with the restriction that the fitness value exceeds a certain threshold at each step) exhibit essentially the same statistics as free random walks. We find that an AR(1) time series is in general approximately self-affine on timescales up to approximately the correlation length. We present an empirical relation between the correlation parameter rho of the AR(1) model and the exponents characterizing self-affinity.

Multimodal far-field acoustic radiation pattern: An approximate equation

NASA Technical Reports Server (NTRS)

Rice, E. J.

1977-01-01

The far-field sound radiation theory for a circular duct was studied for both single mode and multimodal inputs. The investigation was intended to develop a method to determine the acoustic power produced by turbofans as a function of mode cut-off ratio. With reasonable simplifying assumptions the single mode radiation pattern was shown to be reducible to a function of mode cut-off ratio only. With modal cut-off ratio as the dominant variable, multimodal radiation patterns can be reduced to a simple explicit expression. This approximate expression provides excellent agreement with an exact calculation of the sound radiation pattern using equal acoustic power per mode.

Accuracy of the adiabatic-impulse approximation for closed and open quantum systems

NASA Astrophysics Data System (ADS)

Tomka, Michael; Campos Venuti, Lorenzo; Zanardi, Paolo

2018-03-01

We study the adiabatic-impulse approximation (AIA) as a tool to approximate the time evolution of quantum states when driven through a region of small gap. Such small-gap regions are a common situation in adiabatic quantum computing and having reliable approximations is important in this context. The AIA originates from the Kibble-Zurek theory applied to continuous quantum phase transitions. The Kibble-Zurek mechanism was developed to predict the power-law scaling of the defect density across a continuous quantum phase transition. Instead, here we quantify the accuracy of the AIA via the trace norm distance with respect to the exact evolved state. As expected, we find that for short times or fast protocols, the AIA outperforms the simple adiabatic approximation. However, for large times or slow protocols, the situation is actually reversed and the AIA provides a worse approximation. Nevertheless, we found a variation of the AIA that can perform better than the adiabatic one. This counterintuitive modification consists in crossing the region of small gap twice. Our findings are illustrated by several examples of driven closed and open quantum systems.

Is Approximate Number Precision a Stable Predictor of Math Ability?

ERIC Educational Resources Information Center

Libertus, Melissa E.; Feigenson, Lisa; Halberda, Justin

2013-01-01

Previous research shows that children's ability to estimate numbers of items using their Approximate Number System (ANS) predicts later math ability. To more closely examine the predictive role of early ANS acuity on later abilities, we assessed the ANS acuity, math ability, and expressive vocabulary of preschoolers twice, six months apart. We…

A New Closed Form Approximation for BER for Optical Wireless Systems in Weak Atmospheric Turbulence

NASA Astrophysics Data System (ADS)

Kaushik, Rahul; Khandelwal, Vineet; Jain, R. C.

2018-04-01

Weak atmospheric turbulence condition in an optical wireless communication (OWC) is captured by log-normal distribution. The analytical evaluation of average bit error rate (BER) of an OWC system under weak turbulence is intractable as it involves the statistical averaging of Gaussian Q-function over log-normal distribution. In this paper, a simple closed form approximation for BER of OWC system under weak turbulence is given. Computation of BER for various modulation schemes is carried out using proposed expression. The results obtained using proposed expression compare favorably with those obtained using Gauss-Hermite quadrature approximation and Monte Carlo Simulations.

Development of a consensus approach for return of pathology incidental findings in the Genotype-Tissue Expression (GTEx) project.

PubMed

Lockhart, Nicole C; Weil, Carol J; Carithers, Latarsha J; Koester, Susan E; Little, A Roger; Volpi, Simona; Moore, Helen M; Berkman, Benjamin E

2018-06-14

The active debate about the return of incidental or secondary findings in research has primarily focused on return to research participants, or in some cases, family members. Particular attention has been paid to return of genomic findings. Yet, research may generate other types of findings that warrant consideration for return, including findings related to the pathology of donated biospecimens. In the case of deceased biospecimen donors who are also organ and/or tissue transplant donors, pathology incidental findings may be relevant not to family members, but to potential organ or tissue transplant recipients. This paper will describe the ethical implications of pathology incidental findings in the Genotype-Tissue Expression (GTEx) project, the process for developing a consensus approach as to if/when such findings should be returned, possible implications for other research projects collecting postmortem tissues and how the scenario encountered in GTEx fits into the larger return of results/incidental findings debate. © Article author(s) (or their employer(s) unless otherwise stated in the text of the article) 2018. All rights reserved. No commercial use is permitted unless otherwise expressly granted.

Low rank approximation in G 0W 0 calculations

DOE PAGES

Shao, MeiYue; Lin, Lin; Yang, Chao; ...

2016-06-04

The single particle energies obtained in a Kohn-Sham density functional theory (DFT) calculation are generally known to be poor approximations to electron excitation energies that are measured in tr ansport, tunneling and spectroscopic experiments such as photo-emission spectroscopy. The correction to these energies can be obtained from the poles of a single particle Green’s function derived from a many-body perturbation theory. From a computational perspective, the accuracy and efficiency of such an approach depends on how a self energy term that properly accounts for dynamic screening of electrons is approximated. The G 0W 0 approximation is a widely used techniquemore » in which the self energy is expressed as the convolution of a noninteracting Green’s function (G 0) and a screened Coulomb interaction (W 0) in the frequency domain. The computational cost associated with such a convolution is high due to the high complexity of evaluating W 0 at multiple frequencies. In this paper, we discuss how the cost of G 0W 0 calculation can be reduced by constructing a low rank approximation to the frequency dependent part of W 0 . In particular, we examine the effect of such a low rank approximation on the accuracy of the G 0W 0 approximation. We also discuss how the numerical convolution of G 0 and W 0 can be evaluated efficiently and accurately by using a contour deformation technique with an appropriate choice of the contour.« less

Harmonic-phase path-integral approximation of thermal quantum correlation functions

NASA Astrophysics Data System (ADS)

Robertson, Christopher; Habershon, Scott

2018-03-01

We present an approximation to the thermal symmetric form of the quantum time-correlation function in the standard position path-integral representation. By transforming to a sum-and-difference position representation and then Taylor-expanding the potential energy surface of the system to second order, the resulting expression provides a harmonic weighting function that approximately recovers the contribution of the phase to the time-correlation function. This method is readily implemented in a Monte Carlo sampling scheme and provides exact results for harmonic potentials (for both linear and non-linear operators) and near-quantitative results for anharmonic systems for low temperatures and times that are likely to be relevant to condensed phase experiments. This article focuses on one-dimensional examples to provide insights into convergence and sampling properties, and we also discuss how this approximation method may be extended to many-dimensional systems.

Mean-field approximation for spacing distribution functions in classical systems.

PubMed

González, Diego Luis; Pimpinelli, Alberto; Einstein, T L

2012-01-01

We propose a mean-field method to calculate approximately the spacing distribution functions p((n))(s) in one-dimensional classical many-particle systems. We compare our method with two other commonly used methods, the independent interval approximation and the extended Wigner surmise. In our mean-field approach, p((n))(s) is calculated from a set of Langevin equations, which are decoupled by using a mean-field approximation. We find that in spite of its simplicity, the mean-field approximation provides good results in several systems. We offer many examples illustrating that the three previously mentioned methods give a reasonable description of the statistical behavior of the system. The physical interpretation of each method is also discussed. © 2012 American Physical Society

Approximate Dynamic Programming: Combining Regional and Local State Following Approximations.

PubMed

Deptula, Patryk; Rosenfeld, Joel A; Kamalapurkar, Rushikesh; Dixon, Warren E

2018-06-01

An infinite-horizon optimal regulation problem for a control-affine deterministic system is solved online using a local state following (StaF) kernel and a regional model-based reinforcement learning (R-MBRL) method to approximate the value function. Unlike traditional methods such as R-MBRL that aim to approximate the value function over a large compact set, the StaF kernel approach aims to approximate the value function in a local neighborhood of the state that travels within a compact set. In this paper, the value function is approximated using a state-dependent convex combination of the StaF-based and the R-MBRL-based approximations. As the state enters a neighborhood containing the origin, the value function transitions from being approximated by the StaF approach to the R-MBRL approach. Semiglobal uniformly ultimately bounded (SGUUB) convergence of the system states to the origin is established using a Lyapunov-based analysis. Simulation results are provided for two, three, six, and ten-state dynamical systems to demonstrate the scalability and performance of the developed method.

Correlation energy functional within the GW -RPA: Exact forms, approximate forms, and challenges

NASA Astrophysics Data System (ADS)

Ismail-Beigi, Sohrab

2010-05-01

In principle, the Luttinger-Ward Green’s-function formalism allows one to compute simultaneously the total energy and the quasiparticle band structure of a many-body electronic system from first principles. We present approximate and exact expressions for the correlation energy within the GW -random-phase approximation that are more amenable to computation and allow for developing efficient approximations to the self-energy operator and correlation energy. The exact form is a sum over differences between plasmon and interband energies. The approximate forms are based on summing over screened interband transitions. We also demonstrate that blind extremization of such functionals leads to unphysical results: imposing physical constraints on the allowed solutions (Green’s functions) is necessary. Finally, we present some relevant numerical results for atomic systems.

Direct application of Padé approximant for solving nonlinear differential equations.

PubMed

Vazquez-Leal, Hector; Benhammouda, Brahim; Filobello-Nino, Uriel; Sarmiento-Reyes, Arturo; Jimenez-Fernandez, Victor Manuel; Garcia-Gervacio, Jose Luis; Huerta-Chua, Jesus; Morales-Mendoza, Luis Javier; Gonzalez-Lee, Mario

2014-01-01

This work presents a direct procedure to apply Padé method to find approximate solutions for nonlinear differential equations. Moreover, we present some cases study showing the strength of the method to generate highly accurate rational approximate solutions compared to other semi-analytical methods. The type of tested nonlinear equations are: a highly nonlinear boundary value problem, a differential-algebraic oscillator problem, and an asymptotic problem. The high accurate handy approximations obtained by the direct application of Padé method shows the high potential if the proposed scheme to approximate a wide variety of problems. What is more, the direct application of the Padé approximant aids to avoid the previous application of an approximative method like Taylor series method, homotopy perturbation method, Adomian Decomposition method, homotopy analysis method, variational iteration method, among others, as tools to obtain a power series solutions to post-treat with the Padé approximant. 34L30.

Effective hepatocyte transplantation using rat hepatocytes with low asialoglycoprotein receptor expression.

PubMed

Ise, Hirohiko; Nikaido, Toshio; Negishi, Naoki; Sugihara, Nobuhiro; Suzuki, Fumitaka; Akaike, Toshihiro; Ikeda, Uichi

2004-08-01

Development of a reliable method of isolating highly proliferative potential hepatocytes provides information crucial to progress in the field of hepatocyte transplantation. The aim of this study was to develop reliable hepatocyte transplantation using highly proliferative, eg, progenitor-like hepatocytes, based on asialoglycoprotein receptor (ASGPR) expression levels for hepatocyte transplantation. We have previously reported that mouse hepatocytes with low ASGPR expression levels have highly proliferative potential and can be used as progenitor-like hepatocytes. We therefore fractionated F344 male rat hepatocytes expressing low and high levels of ASGPR and determined the liver repopulation capacity of hepatocytes according to low and high ASGPR expression in the liver. Next, 2 x 10(5) cells of each type were transplanted into female liver regenerative model dipeptidyl peptidase-deficient rats, and we estimated the rate of liver repopulation by the transplanted hepatocytes in the host liver, as determined by recognition of the Sry gene on the Y-chromosome. At 60 days after hepatocyte transplantation, the transplanted hepatocytes occupied approximately 76% of the total hepatocyte mass in the case of the transplantation of hepatocytes with low ASGPR expression, but accounted for approximately 12% and 17% of the mass in the case of the transplantation of hepatocytes with high ASGPR expression and unfractionated hepatocytes, respectively. In conclusion, these findings suggest that hepatocytes with low ASGPR expression can result in normal liver function and a high repopulation capacity in vivo. These results provide insight into development of a strategy for effective liver repopulation using transplanted hepatocytes.

Countably QC-Approximating Posets

PubMed Central

Mao, Xuxin; Xu, Luoshan

2014-01-01

As a generalization of countably C-approximating posets, the concept of countably QC-approximating posets is introduced. With the countably QC-approximating property, some characterizations of generalized completely distributive lattices and generalized countably approximating posets are given. The main results are as follows: (1) a complete lattice is generalized completely distributive if and only if it is countably QC-approximating and weakly generalized countably approximating; (2) a poset L having countably directed joins is generalized countably approximating if and only if the lattice σ c(L)op of all σ-Scott-closed subsets of L is weakly generalized countably approximating. PMID:25165730

Approximations for column effect in airplane wing spars

NASA Technical Reports Server (NTRS)

Warner, Edward P; Short, Mac

1927-01-01

The significance attaching to "column effect" in airplane wing spars has been increasingly realized with the passage of time, but exact computations of the corrections to bending moment curves resulting from the existence of end loads are frequently omitted because of the additional labor involved in an analysis by rigorously correct methods. The present report represents an attempt to provide for approximate column effect corrections that can be graphically or otherwise expressed so as to be applied with a minimum of labor. Curves are plotted giving approximate values of the correction factors for single and two bay trusses of varying proportions and with various relationships between axial and lateral loads. It is further shown from an analysis of those curves that rough but useful approximations can be obtained from Perry's formula for corrected bending moment, with the assumed distance between points of inflection arbitrarily modified in accordance with rules given in the report. The discussion of general rules of variation of bending stress with axial load is accompanied by a study of the best distribution of the points of support along a spar for various conditions of loading.

Exponential Approximations Using Fourier Series Partial Sums

NASA Technical Reports Server (NTRS)

Banerjee, Nana S.; Geer, James F.

1997-01-01

The problem of accurately reconstructing a piece-wise smooth, 2(pi)-periodic function f and its first few derivatives, given only a truncated Fourier series representation of f, is studied and solved. The reconstruction process is divided into two steps. In the first step, the first 2N + 1 Fourier coefficients of f are used to approximate the locations and magnitudes of the discontinuities in f and its first M derivatives. This is accomplished by first finding initial estimates of these quantities based on certain properties of Gibbs phenomenon, and then refining these estimates by fitting the asymptotic form of the Fourier coefficients to the given coefficients using a least-squares approach. It is conjectured that the locations of the singularities are approximated to within O(N(sup -M-2), and the associated jump of the k(sup th) derivative of f is approximated to within O(N(sup -M-l+k), as N approaches infinity, and the method is robust. These estimates are then used with a class of singular basis functions, which have certain 'built-in' singularities, to construct a new sequence of approximations to f. Each of these new approximations is the sum of a piecewise smooth function and a new Fourier series partial sum. When N is proportional to M, it is shown that these new approximations, and their derivatives, converge exponentially in the maximum norm to f, and its corresponding derivatives, except in the union of a finite number of small open intervals containing the points of singularity of f. The total measure of these intervals decreases exponentially to zero as M approaches infinity. The technique is illustrated with several examples.

Capacitor-Chain Successive-Approximation ADC

NASA Technical Reports Server (NTRS)

Cunningham, Thomas

2003-01-01

A proposed successive-approximation analog-to-digital converter (ADC) would contain a capacitively terminated chain of identical capacitor cells. Like a conventional successive-approximation ADC containing a bank of binary-scaled capacitors, the proposed ADC would store an input voltage on a sample-and-hold capacitor and would digitize the stored input voltage by finding the closest match between this voltage and a capacitively generated sum of binary fractions of a reference voltage (Vref). However, the proposed capacitor-chain ADC would offer two major advantages over a conventional binary-scaled-capacitor ADC: (1) In a conventional ADC that digitizes to n bits, the largest capacitor (representing the most significant bit) must have 2(exp n-1) times as much capacitance, and hence, approximately 2(exp n-1) times as much area as does the smallest capacitor (representing the least significant bit), so that the total capacitor area must be 2(exp n) times that of the smallest capacitor. In the proposed capacitor-chain ADC, there would be three capacitors per cell, each approximately equal to the smallest capacitor in the conventional ADC, and there would be one cell per bit. Therefore, the total capacitor area would be only about 3(exp n) times that of the smallest capacitor. The net result would be that the proposed ADC could be considerably smaller than the conventional ADC. (2) Because of edge effects, parasitic capacitances, and manufacturing tolerances, it is difficult to make capacitor banks in which the values of capacitance are scaled by powers of 2 to the required precision. In contrast, because all the capacitors in the proposed ADC would be identical, the problem of precise binary scaling would not arise.

Compensation for intracellular environment in expression levels of mammalian circadian clock genes

PubMed Central

Matsumura, Ritsuko; Okamoto, Akihiko; Node, Koichi; Akashi, Makoto

2014-01-01

The circadian clock is driven by transcriptional oscillation of clock genes in almost all body cells. To investigate the effect of cell type-specific intracellular environment on the circadian machinery, we examined gene expression profiles in five peripheral tissues. As expected, the phase relationship between expression rhythms of nine clock genes was similar in all tissues examined. We also compared relative expression levels of clock genes among tissues, and unexpectedly found that quantitative variation remained within an approximately three-fold range, which was substantially smaller than that of metabolic housekeeping genes. Interestingly, circadian gene expression was little affected even when fibroblasts were cultured with different concentrations of serum. Together, these findings support a hypothesis that expression levels of clock genes are quantitatively compensated for the intracellular environment, such as redox potential and metabolite composition. However, more comprehensive studies are required to reach definitive conclusions. PMID:24504324

Posed versus spontaneous facial expressions are modulated by opposite cerebral hemispheres.

PubMed

Ross, Elliott D; Pulusu, Vinay K

2013-05-01

Clinical research has indicated that the left face is more expressive than the right face, suggesting that modulation of facial expressions is lateralized to the right hemisphere. The findings, however, are controversial because the results explain, on average, approximately 4% of the data variance. Using high-speed videography, we sought to determine if movement-onset asymmetry was a more powerful research paradigm than terminal movement asymmetry. The results were very robust, explaining up to 70% of the data variance. Posed expressions began overwhelmingly on the right face whereas spontaneous expressions began overwhelmingly on the left face. This dichotomy was most robust for upper facial expressions. In addition, movement-onset asymmetries did not predict terminal movement asymmetries, which were not significantly lateralized. The results support recent neuroanatomic observations that upper versus lower facial movements have different forebrain motor representations and recent behavioral constructs that posed versus spontaneous facial expressions are modulated preferentially by opposite cerebral hemispheres and that spontaneous facial expressions are graded rather than non-graded movements. Published by Elsevier Ltd.

ERK1/2 mediates glucose-regulated POMC gene expression in hypothalamic neurons.

PubMed

Zhang, Juan; Zhou, Yunting; Chen, Cheng; Yu, Feiyuan; Wang, Yun; Gu, Jiang; Ma, Lian; Ho, Guyu

2015-04-01

Hypothalamic glucose-sensing neurons regulate the expression of genes encoding feeding-related neuropetides POMC, AgRP, and NPY - the key components governing metabolic homeostasis. AMP-activated protein kinase (AMPK) is postulated to be the molecular mediator relaying glucose signals to regulate the expression of these neuropeptides. Whether other signaling mediator(s) plays a role is not clear. In this study, we investigated the role of ERK1/2 using primary hypothalamic neurons as the model system. The primary neurons were differentiated from hypothalamic progenitor cells. The differentiated neurons possessed the characteristic neuronal cell morphology and expressed neuronal post-mitotic markers as well as leptin-regulated orexigenic POMC and anorexigenic AgRP/NPY genes. Treatment of cells with glucose dose-dependently increased POMC and decreased AgRP/NPY expression with a concurrent suppression of AMPK phosphorylation. In addition, glucose treatment dose-dependently increased the ERK1/2 phosphorylation. Blockade of ERK1/2 activity with its specific inhibitor PD98059 partially (approximately 50%) abolished glucose-induced POMC expression, but had little effect on AgRP/NPY expression. Conversely, blockade of AMPK activity with its specific inhibitor produced a partial (approximately 50%) reversion of low-glucose-suppressed POMC expression, but almost completely blunted the low-glucose-induced AgRP/NPY expression. The results indicate that ERK1/2 mediated POMC but not AgRP/NPY expression. Confirming the in vitro findings, i.c.v. administration of PD98059 in rats similarly attenuated glucose-induced POMC expression in the hypothalamus, but again had little effect on AgRP/NPY expression. The results are indicative of a novel role of ERK1/2 in glucose-regulated POMC expression and offer new mechanistic insights into hypothalamic glucose sensing. © 2015 Society for Endocrinology.

A Space-Saving Approximation Algorithm for Grammar-Based Compression

NASA Astrophysics Data System (ADS)

Sakamoto, Hiroshi; Maruyama, Shirou; Kida, Takuya; Shimozono, Shinichi

A space-efficient approximation algorithm for the grammar-based compression problem, which requests for a given string to find a smallest context-free grammar deriving the string, is presented. For the input length n and an optimum CFG size g, the algorithm consumes only O(g log g) space and O(n log*n) time to achieve O((log*n)log n) approximation ratio to the optimum compression, where log*n is the maximum number of logarithms satisfying log log…log n > 1. This ratio is thus regarded to almost O(log n), which is the currently best approximation ratio. While g depends on the string, it is known that g =Ω(log n) and g=\\\\Omega(\\\\log n) and g=O\\\\left(\\\\frac{n}{log_kn}\\\\right) for strings from k-letter alphabet[12].

Phonons in random alloys: The itinerant coherent-potential approximation

NASA Astrophysics Data System (ADS)

Ghosh, Subhradip; Leath, P. L.; Cohen, Morrel H.

2002-12-01

We present the itinerant coherent-potential approximation (ICPA), an analytic, translationally invariant, and tractable form of augmented-space-based multiple-scattering theory18 in a single-site approximation for harmonic phonons in realistic random binary alloys with mass and force-constant disorder. We provide expressions for quantities needed for comparison with experimental structure factors such as partial and average spectral functions and derive the sum rules associated with them. Numerical results are presented for Ni55Pd45 and Ni50Pt50 alloys which serve as test cases, the former for weak force-constant disorder and the latter for strong. We present results on dispersion curves and disorder-induced widths. Direct comparisons with the single-site coherent potential approximation (CPA) and experiment are made which provide insight into the physics of force-constant changes in random alloys. The CPA accounts well for the weak force-constant disorder case but fails for strong force-constant disorder where the ICPA succeeds.

The TSP-approach to approximate solving the m-Cycles Cover Problem

NASA Astrophysics Data System (ADS)

Gimadi, Edward Kh.; Rykov, Ivan; Tsidulko, Oxana

2016-10-01

In the m-Cycles Cover problem it is required to find a collection of m vertex-disjoint cycles that covers all vertices of the graph and the total weight of edges in the cover is minimum (or maximum). The problem is a generalization of the Traveling salesmen problem. It is strongly NP-hard. We discuss a TSP-approach that gives polynomial approximate solutions for this problem. It transforms an approximation TSP algorithm into an approximation m-CCP algorithm. In this paper we present a number of successful transformations with proven performance guarantees for the obtained solutions.

Estimating proportions in petrographic mixing equations by least-squares approximation.

PubMed

Bryan, W B; Finger, L W; Chayes, F

1969-02-28

Petrogenetic hypotheses involving fractional crystallization, assimilation, or mixing of magmas may be expressed and tested as problems in leastsquares approximation. The calculation uses all of the data and yields a unique solution for each model, thus avoiding the ambiguity inherent in graphical or trial-and-error procedures. The compositional change in the 1960 lavas of Kilauea Volcano, Hawaii, is used to illustrate the method of calculation.

Approximate symmetries of Hamiltonians

NASA Astrophysics Data System (ADS)

Chubb, Christopher T.; Flammia, Steven T.

2017-08-01

We explore the relationship between approximate symmetries of a gapped Hamiltonian and the structure of its ground space. We start by considering approximate symmetry operators, defined as unitary operators whose commutators with the Hamiltonian have norms that are sufficiently small. We show that when approximate symmetry operators can be restricted to the ground space while approximately preserving certain mutual commutation relations. We generalize the Stone-von Neumann theorem to matrices that approximately satisfy the canonical (Heisenberg-Weyl-type) commutation relations and use this to show that approximate symmetry operators can certify the degeneracy of the ground space even though they only approximately form a group. Importantly, the notions of "approximate" and "small" are all independent of the dimension of the ambient Hilbert space and depend only on the degeneracy in the ground space. Our analysis additionally holds for any gapped band of sufficiently small width in the excited spectrum of the Hamiltonian, and we discuss applications of these ideas to topological quantum phases of matter and topological quantum error correcting codes. Finally, in our analysis, we also provide an exponential improvement upon bounds concerning the existence of shared approximate eigenvectors of approximately commuting operators under an added normality constraint, which may be of independent interest.

Optimal causal inference: estimating stored information and approximating causal architecture.

PubMed

Still, Susanne; Crutchfield, James P; Ellison, Christopher J

2010-09-01

We introduce an approach to inferring the causal architecture of stochastic dynamical systems that extends rate-distortion theory to use causal shielding--a natural principle of learning. We study two distinct cases of causal inference: optimal causal filtering and optimal causal estimation. Filtering corresponds to the ideal case in which the probability distribution of measurement sequences is known, giving a principled method to approximate a system's causal structure at a desired level of representation. We show that in the limit in which a model-complexity constraint is relaxed, filtering finds the exact causal architecture of a stochastic dynamical system, known as the causal-state partition. From this, one can estimate the amount of historical information the process stores. More generally, causal filtering finds a graded model-complexity hierarchy of approximations to the causal architecture. Abrupt changes in the hierarchy, as a function of approximation, capture distinct scales of structural organization. For nonideal cases with finite data, we show how the correct number of the underlying causal states can be found by optimal causal estimation. A previously derived model-complexity control term allows us to correct for the effect of statistical fluctuations in probability estimates and thereby avoid overfitting.

A series approximation model for optical light transport and output intensity field distribution in large aspect ratio cylindrical scintillation crystals

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

Tobias, Benjamin John

A series approximation has been derived for the transport of optical photons within a cylindrically symmetric light pipe and applied to the task of evaluating both the origin and angular distribution of light reaching the output plane. This analytic expression finds particular utility in first-pass photonic design applications since it may be evaluated at a very modest computational cost and is readily parameterized for relevant design constraints. It has been applied toward quantitative exploration of various scintillation crystal preparations and their impact on both quantum efficiency and noise, reproducing sensible dependencies and providing physical justification for certain gamma ray cameramore » design choices.« less

Description of high-power laser radiation in the paraxial approximation

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

Milant'ev, V P; Karnilovich, S P; Shaar, Ya N

2015-11-30

We consider the feasibility of an adequate description of a laser pulse of arbitrary shape within the framework of the paraxial approximation. In this approximation, using a parabolic equation and an expansion in the small parameter, expressions are obtained for the field of a sufficiently intense laser radiation given in the form of axially symmetric Hermite – Gaussian beams of arbitrary mode and arbitrary polarisation. It is shown that in the case of sufficiently short pulses, corrections to the transverse components of the laser field are the first-order rather than the secondorder quantities in the expansion in the small parameter.more » The peculiarities of the description of higher-mode Hermite – Gaussian beams are outlined. (light wave transformation)« less

Collisional excitation of the highly excited hydrogen atoms in the dipole form of the semiclassical impact parameter and Born approximations

NASA Technical Reports Server (NTRS)

Omidvar, K.

1971-01-01

Expressions for the excitation cross section of the highly excited states of the hydrogenlike atoms by fast charged particles have been derived in the dipole approximation of the semiclassical impact parameter and the Born approximations, making use of a formula for the asymptotic expansion of the oscillator strength of the hydrogenlike atoms given by Menzel. When only the leading term in the asymptotic expansion is retained, the expression for the cross section becomes identical to the expression obtained by the method of the classical collision and correspondence principle given by Percival and Richards. Comparisons are made between the Bethe coefficients obtained here and the Bethe coefficients of the Born approximation for transitions where the Born calculation is available. Satisfactory agreement is obtained only for n yields n + 1 transitions, with n the principal quantum number of the excited state.

Plasma Physics Approximations in Ares

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

Managan, R. A.

Lee & More derived analytic forms for the transport properties of a plasma. Many hydro-codes use their formulae for electrical and thermal conductivity. The coefficients are complex functions of Fermi-Dirac integrals, F n( μ/θ ), the chemical potential, μ or ζ = ln(1+e μ/θ ), and the temperature, θ = kT. Since these formulae are expensive to compute, rational function approximations were fit to them. Approximations are also used to find the chemical potential, either μ or ζ . The fits use ζ as the independent variable instead of μ/θ . New fits are provided for A α (ζmore » ),A β (ζ ), ζ, f(ζ ) = (1 + e -μ/θ)F 1/2(μ/θ), F 1/2'/F 1/2, F c α, and F c β. In each case the relative error of the fit is minimized since the functions can vary by many orders of magnitude. The new fits are designed to exactly preserve the limiting values in the non-degenerate and highly degenerate limits or as ζ→ 0 or ∞. The original fits due to Lee & More and George Zimmerman are presented for comparison.« less

Thermodynamic properties and static structure factor for a Yukawa fluid in the mean spherical approximation.

PubMed

Montes-Perez, J; Cruz-Vera, A; Herrera, J N

2011-12-01

This work presents the full analytic expressions for the thermodynamic properties and the static structure factor for a hard sphere plus 1-Yukawa fluid within the mean spherical approximation. To obtain these properties of the fluid type Yukawa analytically it was necessary to solve an equation of fourth order for the scaling parameter on a large scale. The physical root of this equation was determined by imposing physical conditions. The results of this work are obtained from seminal papers of Blum and Høye. We show that is not necessary the use the series expansion to solve the equation for the scaling parameter. We applied our theoretical result to find the thermodynamic and the static structure factor for krypton. Our results are in good agreement with those obtained in an experimental form or by simulation using the Monte Carlo method.

Thermally Driven One-Fluid Electron-Proton Solar Wind: Eight-Moment Approximation

NASA Astrophysics Data System (ADS)

Olsen, Espen Lyngdal; Leer, Egil

1996-05-01

In an effort to improve the "classical" solar wind model, we study an eight-moment approximation hydrodynamic solar wind model, in which the full conservation equation for the heat conductive flux is solved together with the conservation equations for mass, momentum, and energy. We consider two different cases: In one model the energy flux needed to drive the solar wind is supplied as heat flux from a hot coronal base, where both the density and temperature are specified. In the other model, the corona is heated. In that model, the coronal base density and temperature are also specified, but the temperature increases outward from the coronal base due to a specified energy flux that is dissipated in the corona. The eight-moment approximation solutions are compared with the results from a "classical" solar wind model in which the collision-dominated gas expression for the heat conductive flux is used. It is shown that the "classical" expression for the heat conductive flux is generally not valid in the solar wind. In collisionless regions of the flow, the eight-moment approximation gives a larger thermalization of the heat conductive flux than the models using the collision-dominated gas approximation for the heat flux, but the heat flux is still larger than the "saturation heat flux." This leads to a breakdown of the electron distribution function, which turns negative in the collisionless region of the flow. By increasing the interaction between the electrons, the heat flux is reduced, and a reasonable shape is obtained on the distribution function. By solving the full set of equations consistent with the eight-moment distribution function for the electrons, we are thus able to draw inferences about the validity of the eight-moment description of the solar wind as well as the validity of the very commonly used collision-dominated gas approximation for the heat conductive flux in the solar wind.

Localization and stationary phase approximation on supermanifolds

NASA Astrophysics Data System (ADS)

Zakharevich, Valentin

2017-08-01

Given an odd vector field Q on a supermanifold M and a Q-invariant density μ on M, under certain compactness conditions on Q, the value of the integral ∫Mμ is determined by the value of μ on any neighborhood of the vanishing locus N of Q. We present a formula for the integral in the case where N is a subsupermanifold which is appropriately non-degenerate with respect to Q. In the process, we discuss the linear algebra necessary to express our result in a coordinate independent way. We also extend the stationary phase approximation and the Morse-Bott lemma to supermanifolds.

Serial transmission in rodents of neurodegeneration from transgenic mice expressing mutant prion protein.

PubMed Central

Hsiao, K K; Groth, D; Scott, M; Yang, S L; Serban, H; Rapp, D; Foster, D; Torchia, M; Dearmond, S J; Prusiner, S B

1994-01-01

Two lines of transgenic (Tg) mice expressing high (H) levels of the mutant P101L prion protein (PrP) developed a neurologic illness and central nervous system pathology indistinguishable from experimental murine scrapie; these mice were designated Tg(MoPrP-P101L)H. Brain homogenates from Tg(MoPrP-P101L)H mice were inoculated intracerebrally into CD-1 Swiss mice, Syrian hamsters, and Tg196 mice, Tg mice expressing the MoPrP-P101L transgene at low levels. None of the CD-1 mice developed central nervous system dysfunction, whereas approximately 10% of hamsters and approximately 40% of the Tg196 mice manifested neurologic signs between 117 and 639 days after inoculation. Serial transmission of neurodegeneration in Tg196 mice and Syrian hamsters was initiated with brain extracts, producing incubation times of approximately 400 and approximately 75 days, respectively. Although the Tg(MoPrP-P101L)H mice appear to accumulate only low levels of infections prions in their brains, the serial transmission of disease to inoculated recipients argues that prion formation occurs de novo in the brains of these uninoculated animals. These Tg mouse studies, taken together with similar findings in humans dying of inherited prion diseases, provide additional evidence that prions lack a foreign nucleic acid. Images PMID:7916462

Dissociation between exact and approximate addition in developmental dyslexia.

PubMed

Yang, Xiujie; Meng, Xiangzhi

2016-09-01

Previous research has suggested that number sense and language are involved in number representation and calculation, in which number sense supports approximate arithmetic, and language permits exact enumeration and calculation. Meanwhile, individuals with dyslexia have a core deficit in phonological processing. Based on these findings, we thus hypothesized that children with dyslexia may exhibit exact calculation impairment while doing mental arithmetic. The reaction time and accuracy while doing exact and approximate addition with symbolic Arabic digits and non-symbolic visual arrays of dots were compared between typically developing children and children with dyslexia. Reaction time analyses did not reveal any differences across two groups of children, the accuracies, interestingly, revealed a distinction of approximation and exact addition across two groups of children. Specifically, two groups of children had no differences in approximation. Children with dyslexia, however, had significantly lower accuracy in exact addition in both symbolic and non-symbolic tasks than that of typically developing children. Moreover, linguistic performances were selectively associated with exact calculation across individuals. These results suggested that children with dyslexia have a mental arithmetic deficit specifically in the realm of exact calculation, while their approximation ability is relatively intact. Copyright © 2016 Elsevier Ltd. All rights reserved.

An asymptotically consistent approximant for the equatorial bending angle of light due to Kerr black holes

NASA Astrophysics Data System (ADS)

Barlow, Nathaniel S.; Weinstein, Steven J.; Faber, Joshua A.

2017-07-01

An accurate closed-form expression is provided to predict the bending angle of light as a function of impact parameter for equatorial orbits around Kerr black holes of arbitrary spin. This expression is constructed by assuring that the weak- and strong-deflection limits are explicitly satisfied while maintaining accuracy at intermediate values of impact parameter via the method of asymptotic approximants (Barlow et al 2017 Q. J. Mech. Appl. Math. 70 21-48). To this end, the strong deflection limit for a prograde orbit around an extremal black hole is examined, and the full non-vanishing asymptotic behavior is determined. The derived approximant may be an attractive alternative to computationally expensive elliptical integrals used in black hole simulations.

Mars Express recent findings and future plans

NASA Astrophysics Data System (ADS)

Titov, Dmitrij; Bibring, Jean-Pierre; Cardesin, Alejandro; Duxbury, Thomas; Forget, Francois; Giuranna, Marco; González-Galindo, Francisco; Holmström, Mats; Jaumann, Ralf; Määttänen, Anni; Martin, Patrick; Montmessin, Franck; Orosei, Roberto; Pätzold, Martin; Plaut, Jeffrey; MEx SGS Team

2017-04-01

Mars Express remains one of ESA's most scientifically productive missions whose publication record now exceeds 1000 papers. Characterization of geological processes on a local-to-regional scale by HRSC, OMEGA and partner experiments on NASA spacecraft has allowed constraining land-forming processes in space and time. Recent results suggest episodic geological activity as well as the presence of large bodies of liquid water in several provinces (e.g. Eridania Planum, Terra Chimeria) in the early and middle Amazonian epoch and formation of vast sedimentary plains north of the Hellas basin. Mars Express observations and experimental teams provided an essential contribution to the selection of the Mars-2020 landing sites. More than a decade-long record of atmospheric parameters such as temperature, dust loading, water vapor and ozone abundance, water ice and CO2 clouds distribution, collected by SPICAM, PFS and OMEGA spectrometers as well as subsequent modeling have provided key contributions to our understanding of the martian climate. The ASPERA-3 observations of ion escape covering a complete solar cycle have revealed important dependencies of the atmospheric erosion rate on parameters of the solar wind and EUV flux. Structure of the ionosphere sounded by the MARSIS radar and the MaRS radio science experiment was found to be significantly affected by the solar activity, the crustal magnetic field, as well as by the influx of meteorite and cometary dust. MARSIS and ASPERA-3 observations suggest that the sunlit ionosphere over the regions with strong crustal fields is denser and extends to higher altitudes as compared to the regions with no crustal anomalies. The ionospheric plasma expands to higher altitudes where it contacts with the solar wind plasma. Reconnection of solar magnetic field lines carried by the solar wind with field lines of crustal origin opens channels through which the ionospheric plasma escapes to space, producing strong and narrow cavities in the

Pediatric acute myeloid leukemia with NPM1 mutations is characterized by a gene expression profile with dysregulated HOX gene expression distinct from MLL-rearranged leukemias.

PubMed

Mullighan, C G; Kennedy, A; Zhou, X; Radtke, I; Phillips, L A; Shurtleff, S A; Downing, J R

2007-09-01

Somatic mutations in nucleophosmin (NPM1) occur in approximately 35% of adult acute myeloid leukemia (AML). To assess the frequency of NPM1 mutations in pediatric AML, we sequenced NPM1 in the diagnostic blasts from 93 pediatric AML patients. Six cases harbored NPM1 mutations, with each case lacking common cytogenetic abnormalities. To explore the phenotype of the AMLs with NPM1 mutations, gene expression profiles were obtained using Affymetrix U133A microarrays. NPM1 mutations were associated with increased expression of multiple homeobox genes including HOXA9, A10, B2, B6 and MEIS1. As dysregulated homeobox gene expression is also a feature of MLL-rearranged leukemia, the gene expression signatures of NPM1-mutated and MLL-rearranged leukemias were compared. Significant differences were identified between these leukemia subtypes including the expression of different HOX genes, with NPM1-mutated AML showing higher levels of expression of HOXB2, B3, B6 and D4. These results confirm recent reports of perturbed HOX expression in NPM1-mutated adult AML, and provide the first evidence that the NPM1-mutated signature is distinct from MLL-rearranged AML. These findings suggest that mutated NPM1 leads to dysregulated HOX expression via a different mechanism than MLL rearrangement.

An Extension of the Krieger-Li-Iafrate Approximation to the Optimized-Effective-Potential Method

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

Wilson, B.G.

1999-11-11

The Krieger-Li-Iafrate approximation can be expressed as the zeroth order result of an unstable iterative method for solving the integral equation form of the optimized-effective-potential method. By pre-conditioning the iterate a first order correction can be obtained which recovers the bulk of quantal oscillations missing in the zeroth order approximation. A comparison of calculated total energies are given with Krieger-Li-Iafrate, Local Density Functional, and Hyper-Hartree-Fock results for non-relativistic atoms and ions.

Combinatorial approximation algorithms for MAXCUT using random walks.

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

Seshadhri, Comandur; Kale, Satyen

We give the first combinatorial approximation algorithm for MaxCut that beats the trivial 0.5 factor by a constant. The main partitioning procedure is very intuitive, natural, and easily described. It essentially performs a number of random walks and aggregates the information to provide the partition. We can control the running time to get an approximation factor-running time tradeoff. We show that for any constant b > 1.5, there is an {tilde O}(n{sup b}) algorithm that outputs a (0.5 + {delta})-approximation for MaxCut, where {delta} = {delta}(b) is some positive constant. One of the components of our algorithm is a weakmore » local graph partitioning procedure that may be of independent interest. Given a starting vertex i and a conductance parameter {phi}, unless a random walk of length {ell} = O(log n) starting from i mixes rapidly (in terms of {phi} and {ell}), we can find a cut of conductance at most {phi} close to the vertex. The work done per vertex found in the cut is sublinear in n.« less

Thermal effects and sudden decay approximation in the curvaton scenario

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

Kitajima, Naoya; Takesako, Tomohiro; Yokoyama, Shuichiro

2014-10-01

We study the impact of a temperature-dependent curvaton decay rate on the primordial curvature perturbation generated in the curvaton scenario. Using the familiar sudden decay approximation, we obtain an analytical expression for the curvature perturbation after the decay of the curvaton. We then investigate numerically the evolution of the background and of the perturbations during the decay. We first show that the instantaneous transfer coefficient, related to the curvaton energy fraction at the decay, can be extended into a more general parameter, which depends on the net transfer of the curvaton energy into radiation energy or, equivalently, on the totalmore » entropy ratio after the complete curvaton decay. We then compute the curvature perturbation and compare this result with the sudden decay approximation prediction.« less

Approximate inference on planar graphs using loop calculus and belief progagation

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

Chertkov, Michael; Gomez, Vicenc; Kappen, Hilbert

We introduce novel results for approximate inference on planar graphical models using the loop calculus framework. The loop calculus (Chertkov and Chernyak, 2006b) allows to express the exact partition function Z of a graphical model as a finite sum of terms that can be evaluated once the belief propagation (BP) solution is known. In general, full summation over all correction terms is intractable. We develop an algorithm for the approach presented in Chertkov et al. (2008) which represents an efficient truncation scheme on planar graphs and a new representation of the series in terms of Pfaffians of matrices. We analyzemore » in detail both the loop series and the Pfaffian series for models with binary variables and pairwise interactions, and show that the first term of the Pfaffian series can provide very accurate approximations. The algorithm outperforms previous truncation schemes of the loop series and is competitive with other state-of-the-art methods for approximate inference.« less

Improvements in the Approximate Formulae for the Period of the Simple Pendulum

ERIC Educational Resources Information Center

Turkyilmazoglu, M.

2010-01-01

This paper is concerned with improvements in some exact formulae for the period of the simple pendulum problem. Two recently presented formulae are re-examined and refined rationally, yielding more accurate approximate periods. Based on the improved expressions here, a particular new formula is proposed for the period. It is shown that the derived…

Approximate spatial reasoning

NASA Technical Reports Server (NTRS)

Dutta, Soumitra

1988-01-01

Much of human reasoning is approximate in nature. Formal models of reasoning traditionally try to be precise and reject the fuzziness of concepts in natural use and replace them with non-fuzzy scientific explicata by a process of precisiation. As an alternate to this approach, it has been suggested that rather than regard human reasoning processes as themselves approximating to some more refined and exact logical process that can be carried out with mathematical precision, the essence and power of human reasoning is in its capability to grasp and use inexact concepts directly. This view is supported by the widespread fuzziness of simple everyday terms (e.g., near tall) and the complexity of ordinary tasks (e.g., cleaning a room). Spatial reasoning is an area where humans consistently reason approximately with demonstrably good results. Consider the case of crossing a traffic intersection. We have only an approximate idea of the locations and speeds of various obstacles (e.g., persons and vehicles), but we nevertheless manage to cross such traffic intersections without any harm. The details of our mental processes which enable us to carry out such intricate tasks in such apparently simple manner are not well understood. However, it is that we try to incorporate such approximate reasoning techniques in our computer systems. Approximate spatial reasoning is very important for intelligent mobile agents (e.g., robots), specially for those operating in uncertain or unknown or dynamic domains.

Weber's gravitational force as static weak field approximation

NASA Astrophysics Data System (ADS)

Tiandho, Yuant

2016-02-01

Weber's gravitational force (WGF) is one of gravitational model that can accommodate a non-static system because it depends not only on the distance but also on the velocity and the acceleration. Unlike Newton's law of gravitation, WGF can predict the anomalous of Mercury and gravitational bending of light near massive object very well. Then, some researchers use WGF as an alternative model of gravitation and propose a new mechanics theory namely the relational mechanics theory. However, currently we have known that the theory of general relativity which proposed by Einstein can explain gravity with very accurate. Through the static weak field approximation for the non-relativistic object, we also have known that the theory of general relativity will reduce to Newton's law of gravity. In this work, we expand the static weak field approximation that compatible with relativistic object and we obtain a force equation which correspond to WGF. Therefore, WGF is more precise than Newton's gravitational law. The static-weak gravitational field that we used is a solution of the Einstein's equation in the vacuum that satisfy the linear field approximation. The expression of WGF with ξ = 1 and satisfy the requirement of energy conservation are obtained after resolving the geodesic equation. By this result, we can conclude that WGF can be derived from the general relativity.

Differential TRPV1 and TRPV2 Channel Expression in Dental Pulp

PubMed Central

Gibbs, J.L.; Melnyk, J.L.; Basbaum, A.I.

2011-01-01

Hypersensitivity to thermal and mechanical stimuli can occur in painful pulpitis. To explore the neuro-anatomical basis of heat and mechanical sensitivity, we evaluated expression of TRPV1 (heat) and TRPV2 (heat/mechanical) channels in the cell bodies and terminal arborizations of neurons that innervate the dental pulp (DP) and periodontal tissues (PDL). We report that ~50% of trigeminal ganglion (TG) neurons retrogradely labeled from the DP express TRPV2, and this was significantly greater than the general expression of this channel in the TG (15%) and slightly more than what is expressed in the PDL by retrograde labeling (40%). The TRPV1 receptor, however, was less prevalent in neurons innervating the DP than their general expression in the TG (17% vs. 26%) and was more extensively expressed in neurons innervating the PDL (26%). Co-labeling studies showed that 70% of neurons that innervate the DP are myelinated. Approximately 1/3 of the retrogradely labeled neurons from the DP were calcitonin-gene-related-peptide-positive (peptide-expressing), but very few expressed the IB4 marker of non-peptidergic unmyelinated afferents. These findings suggest that the DP has a unique neurochemical innervation with regard to TRP receptor expression, which has significant implications for the mechanisms contributing to odontogenic pain and management strategies. PMID:21406609

Differential TRPV1 and TRPV2 channel expression in dental pulp.

PubMed

Gibbs, J L; Melnyk, J L; Basbaum, A I

2011-06-01

Hypersensitivity to thermal and mechanical stimuli can occur in painful pulpitis. To explore the neuro-anatomical basis of heat and mechanical sensitivity, we evaluated expression of TRPV1 (heat) and TRPV2 (heat/mechanical) channels in the cell bodies and terminal arborizations of neurons that innervate the dental pulp (DP) and periodontal tissues (PDL). We report that ~50% of trigeminal ganglion (TG) neurons retrogradely labeled from the DP express TRPV2, and this was significantly greater than the general expression of this channel in the TG (15%) and slightly more than what is expressed in the PDL by retrograde labeling (40%). The TRPV1 receptor, however, was less prevalent in neurons innervating the DP than their general expression in the TG (17% vs. 26%) and was more extensively expressed in neurons innervating the PDL (26%). Co-labeling studies showed that 70% of neurons that innervate the DP are myelinated. Approximately 1/3 of the retrogradely labeled neurons from the DP were calcitonin-gene-related-peptide-positive (peptide-expressing), but very few expressed the IB4 marker of non-peptidergic unmyelinated afferents. These findings suggest that the DP has a unique neurochemical innervation with regard to TRP receptor expression, which has significant implications for the mechanisms contributing to odontogenic pain and management strategies.

Maximizing capture of gene co-expression relationships through pre-clustering of input expression samples: an Arabidopsis case study.

PubMed

Feltus, F Alex; Ficklin, Stephen P; Gibson, Scott M; Smith, Melissa C

2013-06-05

In genomics, highly relevant gene interaction (co-expression) networks have been constructed by finding significant pair-wise correlations between genes in expression datasets. These networks are then mined to elucidate biological function at the polygenic level. In some cases networks may be constructed from input samples that measure gene expression under a variety of different conditions, such as for different genotypes, environments, disease states and tissues. When large sets of samples are obtained from public repositories it is often unmanageable to associate samples into condition-specific groups, and combining samples from various conditions has a negative effect on network size. A fixed significance threshold is often applied also limiting the size of the final network. Therefore, we propose pre-clustering of input expression samples to approximate condition-specific grouping of samples and individual network construction of each group as a means for dynamic significance thresholding. The net effect is increase sensitivity thus maximizing the total co-expression relationships in the final co-expression network compendium. A total of 86 Arabidopsis thaliana co-expression networks were constructed after k-means partitioning of 7,105 publicly available ATH1 Affymetrix microarray samples. We term each pre-sorted network a Gene Interaction Layer (GIL). Random Matrix Theory (RMT), an un-supervised thresholding method, was used to threshold each of the 86 networks independently, effectively providing a dynamic (non-global) threshold for the network. The overall gene count across all GILs reached 19,588 genes (94.7% measured gene coverage) and 558,022 unique co-expression relationships. In comparison, network construction without pre-sorting of input samples yielded only 3,297 genes (15.9%) and 129,134 relationships. in the global network. Here we show that pre-clustering of microarray samples helps approximate condition-specific networks and allows for dynamic

Maximizing capture of gene co-expression relationships through pre-clustering of input expression samples: an Arabidopsis case study

PubMed Central

2013-01-01

Background In genomics, highly relevant gene interaction (co-expression) networks have been constructed by finding significant pair-wise correlations between genes in expression datasets. These networks are then mined to elucidate biological function at the polygenic level. In some cases networks may be constructed from input samples that measure gene expression under a variety of different conditions, such as for different genotypes, environments, disease states and tissues. When large sets of samples are obtained from public repositories it is often unmanageable to associate samples into condition-specific groups, and combining samples from various conditions has a negative effect on network size. A fixed significance threshold is often applied also limiting the size of the final network. Therefore, we propose pre-clustering of input expression samples to approximate condition-specific grouping of samples and individual network construction of each group as a means for dynamic significance thresholding. The net effect is increase sensitivity thus maximizing the total co-expression relationships in the final co-expression network compendium. Results A total of 86 Arabidopsis thaliana co-expression networks were constructed after k-means partitioning of 7,105 publicly available ATH1 Affymetrix microarray samples. We term each pre-sorted network a Gene Interaction Layer (GIL). Random Matrix Theory (RMT), an un-supervised thresholding method, was used to threshold each of the 86 networks independently, effectively providing a dynamic (non-global) threshold for the network. The overall gene count across all GILs reached 19,588 genes (94.7% measured gene coverage) and 558,022 unique co-expression relationships. In comparison, network construction without pre-sorting of input samples yielded only 3,297 genes (15.9%) and 129,134 relationships. in the global network. Conclusions Here we show that pre-clustering of microarray samples helps approximate condition

Quantum cognition based on an ambiguous representation derived from a rough set approximation.

PubMed

Gunji, Yukio-Pegio; Sonoda, Kohei; Basios, Vasileios

2016-03-01

Over the last years, in a series papers by Arecchi and others, a model for the cognitive processes involved in decision making has been proposed and investigated. The key element of this model is the expression of apprehension and judgment, basic cognitive process of decision making, as an inverse Bayes inference classifying the information content of neuron spike trains. It has been shown that for successive plural stimuli this inference, equipped with basic non-algorithmic jumps, is affected by quantum-like characteristics. We show here that such a decision making process is related consistently with an ambiguous representation by an observer within a universe of discourse. In our work the ambiguous representation of an object or a stimuli is defined as a pair of maps from objects of a set to their representations, where these two maps are interrelated in a particular structure. The a priori and a posteriori hypotheses in Bayes inference are replaced by the upper and lower approximations, correspondingly, for the initial data sets that are derived with respect to each map. Upper and lower approximations herein are defined in the context of "rough set" analysis. The inverse Bayes inference is implemented by the lower approximations with respect to the one map and for the upper approximation with respect to the other map for a given data set. We show further that, due to the particular structural relation between the two maps, the logical structure of such combined approximations can only be expressed as an orthomodular lattice and therefore can be represented by a quantum rather than a Boolean logic. To our knowledge, this is the first investigation aiming to reveal the concrete logic structure of inverse Bayes inference in cognitive processes. Copyright © 2016. Published by Elsevier Ireland Ltd.

Strong washout approximation to resonant leptogenesis

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

Garbrecht, Björn; Gautier, Florian; Klaric, Juraj, E-mail: garbrecht@tum.de, E-mail: florian.gautier@tum.de, E-mail: juraj.klaric@tum.de

We show that the effective decay asymmetry for resonant Leptogenesis in the strong washout regime with two sterile neutrinos and a single active flavour can in wide regions of parameter space be approximated by its late-time limit ε=Xsin(2φ)/(X{sup 2}+sin{sup 2}φ), where X=8πΔ/(|Y{sub 1}|{sup 2}+|Y{sub 2}|{sup 2}), Δ=4(M{sub 1}-M{sub 2})/(M{sub 1}+M{sub 2}), φ=arg(Y{sub 2}/Y{sub 1}), and M{sub 1,2}, Y{sub 1,2} are the masses and Yukawa couplings of the sterile neutrinos. This approximation in particular extends to parametric regions where |Y{sub 1,2}|{sup 2}>> Δ, i.e. where the width dominates the mass splitting. We generalise the formula for the effective decay asymmetry to themore » case of several flavours of active leptons and demonstrate how this quantity can be used to calculate the lepton asymmetry for phenomenological scenarios that are in agreement with the observed neutrino oscillations. We establish analytic criteria for the validity of the late-time approximation for the decay asymmetry and compare these with numerical results that are obtained by solving for the mixing and the oscillations of the sterile neutrinos. For phenomenologically viable models with two sterile neutrinos, we find that the flavoured effective late-time decay asymmetry can be applied throughout parameter space.« less

Exact and Approximate Statistical Inference for Nonlinear Regression and the Estimating Equation Approach.

PubMed

Demidenko, Eugene

2017-09-01

The exact density distribution of the nonlinear least squares estimator in the one-parameter regression model is derived in closed form and expressed through the cumulative distribution function of the standard normal variable. Several proposals to generalize this result are discussed. The exact density is extended to the estimating equation (EE) approach and the nonlinear regression with an arbitrary number of linear parameters and one intrinsically nonlinear parameter. For a very special nonlinear regression model, the derived density coincides with the distribution of the ratio of two normally distributed random variables previously obtained by Fieller (1932), unlike other approximations previously suggested by other authors. Approximations to the density of the EE estimators are discussed in the multivariate case. Numerical complications associated with the nonlinear least squares are illustrated, such as nonexistence and/or multiple solutions, as major factors contributing to poor density approximation. The nonlinear Markov-Gauss theorem is formulated based on the near exact EE density approximation.

Computational Modeling of Proteins based on Cellular Automata: A Method of HP Folding Approximation.

PubMed

Madain, Alia; Abu Dalhoum, Abdel Latif; Sleit, Azzam

2018-06-01

The design of a protein folding approximation algorithm is not straightforward even when a simplified model is used. The folding problem is a combinatorial problem, where approximation and heuristic algorithms are usually used to find near optimal folds of proteins primary structures. Approximation algorithms provide guarantees on the distance to the optimal solution. The folding approximation approach proposed here depends on two-dimensional cellular automata to fold proteins presented in a well-studied simplified model called the hydrophobic-hydrophilic model. Cellular automata are discrete computational models that rely on local rules to produce some overall global behavior. One-third and one-fourth approximation algorithms choose a subset of the hydrophobic amino acids to form H-H contacts. Those algorithms start with finding a point to fold the protein sequence into two sides where one side ignores H's at even positions and the other side ignores H's at odd positions. In addition, blocks or groups of amino acids fold the same way according to a predefined normal form. We intend to improve approximation algorithms by considering all hydrophobic amino acids and folding based on the local neighborhood instead of using normal forms. The CA does not assume a fixed folding point. The proposed approach guarantees one half approximation minus the H-H endpoints. This lower bound guaranteed applies to short sequences only. This is proved as the core and the folds of the protein will have two identical sides for all short sequences.

High-Dimensional Function Approximation With Neural Networks for Large Volumes of Data.

PubMed

Andras, Peter

2018-02-01

Approximation of high-dimensional functions is a challenge for neural networks due to the curse of dimensionality. Often the data for which the approximated function is defined resides on a low-dimensional manifold and in principle the approximation of the function over this manifold should improve the approximation performance. It has been show that projecting the data manifold into a lower dimensional space, followed by the neural network approximation of the function over this space, provides a more precise approximation of the function than the approximation of the function with neural networks in the original data space. However, if the data volume is very large, the projection into the low-dimensional space has to be based on a limited sample of the data. Here, we investigate the nature of the approximation error of neural networks trained over the projection space. We show that such neural networks should have better approximation performance than neural networks trained on high-dimensional data even if the projection is based on a relatively sparse sample of the data manifold. We also find that it is preferable to use a uniformly distributed sparse sample of the data for the purpose of the generation of the low-dimensional projection. We illustrate these results considering the practical neural network approximation of a set of functions defined on high-dimensional data including real world data as well.

ɛ-connectedness, finite approximations, shape theory and coarse graining in hyperspaces

NASA Astrophysics Data System (ADS)

Alonso-Morón, Manuel; Cuchillo-Ibanez, Eduardo; Luzón, Ana

2008-12-01

We use upper semifinite hyperspaces of compacta to describe ε-connectedness and to compute homology from finite approximations. We find a new connection between ε-connectedness and the so-called Shape Theory. We construct a geodesically complete R-tree, by means of ε-components at different resolutions, whose behavior at infinite captures the topological structure of the space of components of a given compact metric space. We also construct inverse sequences of finite spaces using internal finite approximations of compact metric spaces. These sequences can be converted into inverse sequences of polyhedra and simplicial maps by means of what we call the Alexandroff-McCord correspondence. This correspondence allows us to relate upper semifinite hyperspaces of finite approximation with the Vietoris-Rips complexes of such approximations at different resolutions. Two motivating examples are included in the introduction. We propose this procedure as a different mathematical foundation for problems on data analysis. This process is intrinsically related to the methodology of shape theory. This paper reinforces Robins’s idea of using methods from shape theory to compute homology from finite approximations.

Structural Reliability Analysis and Optimization: Use of Approximations

NASA Technical Reports Server (NTRS)

Grandhi, Ramana V.; Wang, Liping

1999-01-01

This report is intended for the demonstration of function approximation concepts and their applicability in reliability analysis and design. Particularly, approximations in the calculation of the safety index, failure probability and structural optimization (modification of design variables) are developed. With this scope in mind, extensive details on probability theory are avoided. Definitions relevant to the stated objectives have been taken from standard text books. The idea of function approximations is to minimize the repetitive use of computationally intensive calculations by replacing them with simpler closed-form equations, which could be nonlinear. Typically, the approximations provide good accuracy around the points where they are constructed, and they need to be periodically updated to extend their utility. There are approximations in calculating the failure probability of a limit state function. The first one, which is most commonly discussed, is how the limit state is approximated at the design point. Most of the time this could be a first-order Taylor series expansion, also known as the First Order Reliability Method (FORM), or a second-order Taylor series expansion (paraboloid), also known as the Second Order Reliability Method (SORM). From the computational procedure point of view, this step comes after the design point identification; however, the order of approximation for the probability of failure calculation is discussed first, and it is denoted by either FORM or SORM. The other approximation of interest is how the design point, or the most probable failure point (MPP), is identified. For iteratively finding this point, again the limit state is approximated. The accuracy and efficiency of the approximations make the search process quite practical for analysis intensive approaches such as the finite element methods; therefore, the crux of this research is to develop excellent approximations for MPP identification and also different

The evolution of voids in the adhesion approximation

NASA Astrophysics Data System (ADS)

Sahni, Varun; Sathyaprakah, B. S.; Shandarin, Sergei F.

1994-08-01

We apply the adhesion approximation to study the formation and evolution of voids in the universe. Our simulations-carried out using 1283 particles in a cubical box with side 128 Mpc-indicate that the void spectrum evolves with time and that the mean void size in the standard Cosmic Background Explorer Satellite (COBE)-normalized cold dark matter (CDM) model with H50 = 1 scals approximately as bar D(z) = bar Dzero/(1+2)1/2, where bar Dzero approximately = 10.5 Mpc. Interestingly, we find a strong correlation between the sizes of voids and the value of the primordial gravitational potential at void centers. This observation could in principle, pave the way toward reconstructing the form of the primordial potential from a knowledge of the observed void spectrum. Studying the void spectrum at different cosmological epochs, for spectra with a built in k-space cutoff we find that the number of voids in a representative volume evolves with time. The mean number of voids first increases until a maximum value is reached (indicating that the formation of cellular structure is complete), and then begins to decrease as clumps and filaments erge leading to hierarchical clustering and the subsequent elimination of small voids. The cosmological epoch characterizing the completion of cellular structure occurs when the length scale going nonlinear approaches the mean distance between peaks of the gravitaional potential. A central result of this paper is that voids can be populated by substructure such as mini-sheets and filaments, which run through voids. The number of such mini-pancakes that pass through a given void can be measured by the genus characteristic of an individual void which is an indicator of the topology of a given void in intial (Lagrangian) space. Large voids have on an average a larger measure than smaller voids indicating more substructure within larger voids relative to smaller ones. We find that the topology of individual voids is strongly epoch dependent

Multi-variate joint PDF for non-Gaussianities: exact formulation and generic approximations

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

Verde, Licia; Jimenez, Raul; Alvarez-Gaume, Luis

2013-06-01

We provide an exact expression for the multi-variate joint probability distribution function of non-Gaussian fields primordially arising from local transformations of a Gaussian field. This kind of non-Gaussianity is generated in many models of inflation. We apply our expression to the non-Gaussianity estimation from Cosmic Microwave Background maps and the halo mass function where we obtain analytical expressions. We also provide analytic approximations and their range of validity. For the Cosmic Microwave Background we give a fast way to compute the PDF which is valid up to more than 7σ for f{sub NL} values (both true and sampled) not ruledmore » out by current observations, which consists of expressing the PDF as a combination of bispectrum and trispectrum of the temperature maps. The resulting expression is valid for any kind of non-Gaussianity and is not limited to the local type. The above results may serve as the basis for a fully Bayesian analysis of the non-Gaussianity parameter.« less

Validity of the Born approximation for beyond Gaussian weak lensing observables

DOE PAGES

Petri, Andrea; Haiman, Zoltan; May, Morgan

2017-06-06

Accurate forward modeling of weak lensing (WL) observables from cosmological parameters is necessary for upcoming galaxy surveys. Because WL probes structures in the nonlinear regime, analytical forward modeling is very challenging, if not impossible. Numerical simulations of WL features rely on ray tracing through the outputs of N-body simulations, which requires knowledge of the gravitational potential and accurate solvers for light ray trajectories. A less accurate procedure, based on the Born approximation, only requires knowledge of the density field, and can be implemented more efficiently and at a lower computational cost. In this work, we use simulations to show thatmore » deviations of the Born-approximated convergence power spectrum, skewness and kurtosis from their fully ray-traced counterparts are consistent with the smallest nontrivial O(Φ 3) post-Born corrections (so-called geodesic and lens-lens terms). Our results imply a cancellation among the larger O(Φ 4) (and higher order) terms, consistent with previous analytic work. We also find that cosmological parameter bias induced by the Born-approximated power spectrum is negligible even for a LSST-like survey, once galaxy shape noise is considered. When considering higher order statistics such as the κ skewness and kurtosis, however, we find significant bias of up to 2.5σ. Using the LensTools software suite, we show that the Born approximation saves a factor of 4 in computing time with respect to the full ray tracing in reconstructing the convergence.« less

Validity of the Born approximation for beyond Gaussian weak lensing observables

NASA Astrophysics Data System (ADS)

Petri, Andrea; Haiman, Zoltán; May, Morgan

2017-06-01

Accurate forward modeling of weak lensing (WL) observables from cosmological parameters is necessary for upcoming galaxy surveys. Because WL probes structures in the nonlinear regime, analytical forward modeling is very challenging, if not impossible. Numerical simulations of WL features rely on ray tracing through the outputs of N -body simulations, which requires knowledge of the gravitational potential and accurate solvers for light ray trajectories. A less accurate procedure, based on the Born approximation, only requires knowledge of the density field, and can be implemented more efficiently and at a lower computational cost. In this work, we use simulations to show that deviations of the Born-approximated convergence power spectrum, skewness and kurtosis from their fully ray-traced counterparts are consistent with the smallest nontrivial O (Φ3) post-Born corrections (so-called geodesic and lens-lens terms). Our results imply a cancellation among the larger O (Φ4) (and higher order) terms, consistent with previous analytic work. We also find that cosmological parameter bias induced by the Born-approximated power spectrum is negligible even for a LSST-like survey, once galaxy shape noise is considered. When considering higher order statistics such as the κ skewness and kurtosis, however, we find significant bias of up to 2.5 σ . Using the LensTools software suite, we show that the Born approximation saves a factor of 4 in computing time with respect to the full ray tracing in reconstructing the convergence.

Noise in gene expression is coupled to growth rate

PubMed Central

Keren, Leeat; van Dijk, David; Weingarten-Gabbay, Shira; Davidi, Dan; Jona, Ghil; Weinberger, Adina; Milo, Ron; Segal, Eran

2015-01-01

Genetically identical cells exposed to the same environment display variability in gene expression (noise), with important consequences for the fidelity of cellular regulation and biological function. Although population average gene expression is tightly coupled to growth rate, the effects of changes in environmental conditions on expression variability are not known. Here, we measure the single-cell expression distributions of approximately 900 Saccharomyces cerevisiae promoters across four environmental conditions using flow cytometry, and find that gene expression noise is tightly coupled to the environment and is generally higher at lower growth rates. Nutrient-poor conditions, which support lower growth rates, display elevated levels of noise for most promoters, regardless of their specific expression values. We present a simple model of noise in expression that results from having an asynchronous population, with cells at different cell-cycle stages, and with different partitioning of the cells between the stages at different growth rates. This model predicts non-monotonic global changes in noise at different growth rates as well as overall higher variability in expression for cell-cycle–regulated genes in all conditions. The consistency between this model and our data, as well as with noise measurements of cells growing in a chemostat at well-defined growth rates, suggests that cell-cycle heterogeneity is a major contributor to gene expression noise. Finally, we identify gene and promoter features that play a role in gene expression noise across conditions. Our results show the existence of growth-related global changes in gene expression noise and suggest their potential phenotypic implications. PMID:26355006

On the validity of the use of a localized approximation for helical beams. I. Formal aspects

NASA Astrophysics Data System (ADS)

Gouesbet, Gérard; André Ambrosio, Leonardo

2018-03-01

The description of an electromagnetic beam for use in light scattering theories may be carried out by using an expansion over vector spherical wave functions with expansion coefficients expressed in terms of Beam Shape Coefficients (BSCs). A celebrated method to evaluate these BSCs has been the use of localized approximations (with several existing variants). We recently established that the use of any existing localized approximation is of limited validity in the case of Bessel and Mathieu beams. In the present paper, we address a warning against the use of any existing localized approximation in the case of helical beams. More specifically, we demonstrate that a procedure used to validate any existing localized approximation fails in the case of helical beams. Numerical computations in a companion paper will confirm that existing localized approximations are of limited validity in the case of helical beams.

Relaxation approximation in the theory of shear turbulence

NASA Technical Reports Server (NTRS)

Rubinstein, Robert

1995-01-01

Leslie's perturbative treatment of the direct interaction approximation for shear turbulence (Modern Developments in the Theory of Turbulence, 1972) is applied to derive a time dependent model for the Reynolds stresses. The stresses are decomposed into tensor components which satisfy coupled linear relaxation equations; the present theory therefore differs from phenomenological Reynolds stress closures in which the time derivatives of the stresses are expressed in terms of the stresses themselves. The theory accounts naturally for the time dependence of the Reynolds normal stress ratios in simple shear flow. The distortion of wavenumber space by the mean shear plays a crucial role in this theory.

Legendre-Tau approximation for functional differential equations. Part 3: Eigenvalue approximations and uniform stability

NASA Technical Reports Server (NTRS)

Ito, K.

1984-01-01

The stability and convergence properties of the Legendre-tau approximation for hereditary differential systems are analyzed. A charactristic equation is derived for the eigenvalues of the resulting approximate system. As a result of this derivation the uniform exponential stability of the solution semigroup is preserved under approximation. It is the key to obtaining the convergence of approximate solutions of the algebraic Riccati equation in trace norm.

Some remarks concerning the centrifugal term approximation

NASA Astrophysics Data System (ADS)

Ferreira, F. J. S.; Bezerra, V. B.

2017-10-01

We generalize the Pekeris approximation for the centrifugal term potential, l/(l +1 ) r2 , and use this to obtain the solutions of the radial Schrödinger equation for the arbitrary angular quantum number, l, of the Hulthén potential. We also obtain the expressions for the bound state energies corresponding to this potential and calculate their values for different states and compare with other results presented in the literature. We also consider some models of physical potentials, namely, the Eckart potential, the Poschl-Teller potentials, the Rosen-Morse potential, the Woods-Saxon potential, and the Manning-Rosen potential. Thus, following straightforward the example corresponding to the Hulthén potential, we show what the radial solutions and the energy spectra for these potentials are.

Analytical approximation for the Einstein-dilaton-Gauss-Bonnet black hole metric

NASA Astrophysics Data System (ADS)

Kokkotas, K. D.; Konoplya, R. A.; Zhidenko, A.

2017-09-01

We construct an analytical approximation for the numerical black hole metric of P. Kanti et al. [Phys. Rev. D 54, 5049 (1996), 10.1103/PhysRevD.54.5049] in the four-dimensional Einstein-dilaton-Gauss-Bonnet (EdGB) theory. The continued fraction expansion in terms of a compactified radial coordinate, used here, converges slowly when the dilaton coupling approaches its extremal values, but for a black hole far from the extremal state, the analytical formula has a maximal relative error of a fraction of one percent already within the third order of the continued fraction expansion. The suggested analytical representation of the numerical black hole metric is relatively compact and a good approximation in the whole space outside the black hole event horizon. Therefore, it can serve in the same way as an exact solution when analyzing particles' motion, perturbations, quasinormal modes, Hawking radiation, accreting disks, and many other problems in the vicinity of a black hole. In addition, we construct the approximate analytical expression for the dilaton field.

Approximate Dispersion Relations for Waves on Arbitrary Shear Flows

NASA Astrophysics Data System (ADS)

Ellingsen, S. À.; Li, Y.

2017-12-01

An approximate dispersion relation is derived and presented for linear surface waves atop a shear current whose magnitude and direction can vary arbitrarily with depth. The approximation, derived to first order of deviation from potential flow, is shown to produce good approximations at all wavelengths for a wide range of naturally occuring shear flows as well as widely used model flows. The relation reduces in many cases to a 3-D generalization of the much used approximation by Skop (1987), developed further by Kirby and Chen (1989), but is shown to be more robust, succeeding in situations where the Kirby and Chen model fails. The two approximations incur the same numerical cost and difficulty. While the Kirby and Chen approximation is excellent for a wide range of currents, the exact criteria for its applicability have not been known. We explain the apparently serendipitous success of the latter and derive proper conditions of applicability for both approximate dispersion relations. Our new model has a greater range of applicability. A second order approximation is also derived. It greatly improves accuracy, which is shown to be important in difficult cases. It has an advantage over the corresponding second-order expression proposed by Kirby and Chen that its criterion of accuracy is explicitly known, which is not currently the case for the latter to our knowledge. Our second-order term is also arguably significantly simpler to implement, and more physically transparent, than its sibling due to Kirby and Chen.