Onsager model for a variable dielectric permittivity near an interface

Using a generalisation of an Onsager type approach, we are able to predict a dielectric permittivity profile of an inhomogeneous dipolar fluid in the presence of a dielectric interface. The reaction and cavity fields are calculated semi-analytically using bispherical coordinates. An asymptotic expression for the local permittivity is derived as a function of distance from the interface.Comment: 20 pages, 4 figures, submitted to Molecular Physic

Penalized estimation in large-scale generalized linear array models

Large-scale generalized linear array models (GLAMs) can be challenging to fit. Computation and storage of its tensor product design matrix can be impossible due to time and memory constraints, and previously considered design matrix free algorithms do not scale well with the dimension of the parameter vector. A new design matrix free algorithm is proposed for computing the penalized maximum likelihood estimate for GLAMs, which, in particular, handles nondifferentiable penalty functions. The proposed algorithm is implemented and available via the R package \verb+glamlasso+. It combines several ideas -- previously considered separately -- to obtain sparse estimates while at the same time efficiently exploiting the GLAM structure. In this paper the convergence of the algorithm is treated and the performance of its implementation is investigated and compared to that of \verb+glmnet+ on simulated as well as real data. It is shown that the computation time fo

Multi-scale coarse-graining of diblock copolymer self-assembly: from monomers to ordered micelles

Starting from a microscopic lattice model, we investigate clustering, micellization and micelle ordering in semi-dilute solutions of AB diblock copolymers in a selective solvent. To bridge the gap in length scales, from monomers to ordered micellar structures, we implement a two-step coarse graining strategy, whereby the AB copolymers are mapped onto ``ultrasoft'' dumbells with monomer-averaged effective interactions between the centres of mass of the blocks. Monte Carlo simulations of this coarse-grained model yield clear-cut evidence for self-assembly into micelles with a mean aggregation number n of roughly 100 beyond a critical concentration. At a slightly higher concentration the micelles spontaneously undergo a disorder-order transition to a cubic phase. We determine the effective potential between these micelles from first principles.Comment: 4 pages, 4 figures, submitted to Phys. Rev. Lett

Impact of surface roughness on diffusion of confined fluids

Using event-driven molecular dynamics simulations, we quantify how the self diffusivity of confined hard-sphere fluids depends on the nature of the confining boundaries. We explore systems with featureless confining boundaries that treat particle-boundary collisions in different ways and also various types of physically (i.e., geometrically) rough boundaries. We show that, for moderately dense fluids, the ratio of the self diffusivity of a rough wall system to that of an appropriate smooth-wall reference system is a linear function of the reciprocal wall separation, with the slope depending on the nature of the roughness. We also discuss some simple practical ways to use this information to predict confined hard-sphere fluid behavior in different rough-wall systems

Na(V)1.5 sodium channel window currents contribute to spontaneous firing in olfactory sensory neurons

