Blog

Revealing pseudorotation and ring-opening reactions in colloidal organic molecules | Nature Communications

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Nature Communications volume  12, Article number: 2810 (2021 ) Cite this article 2 Methyl Cyclopentane

Revealing pseudorotation and ring-opening reactions in colloidal organic molecules | Nature Communications

Colloids have a rich history of being used as ‘big atoms’ mimicking real atoms to study crystallization, gelation and the glass transition of condensed matter. Emulating the dynamics of molecules, however, has remained elusive. Recent advances in colloid chemistry allow patchy particles to be synthesized with accurate control over shape, functionality and coordination number. Here, we show that colloidal alkanes, specifically colloidal cyclopentane, assembled from tetrameric patchy particles by critical Casimir forces undergo the same chemical transformations as their atomic counterparts, allowing their dynamics to be studied in real time. We directly observe transitions between chair and twist conformations in colloidal cyclopentane, and we elucidate the interplay of bond bending strain and entropy in the molecular transition states and ring-opening reactions. These results open the door to investigate complex molecular kinetics and molecular reactions in the high-temperature classical limit, in which the colloidal analogue becomes a good model.

Theodore Hueckel, Glen M. Hocky & Stefano Sacanna

Theodore Hueckel, Glen M. Hocky, ... Stefano Sacanna

Mariano E. Brito, Marcelo A. Carignano & Verónica I. Marconi

Colloidal particles have been much used to study phase transitions, such as crystallization, gelation, and the glass transition of condensed matter at experimentally accessible length and time scales1. In this analogy, atoms are represented by homogeneous colloidal spheres with tunable interactions between them. This approach has been successful, and many of the lessons learnt are directly applicable to fields, such as photonics2, opto-electronics3, and bio-materials4.

Compared to these crystalline systems, molecules typically have a much richer structure: atoms are generally bonded directionally into molecules via valence. For example, in aliphatic organic compounds, in which a carbon atom has four bonds, such as methane (CH4), the carbon atom is in an sp3 hybridized state, binding ligands in a tetrahedral arrangement at ~109.5° bond angles. However, although many theoretical and simulation studies exist, achieving directional bonding in colloidal systems is challenging, and reliable experimental systems are scarce5,6,7. Recently developed patchy particles, decorated with patches of specific surface chemistry, with well-defined symmetry, allow reproducing the geometry of valence bonds8,9,10,11,12,13. While these particles mimic the rigid shape of small molecules, emulating molecular reactions and conformational dynamics would require degrees of freedom similar to that of organic molecules, through reversible, specific patch–patch interactions, tunable on the scale of kBT, the thermal energy11. Such interactions would open the door to colloidal molecular chemistry: molecules that consist of colloidal particles instead of atoms, reacting at the nano and micrometre scale, directly observable by microscopy. This in turn could unlock design paths for nanostructured materials14.

Specific, adjustable colloidal interactions arise in binary solvents close to their critical point. This so-called ‘critical Casimir force’ results from the confinement of solvent fluctuations between particle surfaces in a near-critical mixture, giving rise to an effective interaction that is precisely tunable with temperature and the adsorption preferences of the confining surfaces. For isotropic particles, these interactions have provided insight into gas–liquid and liquid–solid phase transitions15,16, as well as gelation17. Combined with ‘patchy particles’ exhibiting surface patches of specific solvent affinity, the interactions become directed, and may allow the reversible assembly of more complex, molecule-like structures, and their direct particle-scale observation.

Here, we show that patchy particles bonded via critical Casimir forces allow direct-space investigation of molecular dynamics, providing insight into thermally activated molecular transition states. We bind dimer and tetramer particles into colloidal alkanes, resembling alkane molecules, such as colloidal (cyclo)butane, butyne, cyclopentane, and cyclohexane, and investigate their dynamics directly in real space, using confocal microscopy. We find that just as molecular cyclopentane, colloidal cyclopentane exhibits envelope and twist conformations that interconvert in time, and we follow the kinetic pathway in great detail. We measure the free energy and corresponding bending energy directly from the observed molecular trajectory; this allows us to elucidate the interplay of bond strain and entropy in thermal and catalytic dissociation reactions of colloidal cyclopentane. These results open up directions for studying the dynamics of molecules, using precisely coordinated patchy particles, and for building complex nano- and micrometre-scale materials1.

