Suspensions of polymer-grafted nanoparticles with added polymers-Structure and effective pair-interactions

We present the results of combined experimental and theoretical (molecular dynamics simulations and integral equation theory) studies of the structure and effective interactions of suspensions of polymer grafted nanoparticles (PGNPs) in the presence of linear polymers. Due to the absence of systematic experimental and theoretical studies of PGNPs, it is widely believed that the structure and effective interactions in such binary mixtures would be very similar to those of an analogous soft colloidal material-star polymers. In our study, polystyrene-grafted gold nanoparticles with functionality f = 70 were mixed with linear polystyrene (PS) of two different molecular weights for obtaining two PGNP: PS size ratios, xi = 0.14 and 2.76 (where, xi = M-g/M-m, M-g and M-m being the molecular weights of grafting and matrix polymers, respectively). The experimental structure factor of PGNPs could be modeled with an effective potential (Model-X), which has been found to be widely applicable for star polymers. Similarly, the structure factor of the blends with xi = 0.14 could be modeled reasonably well, while the structure of blends with xi = 2.76 could not be captured, especially for high density of added polymers. A model (Model-Y) for effective interactions between PGNPs in a melt of matrix polymers also failed to provide good agreement with the experimental data for samples with xi = 2.76 and high density of added polymers. We tentatively attribute this anomaly in modeling the structure factor of blends with xi = 2.76 to the questionable assumption of Model-X in describing the added polymers as star polymers with functionality 2, which gets manifested in both polymer-polymer and polymer-PGNP interactions especially at higher fractions of added polymers. The failure of Model-Y may be due to the neglect of possible many-body interactions among PGNPs mediated by matrix polymers when the fraction of added polymers is high. These observations point to the need for a new framework to understand not only the structural behavior of PGNPs but also possibly their dynamics and thermo-mechanical properties as well. (C) 2015 AIP Publishing LLC.

Does size matter? Comparative population genetics of two butterflies with different wingspans

The dispersal ability of a species is central to its biology, affecting other processes like local adaptation, population and community dynamics, and genetic structure. Among the intrinsic, species-specific factors that affect dispersal ability in butterflies, wingspan was recently shown to explain a high amount of variance in dispersal ability. In this study, a comparative approach was adopted to test whether a difference in wingspan translates into a difference in population genetic structure. Two closely related butterfly species from subfamily Satyrinae, family Nymphalidae, which are similar with respect to all traits that affect dispersal ability except for wingspan, were studied. Melanitis leda (wingspan 60-80 mm) and Ypthima baldus (wingspan 30-40 mm) were collected from the same areas along the Western Ghats of southern India. Amplified fragment length polymorphisms were used to test whether the species with a higher wingspan (M. leda) exhibited a more homogenous population genetic structure, as compared to a species with a shorter wingspan (Y. baldus). In all analyses, Y. baldus exhibited greater degree of population genetic structuring. This study is one of the few adopting a comparative approach to establish the relationship between traits that affect dispersal ability and population genetic structure.

A hybrid method for stochastic response analysis of a vibrating structure

Response analysis of a linear structure with uncertainties in both structural parameters and external excitation is considered here. When such an analysis is carried out using the spectral stochastic finite element method (SSFEM), often the computational cost tends to be prohibitive due to the rapid growth of the number of spectral bases with the number of random variables and the order of expansion. For instance, if the excitation contains a random frequency, or if it is a general random process, then a good approximation of these excitations using polynomial chaos expansion (PCE) involves a large number of terms, which leads to very high cost. To address this issue of high computational cost, a hybrid method is proposed in this work. In this method, first the random eigenvalue problem is solved using the weak formulation of SSFEM, which involves solving a system of deterministic nonlinear algebraic equations to estimate the PCE coefficients of the random eigenvalues and eigenvectors. Then the response is estimated using a Monte Carlo (MC) simulation, where the modal bases are sampled from the PCE of the random eigenvectors estimated in the previous step, followed by a numerical time integration. It is observed through numerical studies that this proposed method successfully reduces the computational burden compared with either a pure SSFEM of a pure MC simulation and more accurate than a perturbation method. The computational gain improves as the problem size in terms of degrees of freedom grows. It also improves as the timespan of interest reduces.

Vegetation structure and composition of tropical evergreen and deciduous forests in Uttara Kannada District, Western Ghats under different disturbance regimes

