Changes in vegetation structure and composition of a lowland mire over a sixty-five-year interval.

Mires are characterized by plant communities of high conservation and societal value, which have experienced a major decline in area in many parts of the world, particularly Europe. Evidence suggests that they may be particularly vulnerable to changes in climate and nutrient addition. Although they have been the focus of extensive paleoecological research, few attempts have been made to examine the dynamics of mire vegetation during the current era of anthropogenic environmental change.To assess long-term change in the spatial structure and composition of a lowland mire community, in 2016 we resurveyed plots first surveyed in 1951. Measures of species richness and composition were compared between the two surveys, and changes in community composition were related to plant traits.Overall, mean species richness declined by 26%. The area of occupancy declined in 37% of species, which were primarily oligotrophic species typical of nutrient-poor bog communities. Conversely, occupancy increased in 21% of species, especially those that were more tolerant of higher nutrient availability. These changes were associated with variation in plant functional traits, as indicated by an increase mean Ellenberg trait values for nitrogen and mean temperature, and a decline in values for precipitation. These results suggest that eutrophication and climate change have been key drivers of floristic change on this site. Synthesis. This investigation provides a rare assessment of the dynamics of a mire community over a multi-decadal interval. Results indicate that substantial change has occurred in the composition of the community, and the distribution of species within it. The investigation provides evidence of the impact of environmental change on the composition and structure of a lowland mire community, and highlights challenges for its future conservation.

The opposing effects of isotropic and anisotropic attraction on association kinetics of proteins and colloids.

The association and dissociation of particles via specific anisotropic interactions is a fundamental process, both in biology (proteins) and in soft matter (colloidal patchy particles). The presence of alternative binding sites can lead to multiple productive states and also to non-productive "decoy" or intermediate states. Besides anisotropic interactions, particles can experience non-specific isotropic interactions. We employ single replica transition interface sampling to investigate how adding a non-productive binding site or a nonspecific isotropic interaction alters the dimerization kinetics of a generic patchy particle model. The addition of a decoy binding site reduces the association rate constant, independent of the site's position, while adding an isotropic interaction increases it due to an increased rebinding probability. Surprisingly, the association kinetics becomes non-monotonic for a tetramer complex formed by multivalent patchy particles. While seemingly identical to two-particle binding with a decoy state, the cooperativity of binding multiple particles leads to a kinetic optimum. Our results are relevant for the understanding and modeling of biochemical networks and self-assembly processes.

The role of multivalency in the association kinetics of patchy particle complexes.

Association and dissociation of particles are elementary steps in many natural and technological relevant processes. For many such processes, the presence of multiple binding sites is essential. For instance, protein complexes and regular structures such as virus shells are formed from elementary building blocks with multiple binding sites. Here we address a fundamental question concerning the role of multivalency of binding sites in the association kinetics of such complexes. Using single replica transition interface sampling simulations, we investigate the influence of the multivalency on the binding kinetics and the association mechanism of patchy particles that form polyhedral clusters. When the individual bond strength is fixed, the kinetics naturally is very dependent on the multivalency, with dissociation rate constants exponentially decreasing with the number of bonds. In contrast, we find that when the total bond energy per particle is kept constant, association and dissociation rate constants turn out rather independent of multivalency, although of course still very dependent on the total energy. The association and dissociation mechanisms, however, depend on the presence and nature of the intermediate states. For instance, pathways that visit intermediate states are less prevalent for particles with five binding sites compared to the case of particles with only three bonds. The presence of intermediate states can lead to kinetic trapping and malformed aggregates. We discuss implications for natural forming complexes such as virus shells and for the design of artificial colloidal patchy particles.

Modelling critical Casimir force induced self-assembly experiments on patchy colloidal dumbbells.

Colloidal particles suspended in a binary liquid mixture can interact via solvent mediated interactions, known as critical Casimir forces. For anisotropic colloids this interaction becomes directional, which leads to rich phase behavior. While experimental imaging and particle tracking techniques allow determination of isotropic effective potentials via Boltzmann inversion, the modeling of effective interaction in anisotropic systems is non-trivial precisely because of this directionality. Here we extract effective interaction potentials for non-spherical dumbbell particles from observed radial and angular distributions, by employing reference interaction site model (RISM) theory and direct Monte Carlo simulations. For colloidal dumbbell particles dispersed in a binary liquid mixture and interacting via induced critical Casimir forces, we determine the effective site-site potentials for a range of experimental temperatures. Using these potentials to simulate the system for strong Casimir forces, we reproduce the experimentally observed collapse, and provide a qualitative explanation for this behavior.

Rotational diffusion affects the dynamical self-assembly pathways of patchy particles.