We fabricate patchy particles from polystyrene (PS) and 3-(trimethoxysilyl)propyl methacrylate (TPM) spheres by colloidal fusion12. The synthesis yields particles with precisely linearly and tetrahedrally coordinated, fluorescently labelled TPM patches in a PS matrix (Fig. 1A). The particles have a diameter of d = 3.7 μm and a patch diameter dp = 0.5 μm; latter is sufficiently small to allow only a single other patch to bind. Di- and tetrapatch particles with number ratio 1: 3 are dispersed at a volume fraction of ϕ = 0.01 in a binary mixture of lutidine and water close to its critical point. We also add 1 mM MgSO4 to screen the electrostatic repulsion and enhance the lutidine adsorption of the patches, see Supplementary Note 1. Upon approaching the solvent demixing temperature Tcx = 33.8 °C, lutidine-rich fluctuations arise between the patches, binding them into a covalent bond analogue by critical Casimir interactions, as illustrated in Fig. 1B. To assemble colloidal molecules, we heat the sample to ΔT = 0.04 K below Tcx, resulting in a predicted binding energy of the patches of ~15kBT (see Supplementary Note 5)18.

A Scanning electron microscope (SEM) image of tetramer patchy particle with patches coloured yellow. Scale bar corresponds to 5 μm. B Schematic of the directed critical Casimir interaction: lutidine fluctuations confined between hydrophobic patches cause patch-to-patch attraction. C–G Examples of colloidal molecules: colloidal butane (C), 2-butyne (D), cyclopentane (E), cyclobutane (F), and methylcyclohexane (G), confocal superimposed onto bright-field microscope images highlighting patches in red (left), schematic reconstruction (middle) and chemical representation (right). In C, the characteristic tetrahedral bond angle α = 109.5° is indicated. In D, the linearly connecting dipatch particles mimic sp hybridized carbon atoms, which in the colloidal analogue link both dipatch and tetrapatch particles. H Probability of occurrence of colloidal molecules as a function of their size. Purple circles (left axis) show all molecules, while orange squares (right axis) show cyclic molecules only. Clearly, larger molecules are increasingly rare, while among the cyclic structures, five-particle compounds, such as cyclopentane, are most frequently observed. Data were gathered from 163,178 clusters in 457 frames. Source data are provided as a Source data file.

After a few hours, we observe bonded structures exhibiting the coordination and bond angles of carbon atoms in organic molecules, as shown in Fig. 1C–G (see Supplementary Note 2 for assembly details). A zigzag chain of four tetrapatch particles exhibiting the carbon backbone structure of butane is shown in Fig. 1D. A colloidal analogue of 2-butyne, consisting of two central dipatch particles with two tetrapatch particles capping both ends, is shown in Fig. 1E. Analogues of cyclic alkanes, ubiquitous in carbon chemistry, are shown in Fig. 1E–G. Rings consisting of four, five, and six tetrapatch particles are observed, which we identify as colloidal cyclobutane, cyclopentane, and cyclohexane. We note that unlike their atomic counterpart, these colloidal molecules have unsaturated patches, lacking their ‘hydrogen atoms’, which may lead to some difference in their reactivity and conformations.

We show size distributions of colloidal molecules in Fig. 1H. Clearly, smaller structures are in the majority, but a significant population of larger structures is present in the sample. The population decreases exponentially with size, as predicted for patchy particle systems19. Cyclic molecules, however, show a clear preference for a certain number of particles, reflecting their compatibility with the tetrahedral bond angle. Colloidal cyclobutane (Fig. 1F) is rarely encountered in our samples. This is not surprising considering its highly strained bond angles. In this configuration, two bonded neighbours make an angle of 90°, far from the ideal angle of 109.5°, causing high bond strain. For atomic cyclobutane, this high bond strain is known to cause the enhanced reactivity of cyclobutane compared to butane, making it much less stable than cyclopentane and cyclohexane that exhibit bond angles much closer to 109.5°. Indeed, we find that colloidal cyclopentane (Fig. 1E) is much more ubiquitous in the sample, and by far the most observed colloidal ring structure. Its bonds are much closer to the ideal 109.5° tetrahedral bond angle, compared to cyclobutane. Curiously, six-membered rings—cyclohexanes—are much less frequently observed, even though they have a lower bond angle strain, probably due to kinetic effects7.

We focus on the most frequently observed cyclic structure, colloidal cyclopentane, and image its three-dimensional bonding arrangement, using confocal microscopy. The three-dimensional reconstructions reveal three basic conformations (Fig. 2A): a planar conformation with all particles in the same plane, an ‘envelope’ conformation, with one particle sticking out, and a ‘twist’ conformation with one particle sticking out above, and an adjacent particle sticking out below the plane. Just as its colloidal counterpart, atomic cyclopentane also shows these conformations; like most cyclic molecules, the internal angles of the ring are not compatible with the 109.5° bonding angle of the tetrahedral symmetry. In planar cyclopentane, the angles are only slightly lower, at 108°; nevertheless, due to energetic considerations—the steric hindrance of hydrogen atoms and torsional strain—atomic cyclopentane is virtually always puckered out of plane20.