Olfactory sensory neurons (OSNs) fire spontaneously as well as in response to odor; both forms of firing are physiologically important. We studied voltage-gated Na+ channels in OSNs to assess their role in spontaneous activity. Whole cell patch-clamp recordings from OSNs demonstrated both tetrodotoxin-sensitive and tetrodotoxin-resistant components of Na+ current. RT-PCR showed mRNAs for five of the nine different Na+ channel α-subunits in olfactory tissue; only one was tetrodotoxin resistant, the so-called cardiac subtype NaV1.5. Immunohistochemical analysis indicated that NaV1.5 is present in the apical knob of OSN dendrites but not in the axon. The NaV1.5 channels in OSNs exhibited two important features: 1) a half-inactivation potential near −100 mV, well below the resting potential, and 2) a window current centered near the resting potential. The negative half-inactivation potential renders most NaV1.5 channels in OSNs inactivated at the resting potential, while the window current indicates that the minor fraction of noninactivated NaV1.5 channels have a small probability of opening spontaneously at the resting potential. When the tetrodotoxin-sensitive Na+ channels were blocked by nanomolar tetrodotoxin at the resting potential, spontaneous firing was suppressed as expected. Furthermore, selectively blocking NaV1.5 channels with Zn2+ in the absence of tetrodotoxin also suppressed spontaneous firing, indicating that NaV1.5 channels are required for spontaneous activity despite resting inactivation. We propose that window currents produced by noninactivated NaV1.5 channels are one source of the generator potentials that trigger spontaneous firing, while the upstroke and propagation of action potentials in OSNs are borne by the tetrodotoxin-sensitive Na+ channel subtypes.This work was aided by support from Boston University, the Rocky Mountain Taste and Smell Center Core for Cellular Visualization and Analysis [National Institute on Deafness and Other Communication Disorders (NIDCD) P30 DC-04657; D. Restrepo, principal investigator], and NIDCD Grants DC-04863 to V. Dionne and DC-006070 to D. Restrepo and T. E. Finger. (Boston University; P30 DC-04657 - Rocky Mountain Taste and Smell Center Core for Cellular Visualization and Analysis [National Institute on Deafness and Other Communication Disorders (NIDCD)]; DC-04863 - Rocky Mountain Taste and Smell Center Core for Cellular Visualization and Analysis [National Institute on Deafness and Other Communication Disorders (NIDCD)]; DC-006070 - Rocky Mountain Taste and Smell Center Core for Cellular Visualization and Analysis [National Institute on Deafness and Other Communication Disorders (NIDCD)])https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4122723/Accepted manuscrip

Cofinement, entropy, and single-particle dynamics of equilibrium hard-sphere mixtures

We use discontinuous molecular dynamics and grand-canonical transition-matrix Monte Carlo simulations to explore how confinement between parallel hard walls modifies the relationships between packing fraction, self-diffusivity, partial molar excess entropy, and total excess entropy for binary hard-sphere mixtures. To accomplish this, we introduce an efficient algorithm to calculate partial molar excess entropies from the transition-matrix Monte Carlo simulation data. We find that the species-dependent self-diffusivities of confined fluids are very similar to those of the bulk mixture if compared at the same, appropriately defined, packing fraction up to intermediate values, but then deviate negatively from the bulk behavior at higher packing fractions. On the other hand, the relationships between self-diffusivity and partial molar excess entropy (or total excess entropy) observed in the bulk fluid are preserved under confinement even at relatively high packing fractions and for different mixture compositions. This suggests that the partial molar excess entropy, calculable from classical density functional theories of inhomogeneous fluids, can be used to predict some of the nontrivial dynamical behaviors of fluid mixtures in confined environments.Comment: submitted to JC

Lasso and probabilistic inequalities for multivariate point processes

Due to its low computational cost, Lasso is an attractive regularization method for high-dimensional statistical settings. In this paper, we consider multivariate counting processes depending on an unknown function parameter to be estimated by linear combinations of a fixed dictionary. To select coefficients, we propose an adaptive $\ell_1$-penalization methodology, where data-driven weights of the penalty are derived from new Bernstein type inequalities for martingales. Oracle inequalities are established under assumptions on the Gram matrix of the dictionary. Nonasymptotic probabilistic results for multivariate Hawkes processes are proven, which allows us to check these assumptions by considering general dictionaries based on histograms, Fourier or wavelet bases. Motivated by problems of neuronal activity inference, we finally carry out a simulation study for multivariate Hawkes processes and compare our methodology with the adaptive Lasso procedure proposed by Zou in (J. Amer. Statist. Assoc. 101 (2006) 1418-1429). We observe an excellent behavior of our procedure. We rely on theoretical aspects for the essential question of tuning our methodology. Unlike adaptive Lasso of (J. Amer. Statist. Assoc. 101 (2006) 1418-1429), our tuning procedure is proven to be robust with respect to all the parameters of the problem, revealing its potential for concrete purposes, in particular in neuroscience.Comment: Published at http://dx.doi.org/10.3150/13-BEJ562 in the Bernoulli (http://isi.cbs.nl/bernoulli/) by the International Statistical Institute/Bernoulli Society (http://isi.cbs.nl/BS/bshome.htm