Search

Recently, a nonlinear Poisson equation has been introduced to model nonlinear and nonlocal hyperpolarization effects in electrostatic solute-solvent interaction for biomolecular solvation analysis. Due to a strong nonlinearity associated with the heterogeneous dielectric media, this Poisson model is difficult to solve numerically, particularly for large protein systems. A new pseudo-transient continuation approach is proposed in this paper to efficiently and stably solve the nonlinear Poisson equation. A Douglas type alternating direction implicit (ADI) method is developed for solving the pseudo-time dependent Poisson equation. Different approximations to the dielectric profile in heterogeneous media are considered in the standard finite difference discretization. The proposed ADI scheme is validated by considering benchmark examples with exact solutions and by solvation analysis of real biomolecules with various sizes. Numerical results are in good agreement with the theoretical prediction, experimental measurements, and those obtained from the boundary value problem approach. Since the time stability of the proposed ADI scheme can be maintained even using very large time increments, it is efficient for electrostatic analysis involving hyperpolarization effects.

Based on Newtonian dynamics, observations show that the luminous masses of astrophysical objects that are the size of a galaxy or larger are not enough to generate the measured motions which they supposedly determine. This is typically attributed to the existence of dark matter, which possesses mass but does not radiate (or absorb radiation). Alternatively, the mismatch can be explained if the underlying dynamics is not Newtonian. Within this conceptual scheme, Modified Newtonian Dynamics (MOND) is a successful theoretical paradigm. MOND is usually expressed in terms of a nonlinear Poisson equation, which is difficult to analyse for arbitrary matter distributions. We study the MONDian gravitational field generated by slightly non-spherically symmetric mass distributions based on the fact that both Newtonian and MONDian fields are conservative (which we refer to as the compatibility condition). As the non-relativistic version of MOND has two different formulations (AQUAL and QuMOND) and the compatibility condition can be expressed in two ways, there are four approaches to the problem in total. The method involves solving a suitably defined linear deformation potential, which generally depends on the choice of MOND interpolation function. However, for some specific form of the deformation potential, the solution is independent of the interpolation function., 30 pages, 2 figures

After a quick review of the Lane-Emden equation and its properties, a composite of two different polytropes is introduced and some of the consequences are explored. The results are used to build a nonlinear electromagnetism with non-singular, solitonic solutions as charged particles.

We investigate the structure of self-gravitating polytropic stellar systems. We present a method which allows to obtain approximate analytical solutions, $\psi_{n+\epsilon}({\bf x})$, of the nonlinear Poisson equation with the polytropic index $n+\epsilon$, given the solution $\psi_n({\bf x})$ with the polytropic index n, for any positive or negative $\epsilon$ such that $|\epsilon|\ll1$. Application of this method to the spherically symmetric stellar polytropes with $n\simeq5$ yields the solutions which describe spatially bound systems if n5. A heuristic approximate expressions for the radial profiles are also presented. Due to the duality between stellar and gas polytropes, our results are valid for gaseous, self-gravitating polytropic systems (e.g., molecular clouds) with the index $\gamma\simeq 6/5$. Stability of such systems and observational consequences for both stellar and gaseous systems are discussed., Comment: LaTeX, emulateapj, 5 pages, 2 PS figures. ApJ,accepted version. Other works are at http://www.cita.utoronto.ca/~medvedev/

We study the zero-mass-electrons limits in the hydrodynamic models for plasmas consisting of electrons and ions. We prove, under suitable assumptions, that the evolution of the ions is governed by the Euler equations coupled with a nonlinear Poisson equation, as the mass of electrons tends to zero.

We consider the steady-state two-dimensional motion of an inviscid incompressible fluid which obeys a nonlinear Poisson equation. By seeking solutions of a specific form, we arrive at some interesting new nonlinear vortex solutions.