A Three-dimensional reconstructions of typical conformations of colloidal cyclopentane: planar conformation (top), ‘twist' (or ‘half-chair') conformation (left), and ‘envelope' (or ‘bend') conformation (right). B The distribution of inter-particle angles in cylcobutane (blue circles) and cyclopentane (red squares). The solid lines are a guide to the eye. C Representation of cyclopentane conformations in polar q–ϕ space. At each 0.05 ⋅  2π interval, a 2D representation of the ring is shown. The three paths (grey circles, orange triangles, and blue squares) show typical puckering routes of a thermally activated ring through q–ϕ space. The four labelled subsequent points along the grey circle path correspond to the snapshots in D. D Time series of three-dimensional configurations showing pseudorotation of colloidal cyclopentane. Snapshots are Δt = 12 s apart, while the typical relaxation time of a conformation is 24 s, see Supplementary Note 4. Symbols + and − indicate particles above and below the average plane, respectively, and arrows indicate particle movement towards the next time step. E Correlation function of out-of-plane displacements \(C({{\Delta }}j)=\langle {z}_{j}\cdot {z}_{j+{{\Delta }}j}\rangle /\langle {z}_{j}^{2}\rangle\) in a perspective view, illustrating average puckered configuration of the ring. Source data are provided as a Source data file.

Conversely, colloidal cyclobutane shows virtually no puckering, as shown by the significantly narrower bond angle distributions in Fig. 2B. Cyclobutane shows a much narrower bond angle distribution, in particular lacking the tail towards larger angles as observed for cyclopentane, which is a signature of its puckering configurations (see Supplementary Note 3 for a comparison of bond angle distributions of different colloidal molecules). Hence, while there is a driving force for out-of plane movement of cyclopentane, and thus a relatively wide angular distribution, cyclobutane bonds are stiffer, forcing the highly strained ring into a narrowly confined configuration (see Supplementary Movie 1).

For each particle in cyclopentane, there are two envelope and two twist conformations, as shown in Fig. 2C. In molecular cyclopentane, these conformations are thought to interconvert continuously. This process is known as pseudorotation, and has been suggested as early as the 40s, but so far has only been confirmed by indirect spectroscopic evidence21,22,23,24. Likewise, we find that in colloidal cyclopentane, envelope and twist conformations interconvert continuously. To follow the pseudorotation directly in real space, we rapidly acquire image stacks every 12 s, roughly two times faster than the typical relaxation time of a configuration, see Supplementary Note 4. Three-dimensional reconstructions reveal the pseudorotation of colloidal cyclopentane in Fig. 2D, and Supplementary Movies 1 and 2. From frame 1 to 2, a particle (left) flips from below to above the plane, while its in-plane neighbour moves down out of plane. From frame 3 to 4, a characteristic change from envelope to a neighbouring twist conformation occurs, corresponding to the upward motion of the (central) in-plane particle. Remarkably, during these transition states, the ring remains always highly puckered. To show this, we correlate the vertical displacements, zi, of nearest- and second-nearest-neighbour particles, and find that they are strongly anticorrelated, as shown in Fig. 2E: the upward movement of a particle typically results in the downward movement of the neighbouring particle, followed by an upward movement of the next-nearest neighbour. These correlations highlight the strongly correlated nature of the pseudorotation of colloidal cyclopentane, and suggest that interesting collective behaviour is to be expected in larger colloidal molecules.

To analyse the pseudorotation in detail, we determine the puckering amplitude q and phase ϕ from the out-of-plane displacements zi of the particles. Together, q and ϕ form a polar phase space describing all possible puckering conformations. Given an average plane through the ring, q is a measure of the resulting amplitude of the out-of-plane displacements, while ϕ tells us in what conformation the ring is, as indicated schematically in Fig. 2C25 (see ‘Methods’ for formal definitions of q and ϕ). A polar plot of the pseudorotation trajectory in (q, ϕ) space is shown in Fig. 2C (see Supplementary Movies 3 and 4 for animations). The coordinated up and down motion leads to pronounced changes in the puckering phase (phase change from frame 1–2 to 3–4). For example, the large phase change from 2 to 3 corresponds to a change from twist to a next-nearest envelope conformation, while that from 3 to 4 corresponds to a transition from twist to envelope, as shown by comparison of the trace in Fig. 2C with the corresponding 2D representations.

To explore the full statistics, we follow more rings, and determine distributions of the reduced puckering phase ϕr and amplitude q as shown in Fig. 3A, B. The reduced puckering phase indicates how far a particular configuration is between an envelope (ϕr = 0) and a twist (ϕr = 0.05) conformation. The flat distribution of ϕr indicates that, just like in atomic cyclopentane, there is no preference for either envelope or twist conformation, nor any conformation in between. This is different for the puckering amplitude q (Fig. 3B): planar conformations with q ~ 0 are almost never observed, while mildly puckered configurations with q ~ 0.06 are most prevalent. As shown by simulations in Supplementary Note 5, the peak position depends on the presence of gravity: without gravity, the maximum of the probability density is shifted towards larger values (qN ≈ 0.15) compared to the experimental measurement, meaning that the colloidal cyclopentane ring is more puckered in a system without gravity, as expected. Nevertheless, the presence of a finite puckering amplitude is surprising from an energetic point of view, as the flat ring (q = 0) has both the lowest bending and gravitational energy.