Predicting the self-assembly kinetics of particles with anisotropic interactions, such as colloidal patchy particles or proteins with multiple binding sites, is important for the design of novel high-tech materials, as well as for understanding biological systems, e.g., viruses or regulatory networks. Often stochastic in nature, such self-assembly processes are fundamentally governed by rotational and translational diffusion. Whereas the rotational diffusion constant of particles is usually considered to be coupled to the translational diffusion via the Stokes-Einstein relation, in the past decade it has become clear that they can be independently altered by molecular crowding agents or via external fields. Because virus capsids naturally assemble in crowded environments such as the cell cytoplasm but also in aqueous solution in vitro, it is important to investigate how varying the rotational diffusion with respect to transitional diffusion alters the kinetic pathways of self-assembly. Kinetic trapping in malformed or intermediate structures often impedes a direct simulation approach of a kinetic network by dramatically slowing down the relaxation to the designed ground state. However, using recently developed path-sampling techniques, we can sample and analyze the entire self-assembly kinetic network of simple patchy particle systems. For assembly of a designed cluster of patchy particles we find that changing the rotational diffusion does not change the equilibrium constants, but significantly affects the dynamical pathways, and enhances (suppresses) the overall relaxation process and the yield of the target structure, by avoiding (encountering) frustrated states. Besides insight, this finding provides a design principle for improved control of nanoparticle self-assembly.

Ultrafast dynamics in iron tetracarbonyl olefin complexes investigated with two-dimensional vibrational spectroscopy.

The dynamics of iron tetracarbonyl olefin complexes has been investigated using two-dimensional infrared (2D-IR) spectroscopy. Cross peaks between all CO-stretching bands show that the CO-stretch modes are coupled, and from the cross-peak anisotropies we can confirm previous assignments of the absorption bands. From the pump-probe delay dependence of the diagonal peaks in the 2D-IR spectrum we obtain a correlation time of ∼3 ps for the spectral fluctuations of the CO-stretch modes. We observe a multi-exponential pump-probe delay dependence of the cross-peak intensities, with rate constants ranging from 0.1 ps(-1) to 0.6 ps(-1). To determine whether this delay dependence originates from fluxionality of the complex or from intramolecular vibrational relaxation (IVR), we modulate the free-energy barrier of fluxional rearrangement by varying the pi-backbonding capacities of the olefin ligand in two iron tetracarbonyl olefin complexes: Fe(CO)(4)(cinnamic acid) and Fe(CO)(4)(dimethyl fumarate). Since the pi-backbonding strongly influences the rate of fluxionality, comparing the dynamics in the two complexes allows us to determine to what extent the observed dynamics is caused by fluxionality. We conclude that on the time scale of our experiments (up to 100 ps) the cross-peak dynamics in the iron complexes is determined by intramolecular vibrational energy relaxation. Hence, in contrast to previously investigated irontricarbonyl and ironpentacarbonyl complexes, iron tetracarbonyl olefin complexes exhibit no fluxionality on the picosecond time scale.

Time-resolved vibrational spectroscopy of a molecular shuttle.

Time-resolved vibrational spectroscopy is used to investigate the inter-component motion of an ultraviolet-triggered two-station molecular shuttle. The operation cycle of this molecular shuttle involves several intermediate species, which are observable in the amide I and amide II regions of the mid-IR spectrum. Using ab initio calculations on specific parts of the rotaxane, and by comparing the transient spectra of the normal rotaxane with that of the N-deuterated version, we can assign the observed vibrational modes of each species occurring during the shuttling cycle in an unambiguous way. The complete time- and frequency-dependent data set is analyzed using singular value decomposition (SVD). Using a kinetic model to describe the time-dependent concentrations of the transient species, we derive the absorption spectra associated with each stage in the operation cycle of the molecular shuttle, including the recombination of the charged species.

Operation mechanism of a molecular machine revealed using time-resolved vibrational spectroscopy.

Rotaxanes comprise macrocycles that can shuttle between docking stations along an axle. We explored the nanosecond shuttling mechanism by reversing the relative binding affinities of two stations through ultraviolet-induced transient reduction. We monitored the ensuing changes in the CO-stretching bands of the two stations and the shuttling macrocycle by means of an infrared probing pulse. Because hydrogen-bond scission and formation at the initial and final stations led to well-resolved changes in the respective CO-stretch frequencies, the departure and arrival of the macrocycle could be observed separately. We found that the shuttling involves two steps: thermally driven escape from the initial station, followed by rapid motion along the track ending either at the initial or final station. By varying the track's length, we found that the rapid motion approximates a biased one-dimensional random walk. However, surprisingly, the direction of the overall motion is opposite that of the bias.