In the present paper, we present the structure and composition of tropical evergreen and deciduous forests in the Western Ghats monitored under a long-term programme involving Indian Institute of Science, Earthwatch and volunteer investigators from HSBC. Currently, there is limited evidence on the status and dynamics of tropical forests in the context of human disturbance and climate change. Observations made in this study show that the `more disturbed' evergreen and one of the deciduous plots have low species diversity compared to the less-disturbed forests. There are also variations in the size class structure in the more and `less disturbed' forests of all the locations. The variation is particularly noticeable in the DBH size class 10 - 15 cm category. When biomass stock estimates are considered, there was no significant difference between evergreen and deciduous forests. The difference in biomass stocks between `less disturbed' and `more disturbed' forests within a forest type is also low. Thus, the biomass and carbon stock has not been impacted despite the dependence of communities on the forests. Periodic and long-term monitoring of the status and dynamics of the forests is necessary in the context of potential increased human pressure and climate change. There is, therefore, a need to inform the communities of the impact of extraction and its effect on regeneration so as to motivate them to adopt what may be termed as ``adaptive resource management'', so as to sustain the flow of forest products.

Use of polydispersity index as control parameter to study melting/freezing of Lennard-Jones system: Comparison among predictions of bifurcation theory with Lindemann criterion, inherent structure analysis and Hansen-Verlet rule

Using polydispersity index as an additional order parameter we investigate freezing/melting transition of Lennard-Jones polydisperse systems (with Gaussian polydispersity in size), especially to gain insight into the origin of the terminal polydispersity. The average inherent structure (IS) energy and root mean square displacement (RMSD) of the solid before melting both exhibit quite similar polydispersity dependence including a discontinuity at solid-liquid transition point. Lindemann ratio, obtained from RMSD, is found to be dependent on temperature. At a given number density, there exists a value of polydispersity index (delta (P)) above which no crystalline solid is stable. This transition value of polydispersity(termed as transition polydispersity, delta (P) ) is found to depend strongly on temperature, a feature missed in hard sphere model systems. Additionally, for a particular temperature when number density is increased, delta (P) shifts to higher values. This temperature and number density dependent value of delta (P) saturates surprisingly to a value which is found to be nearly the same for all temperatures, known as terminal polydispersity (delta (TP)). This value (delta (TP) similar to 0.11) is in excellent agreement with the experimental value of 0.12, but differs from hard sphere transition where this limiting value is only 0.048. Terminal polydispersity (delta (TP)) thus has a quasiuniversal character. Interestingly, the bifurcation diagram obtained from non-linear integral equation theories of freezing seems to provide an explanation of the existence of unique terminal polydispersity in polydisperse systems. Global bond orientational order parameter is calculated to obtain further insights into mechanism for melting.

Structure-rheology relationship in a sheared lamellar fluid

The structure-rheology relationship in the shear alignment of a lamellar fluid is studied using a mesoscale model which provides access to the lamellar configurations and the rheology. Based on the equations and free energy functional, the complete set of dimensionless groups that characterize the system are the Reynolds number (rho gamma L-2/mu), the Schmidt number (mu/rho D), the Ericksen number (mu(gamma)/B), the interface sharpness parameter r, the ratio of the viscosities of the hydrophilic and hydrophobic parts mu(r), and the ratio of the system size and layer spacing (L/lambda). Here, rho and mu are the fluid density and average viscosity, (gamma) over dot is the applied strain rate, D is the coefficient of diffusion, B is the compression modulus, mu(r) is the maximum difference in the viscosity of the hydrophilic and hydrophobic parts divided by the average viscosity, and L is the system size in the cross-stream direction. The lattice Boltzmann method is used to solve the concentration and momentum equations for a two dimensional system of moderate size (L/lambda = 32) and for a low Reynolds number, and the other parameters are systematically varied to examine the qualitative features of the structure and viscosity evolution in different regimes. At low Schmidt numbers where mass diffusion is faster than momentum diffusion, there is fast local formation of randomly aligned domains with ``grain boundaries,'' which are rotated by the shear flow to align along the extensional axis as time increases. This configuration offers a high resistance to flow, and the layers do not align in the flow direction even after 1000 strain units, resulting in a viscosity higher than that for an aligned lamellar phase. At high Schmidt numbers where momentum diffusion is fast, the shear flow disrupts layers before they are fully formed by diffusion, and alignment takes place by the breakage and reformation of layers by shear, resulting in defects (edge dislocations) embedded in a background of nearly aligned layers. At high Ericksen number where the viscous forces are large compared to the restoring forces due to layer compression and bending, shear tends to homogenize the concentration field, and the viscosity decreases significantly. At very high Ericksen number, shear even disrupts the layering of the lamellar phase. At low Ericksen number, shear results in the formation of well aligned layers with edge dislocations. However, these edge dislocations take a long time to anneal; the relatively small misalignment due to the defects results in a large increase in viscosity due to high layer stiffness and due to shear localization, because the layers between defects get pinned and move as a plug with no shear. An increase in the viscosity contrast between the hydrophilic and hydrophobic parts does not alter the structural characteristics during alignment. However, there is a significant increase in the viscosity, due to pinning of the layers between defects, which results in a plug flow between defects and a localization of the shear to a part of the domain.