A Normalized frequency of the reduced puckering phase angle, ϕr from envelope (ϕr = 0) to twist conformation (ϕr = 0.05 ⋅  2π rad). The flat distribution reveals no preference for a particular conformation. B Frequency f(qN) of the normalized puckering amplitude qN: experimental data (green circles), Monte Carlo simulations of freely joint particles (blue squares), and particles interacting with an effective critical Casimir potential (purple pentagons), and expected values for atomic cyclopentane as given in ref. 23, assuming the largest reasonable envelope-to-twist energy barrier (solid orange line). Frequency is normalized to result in identical initial slopes. C, D Probability density maps of total bending energy as a function of ϕr and qN, respectively, for the simulated patchy particle rings. The logarithmic colourmap indicates the normalized frequency of observation from <0.01 (blue) to 1.0 (yellow). Red circles indicate the average bending energy as determined from the measured distribution of qN. Source data are provided as a Source data file.

However, the entropy of this state is also the lowest: there is only one way to place the particles into a flat ring; even a small amount of puckering will unlock many configurational microstates, thus increasing the entropy.

To estimate the entropic contribution, we consider a ring of five freely joint particles moving up and down independently. The corresponding probability distribution P(q) increases linearly with q (Fig. 3B), reflecting the increasing number of configurations. The experimental distribution initially follows this trend, but then peaks and diminishes. We associate this decrease with the increasing bending energy cost U(q) suppressing high-puckering amplitude configurations.

Together, entropic and bending energy contributions give rise to the free energy F(q) = U(q) − TS(q), where \(S(q)={k}_{{\rm{B}}}T{\mathrm{ln}}\,P(q)\) is the entropy associated with the ring configurations. In thermal equilibrium, puckering amplitudes follow a Boltzmann distribution \(f(q)\propto \exp (-F(q)/{k}_{{\rm{B}}}T)\) . We invert this relation to determine the average bending energy from the measured distribution f(q), using \(U(q)={k}_{{\rm{B}}}T{\mathrm{ln}}\,(P(q)/f(q))\) ) + U0, where U0 is an arbitrary reference energy, see Supplementary Note 6. The resultant U(q) indeed reveals an increasing bending energy with puckering amplitude, as shown in Fig. 3D (red circles).

These trends are confirmed in Monte Carlo (MC) simulations of five tetrapatch particles interacting with effective critical Casimir potentials, as detailed in Supplementary Note 5. Their puckering distributions are in excellent agreement with the experimental data, as shown in Fig. 3. Furthermore, we compute the total bending energy U directly from the strained bonds for each observed configuration; the resulting probability contour plot shows good agreement with the experimental data (Fig. 3D).

Interestingly, atomic cyclopentane shows a distribution peaked at much higher values of q. Contrary to colloidal cyclopentane, the molecule experiences additional H–H steric repulsions, as well as torsional strain, leading to a stronger degree of puckering. Nevertheless, the subsequent decline indicates occurrence of additional bending energy, similar to the colloidal analogue.

We fully exploit the direct observation of colloidal cyclopentane by studying ring-opening reactions. Rings occasionally break up as shown in Fig. 4A and Supplementary Movie 5; this break up is the thermal ring-opening reaction of colloidal cyclopentane to pentane (formally a pentane bi-radical). Interestingly, we find that just before breaking, rings always exhibit a high-puckering amplitude, as shown in Fig. 4B. Moreover, the bond breaking occurs at a particle strongly involved in the puckering, as shown by comparing the puckering phase just before bond breaking with the breaking point in Fig. 4B. Because a high-puckering amplitude corresponds to high bond strain, a large energy gain is made by breaking the bond. These results highlight the importance of bond-bending strain in the thermal dissociation of colloidal cyclopentane. Furthermore, the high-q values just before breaking suggest that the reactivity of cyclopentane is highest for extreme puckering amplitudes, in the far high-q tail of the distribution i (Fig. 3B). We speculate that a similar mechanism may apply to atomic cyclopentane. Indeed, thermal ring-breaking reactions typically occur in cyclopentane and methyl-cyclopentane26, which unlike higher-symmetry cyclic compounds like cyclohexane exhibit significant bond strain.

A, B Thermal ring-opening of free colloidal cyclopentane. A Three-dimensional reconstruction of a dissociating ring just before (left) and just after breaking the bond between the blue particles (right). Particle diameter is 3.7 μm. B Top: examples of ring configurations just before breaking (circles) and corresponding breaking points (radial lines) in q–ϕ space. Bottom: breaking points in schematic representation, with colours indicating braking position. Colour coding matches the q–ϕ coordinates above. C, D Catalytic ring-opening of colloidal cyclopentane on an attractive (hydrophobic) surface. C Three-dimensional reconstruction of the ring shown before (left) and after breaking the bond between the purple particles (right). The bright patches are attached to the surface. D Time sequence of the surface diffusion and eventual breaking of the bond between the purple particles of the surface-bound ring.

For atomic cyclopentane, a catalyst can greatly accelerate the ring-opening reaction by offering attractive binding sites. For example, noble metal surfaces and mesoporous materials can achieve some selectivity in the ring-opening reaction27,28,29,30, but the mechanism is complex and poorly understood. To investigate such catalytic reactions in colloidal cyclopentane, we assemble rings in the bulk and sediment them onto a hydrophobic surface to which the particle patches bind. We find that this greatly speeds up the ring dissociation: a few minutes after surface adsorption, a ring breaks and opens, as shown in Fig. 4C, D. Upon adsorption, three to four particle patches bind to the substrate (highlighted in Fig. 4C), confining the ring to the surface, and effectively freezing it in a single envelope conformation. At the same time, the ring can still move laterally and diffuse along the substrate, before it opens (Fig. 4D). The binding-induced angular strain, together with the entropically unfavourable locking of configurations, causes the ring to break easily. The hydrophobically treated substrate serving as a catalyst thus allows direct observation of the interplay of conformations, bond geometry and surface interaction in the catalytic dissociation process. This makes our colloidal model very suitable for the study of geometric effects of catalysts, in absence of any electronic or support effects, using designed colloidal crystal surfaces or templates with defined symmetry, lattice constant, and domain size. Fast confocal microscopy imaging can then give detailed insight into energies along the reaction coordinate, identifying the transition states in these catalytic processes.

The direct observation of pseudorotation and catalytic reactions of colloidal cyclopentane offers insight into transition states in a molecular dissociation reaction. As the thermal energy kBT is of the same order as the bond energy, the colloidal molecule corresponds to its high-temperature molecular counterpart. The ratio of 15 for the bond energy to kBT in our experiments would correspond to a temperature of ~640 K for atomic cyclopentane, as can be estimated from its C–C bond energy of ~80 kJ/mol. In this limit, the quantum mechanical energy spectrum of the low-frequency modes of the molecule becomes quasi-continuous, and classical behaviour is expected to emerge, suggesting that the colloidal analogue becomes a good model.

Our results highlight the importance of bond strain in the dissociation process, and suggest specific catalyst design that takes advantage of the puckering amplitude in the ring-breaking process. These results pave the way to the investigation of molecular kinetics by detailed direct observation of colloidal analogues, elucidating transition states, kinetic pathways, and correlations in molecular reactions31. The exquisite temperature control afforded by the critical Casimir interaction creates opportunities for molecular and supramolecular colloidal design, following equilibrium and non-equilibrium routes, and the investigation of dynamic assembly pathways. The demonstrated accurate binding control opens up pathways to ‘colloidal molecular chemistry’, in which bond-stretch and -bend potential-energy functions can be tuned by the experimentalist, and all the reactions, not only ring opening but also synthesis, homogeneous catalysis, and polymerization can be followed in real time, using a conventional microscope to observe the reacting ‘colloidal atoms’.

Patchy particles are made via colloidal fusion12 with a modified recipe. Briefly, negatively charged PS particles are synthesized by seeded emulsion polymerization, using potassium persulfate (>99.0%, Sigma-Aldrich) as initiator. After repeated seeded emulsion polymerization, negatively charged PS particles with a diameter of ~2.2 μm are resuspended into 12.5 mM NaCl solution. Silicone emulsion droplets with a diameter of 1.0 μm are made by hydrolysis and condensation of TPM (98%, Sigma-Aldrich) catalysed by ammonium hydroxide (28 wt.%, Sigma-Aldrich). The droplets are later fluorescently labelled with 3-aminopropyl trimethoxysilane-coupled rhodamine-B isothiocyanate. Pluronic F108 triblock copolymer (average Mn 14,600, Sigma-Aldrich) is added into the silicone emulsion to a final concentration of 0.05 wt.% to further stabilize the droplets. The F108-TPM emulsion is gently washed by two cycles of centrifugation/resuspension to remove free F108 polymer in emulsion, and later resuspended into 12.5 mM NaCl solution. Tetrahedral colloidal clusters are made by mixing the prepared negatively charged PS particles of 2.2 μm and F108-TPM droplets of 1.0 μm at a number ratio ~100:1. To screen the negative charge on both particle surfaces and thus allow for polymer bridging between the two, 12.5 mM NaCl is necessary. These clusters are then separated from excess PS singlets by centrifugation in glycerol/water mixture (~22 wt.% glycerol). The purified clusters are resuspended into aqueous solution containing 0.05 wt.% F108 copolymer and 0.5 wt.% dodecyltrimethylammonium bromide (≥98%, Sigma-Aldrich), to prevent aggregation, and facilitate core extrusion during colloidal fusion, separately. Pure tetrahydrofuran is then added to cluster suspension to reach a final concentration of 30% VTHF/Vcluster. The fusion is allowed to proceed for 3 min before being quenched by deionized water. The resulting patchy particles with liquid patches are then resuspended into deionized water, followed by its polymerization in 80 °C oven for 2 h. After polymerization, patchy particles typically undergo three centrifugation/resuspension cycles and transferred back into deionized water for further experiments. Dipatch particles are the fusion product of trimer colloidal clusters, which are assembled and fused by the same protocol described above, except the TPM emulsion droplets used to make clusters are smaller (diameter being ~0.6 μm).

The patchy particle dimensions were determined using atomic force microscopy (AFM; Supplementary Table 1 and Supplementary Fig. 1) and confirmed by observation of inter-particle distances of bonded tetramer particles (Supplementary Fig. 2). AFM was performed on dried particles on a glass slide. We locate particles where the protruding patch is pointing upwards, so they are suitable for analysis. Very good consistency is observed: the mean of the inter-particle distance distribution at 3.8 μm reflects the AFM-determined particle size (particle diameter 3.7 μm) plus twice the patch height (2 × 48 nm = 96 nm), plus the (short) interaction range, while the standard deviation of σ = 0.11 μm reflects the estimated particle size variation (100 nm) and twice the variation of the patch height (2 × 5 nm = 10 nm).

A mixture of di- and tetrapatch particles (of number ratio 1:3) is dispersed in a binary solvent of 25% 2,6-lutidine (≥99.0%, Sigma-Aldrich) and 75% MilliQ water with 1 mM MgSO4 (≥99.5%, Sigma-Aldrich). The particles are washed several times in the water–lutidine mixture. The resulting particle dispersion is injected into a glass capillary (Vitrotubes, Rectangle Boro Tubing 0.20 × 2.00 mm) and sealed with teflon grease (Krytox GPL-205).

Particles are left to sediment to the bottom of the sample at room temperature before measurements. We choose a particle concentration such that particles cover between 10 and 15% of the surface after sedimentation. We then switch on the critical Casimir attraction to assemble the structures.

We use a well-controlled temperature stage in combination with an objective heating element to obtain a relative temperature accuracy of ~0.01 °C, with minimal temperature gradients.

In an experiment, we heat the sample to 0.04 °C below the phase separation temperature of the lutidine–water mixture, inducing critical Casimir attraction between patches. The structures then grow by two-dimensional diffusion in the plane. No mixing is necessary. After several hours of equilibration, we investigate the structures using a 63× oil-immersion objective.

For the catalytic experiments, we assemble the alkenes in the bulk. This is because in these experiments, the walls are treated to exhibit strong particle–wall interaction. To investigate their ‘catalytic effect’, we assemble the particles in the bulk, and then let the entire assembled structure sediment to the bottom, so that after contact with the wall, the catalytic process can begin. This is done by first letting the particles (without critical Casimir force) sediment on one side of the capillary, turning the capillary upside down, and then switch on the Casimir interaction so that rings form in the bulk (while sedimenting to the other side of the capillary). This process allows us to study the catalytic effect on cyclopentane rings that have assembled in the bulk.

We image the assembled structures using confocal microscope image stacks, alternating with bright-field images. To follow a colloidal molecule in time, we acquire ~100 image stacks and bright-field images during a time interval of 12 min. The bright-field images are processed using particle tracking software (Trackpy32) to determine the centre of the patchy particle in the horizontal plane. The 3D locations of the fluorescent features are determined using the same software. The colloidal ring structures were then detected using an image analysis algorithm, and confirmed by manual inspection of the video.

We use the method by Cremer and Pople25 to find the ring puckering coordinates for each particle configuration. This is done by expressing the positions Rj of the particles in a new coordinate system. In this new coordinate system, the average plane through all particles forms the xy-plane, and the z-displacement zj is the distance of a particle from this plane.

This new coordinate system must satisfy

meaning that the origin of the system is at the centre of mass of the five particles. We further enforce

which means all z-displacement zj in the new coordinate systems average out to zero: the mean plane must be the average of the particles. We choose the unit vector l going through any of the five particles as a convenient x-axis. Unfortunately, Eqs. (1) and (2) do not yet yield a unique average plane through the five particles: multiple planes can be drawn that will all satisfy the above conditions. Therefore, we further impose the conditions

which will fix the plane uniquely. In the case of a regular planar polygon, Eqs. (3) and (4) correspond to the condition that a small displacement zj will not lead to angular momentum. The same conditions may be used more generally for bigger displacements, and any bond lengths or angles. Also, this condition will yield the same plane irrespective of which particle is chosen as j = 1. Combining all these conditions, we can now determine the orientation of the mean plane for the position vectors Rj:

is perpendicular to \({\bf{R}}^{\prime}\) and R″. Thus, we can use n to define the new z-axis. Finally, the y-axis, defined by unit vector m, can now be found by simply taking the crossproduct of l and n. By linearly displacing coordinates to this new coordinate system, we can determine the ring’s puckering coordinates.

The ring puckering amplitude q and phase ϕ for a ring of five particles are defined as

We can simply calculate the puckering coordinates from the transformed coordinates for every snapshot.

The MC calculations are performed with a potential constructed from the universal scaling behaviour of the critical Casimir interactions, and benchmarked onto experimentally measured static and dynamic equilibrium properties of a dipatch system at various temperatures. As the dipatch and tetrapatch particles are synthesized from the same materials, immersed in the same solvent, and differ only in their physical dimensions, their potentials are similar. The bond directionality and flexibility of the patchy particles is incorporated via the switching function S (see Supplementary Note 5).

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request. Source data are provided with this paper.

The codes of the computer simulations and data analysis are available from the corresponding author upon request.

Manoharan, V. N. Colloidal matter: packing, geometry, and entropy. Science 349, 12537511–12537518 (2015).

Diao, Y. Y. & Liu, X. Y. Controlled colloidal assembly: experimental modeling of general crystallization and biomimicking of structural color. Adv. Funct. Mater. 22, 1354–1375 (2012).

Tang, Z., Kotov, N. A. & Giersig, M. Spontaneous organization of single CdTe nanoparticles into luminescent nanowires. Science 297, 237–240 (2002).

Article  ADS  CAS  Google Scholar 

Liu, Y. et al. Lysozyme protein solution with an intermediate range order structure. J. Phys. Chem. B 115, 7238–7247 (2011).

Tracey, D. F., Noya, E. G. & Doye, J. P. K. Programming patchy particles to form complex periodic structures. J. Chem. Phys. 151, 224506 (2019).

Romano, F., Sanz, E. & Sciortino, F. Crystallization of tetrahedral patchy particles in silico. J. Chem. Phys. 134, 174502 (2011).

Noya, E. G., Zubieta, I., Pine, D. J. & Sciortino, F. Assembly of clathrates from tetrahedral patchy colloids with narrow patches. J. Chem. Phys. 151, 094502 (2019).

Wang, Y. et al. Colloids with valence and specific directional bonding. Nature 491, 51–55 (2012).

Article  ADS  CAS  Google Scholar 

Soto, R. & Golestanian, R. Self-assembly of active colloidal molecules with dynamic function. Phys. Rev. E 91, 052304 (2015).

Ni, S., Leemann, J., Buttinoni, I., Isa, L. & Wolf, H. Programmable colloidal molecules from sequential capillarity-assisted particle assembly. Sci. Adv. 2, e1501779 (2016).

Elacqua, E., Zheng, X., Shillingford, C., Liu, M. & Weck, M. Molecular recognition in the colloidal world. Acc. Chem. Res. 50, 2756–2766 (2017).

Gong, Z., Hueckel, T., Yi, G. R. & Sacanna, S. Patchy particles made by colloidal fusion. Nature 550, 234–238 (2017).

Rouet, P.-E., Chomette, C., Duguet, E. & Ravaine, S. Colloidal molecules from valence-endowed nanoparticles by covalent chemistry. Angew. Chem. Int. Ed. 57, 15754–15757 (2018).

Marino, E. et al. Controlling superstructure–property relationships via critical Casimir assembly of quantum dots. J. Phys. Chem. C 123, 13451–13457 (2019).

Nguyen, V. D., Faber, S., Hu, Z., Wegdam, G. H. & Schall, P. Controlling colloidal phase transitions with critical Casimir forces. Nat. Commun. 4, 1584 (2013).

Dang, M. T., Verde, A. V., Nguyen, V. D., Bolhuis, P. G. & Schall, P. Temperature-sensitive colloidal phase behavior induced by critical Casimir forces. J. Chem. Phys. 139, 094903 (2013).

Veen, S. J. et al. Colloidal aggregation in microgravity by critical Casimir forces. Phys. Rev. Lett. 109, 248302 (2012).

G. Stuij, S., Labbé-Laurent, M., E. Kodger, T., Maciołek, A. & Schall, P. Critical Casimir interactions between colloids around the critical point of binary solvents. Soft Matter 13, 5233–5249 (2017).

Sciortino, F., Bianchi, E., Douglas, J. F. & Tartaglia, P. Self-assembly of patchy particles into polymer chains: a parameter-free comparison between Wertheim theory and Monte Carlo simulation. J. Chem. Phys. 126, 194903 (2007).

Hendrickson, J. B. Molecular geometry. I. Machine computation of the common rings. J. Am. Chem. Soc. 83, 4537–4547 (1961).

Kilpatrick, J. E., Pitzer, K. S. & Spitzer, R. The thermodynamics and molecular structure of cyclopentane. J. Am. Chem. Soc. 69, 2483–2488 (1947).

Poupko, R., Luz, Z. & Zimmermann, H. Pseudorotation in cyclopentane. An experimental determination of the puckering amplitude by NMR in oriented solvents. J. Am. Chem. Soc. 104, 5307–5314 (1982).

Ocola, E. J., Bauman, L. E. & Laane, J. Vibrational spectra and structure of cyclopentane and its isotopomers. J. Phys. Chem. A 115, 6531–6542 (2011).

Kowalewski, P., Frey, H.-M., Infanger, D. & Leutwyler, S. Probing the structure, pseudorotation, and radial vibrations of cyclopentane by femtosecond rotational raman coherence spectroscopy. J. Phys. Chem. A 119, 11215–11225 (2015).

Cremer, D. & Pople, J. A. General definition of ring puckering coordinates. J. Am. Chem. Soc. 97, 1354–1358 (1975).

Flaherty, D. W., Uzun, A. & Iglesia, E. Catalytic ring opening of cycloalkanes on Ir Clusters: alkyl substitution effects on the structure and stability of C–C bond cleavage transition states. J. Phys. Chem. C 119, 2597–2613 (2015).

Du, H., Fairbridge, C., Yang, H. & Ring, Z. The chemistry of selective ring-opening catalysts. Appl. Catal. A Gen. 294, 1–21 (2005).

Zhao, Z.-J., Moskaleva, L. V. & Rösch, N. Ring-opening reactions of methylcyclopentane over metal catalysts, M = Pt, Rh, Ir, and Pd: a mechanistic study from first-principles calculations. ACS Catal. 3, 196–205 (2013).

Moskaleva, L., Chiu, C.-c, Genest, A. & Rösch, N. Transformations of organic molecules over metal surfaces: insights from computational catalysis. Chem. Rec. 16, 2388–2404 (2016).

Haddoum, S. et al. Fe-TUD-1 for the preferential rupture of the substituted CC bond of methylcyclopentane (MCP). Catal. Commun. 27, 141–147 (2012).

Newton, A. C., Kools, R., Swenson, D. W. & Bolhuis, P. G. The opposing effects of isotropic and anisotropic attraction on association kinetics of proteins and colloids. J. Chem. Phys. 147, 155101 (2017).

Allan, D., Keim, N., Caswell, T., van der Wel, C. & Trackpy community. Trackpy v0.3.0. Zenodo http://zenodo.org/record/34028 (2014).

P.S. acknowledges support by a Vici Fellowship from the Netherlands Organization for Scientific Research (NWO). P.S. and P.B. acknowledge support by grant 680.91.124 from NWO. We thank Michele Zanini for the AFM data and analysis.

Institute of Physics, University of Amsterdam, Amsterdam, The Netherlands

PJM Swinkels, SG Stuij, N. Ruffino, B. van der Linden & P. ​​Schall

Molecular Design Institute, Department of Chemistry, New York University, New York, NY, USA

Z. Gong & S. Sacanna

Van ’t Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, The Netherlands

H. Jonas, PG Bolhuis & S. Woutersen

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

P.J.M.S., S.G.S., S.W., and P.S. conceived the study. P.J.M.S. performed the experiments with help of N.R. and B.vd.L., and analysed the data. H.J. and P.G.B. performed the simulations. Z.G. and S.S. made the particles. P.J.M.S. and P.S. wrote the manuscript. All authors discussed the data and reviewed the manuscript.

The authors declare no competing interests.

Peer review information Nature Communications thanks Ivo Buttinoni, Lucio Isa and the other anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Swinkels, P.J.M., Stuij, S.G., Gong, Z. et al. Revealing pseudorotation and ring-opening reactions in colloidal organic molecules. Nat Commun 12, 2810 (2021). https://doi.org/10.1038/s41467-021-23144-6

DOI: https://doi.org/10.1038/s41467-021-23144-6

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Nature Communications (Nat Commun) ISSN 2041-1723 (online)

Revealing pseudorotation and ring-opening reactions in colloidal organic molecules | Nature Communications

Cyclopentane Insulation Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.