Skip to main content
zenodoopen

Supplementary material: Picosecond pulse-shaping for strong three-dimensional field-free alignment of generic asymmetric-top molecules

<p><strong>Supplementary material to the manuscript <em>&quot;Picosecond pulse-shaping for strong three-dimensional field-free alignment of generic asymmetric-top molecules&quot;</em> by Terry Mullins, Evangelos T. Karamatskos, Joss Wiese, Jolijn Onvlee, Arnaud Rouz&eacute;e, Andrey Yachmenev, Sebastian Trippel, and Jochen K&uuml;pper, <em>Nat Commun</em> 13, 1431 (2022). <a href="https://doi.org/10.1038/s41467-022-28951-z">https://doi.org/10.1038/s41467-022-28951-z</a>, arXiv: <a href="https://arxiv.org/abs/2009.08157">2009.08157 </a></strong></p> <ul> <li>&nbsp;<em><strong>simulations_part.z*</strong> </em>is split zip archive containing simulations data for indole molecule, such as files with rotational probability density distributions computed at different times <span class="math-tex">\(t=0..1500\)</span> ps during the laser pulse and field-free evolution, and example python scripts for data retrieval.</li> <li><strong><em>rawdata_part.z*</em></strong> is split zip archive containing raw experimental data.</li> <li><strong><em>analysis_scripts.zip</em></strong> is zip archive containing experimental analysis codes.</li> </ul> <p><strong>The <em>simulations_part.zip</em> contains the following files and folders:</strong></p> <ul> <li><em><strong>prob_density_euler_angles</strong></em> contains files <em>rotdens_av_&lt;time&gt;.gz</em> with simulated state-averaged rotational probability density distributions in terms of Euler angles for different times &lt;time&gt;, ranging from the beginning of the alignment laser pulse at&nbsp;<span class="math-tex">\(t=0\)</span> up to&nbsp;<span class="math-tex">\(t=1500\)</span> ps with a time step of 1 ps.<br> Calculations of probability density distributions were done using <a href="https://github.com/CFEL-CMI/richmol">Richmol</a> program.<br> The gzipped ASCII data files <em>rotdens_av_&lt;time&gt;.gz</em> contain in columns the values of the Euler angles&nbsp;<span class="math-tex">\(\phi,\theta,\chi\)</span> followed by the normalized probability density value.</li> <li><em><strong>prob_density_atoms_xyz</strong></em> contains files <em>monte_carlo_av_&lt;time&gt;.h5</em> with state-averaged rotational probability density distributions of all atoms in the indole molecule in terms of their Cartesian coordinates, for different times &lt;time&gt;, ranging from the beginning of the alignment pulse&nbsp;at&nbsp;<span class="math-tex">\(t=0\)</span> up to&nbsp;<span class="math-tex">\(t=1500\)</span> ps with a time step of 1 ps.<br> Structure of <em>monte_carlo_av_&lt;time&gt;.h5</em> HDF5 files:<br> Key&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Description<br> -----&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ----------------<br> &#39;C10&#39;&nbsp;&nbsp;&nbsp;&nbsp; - Cartesian coordinates of carbon atom no. 10<br> &#39;C11&#39;&nbsp;&nbsp;&nbsp;&nbsp; - Cartesian coordinates of carbon atom no. 11<br> &#39;C12&#39;&nbsp;&nbsp;&nbsp;&nbsp; - ...<br> &#39;C14&#39;&nbsp;&nbsp;&nbsp;&nbsp; - ...<br> &#39;C3&#39;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; - ...<br> &#39;C6&#39;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; - ...<br> &#39;C7&#39;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; - ...<br> &#39;C9&#39;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; - ...<br> &#39;N4&#39;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; - ...<br> &#39;H1-C3&#39;&nbsp;&nbsp; - Cartesian coordinates of a vector pointing from carbon atom no. 3 to hydrogen atom no. 1<br> &#39;H13-C11&#39; - ...<br> &#39;H15-C12&#39; - ...<br> &#39;H16-C14&#39; - ...<br> &#39;H2-N4&#39;&nbsp;&nbsp; - ...<br> &#39;H5-C7&#39;&nbsp;&nbsp; - ...<br> &#39;H8-C9&#39;&nbsp;&nbsp; - ...<br> &#39;ref_vectors&#39; - reference molecular-frame Cartesian coordinates of all atoms<br> &#39;x&#39;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; - coordinates of the x-axis of Principal Axes of Inertia Frame<br> &#39;y&#39;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; - coordinates of the y-axis of Principal Axes of Inertia Frame<br> &#39;z&#39;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; - coordinates of the z-axis of Principal Axes of Inertia Frame<br> &#39;pol_x&#39;&nbsp;&nbsp; - coordinates of the x-axis of Principal Axes of Polarizability Frame<br> &#39;pol_y&#39;&nbsp;&nbsp; - coordinates of the y-axis of Principal Axes of Polarizability Frame<br> &#39;pol_z&#39;&nbsp;&nbsp; - coordinates of the z-axis of Principal Axes of Polarizability Frame</li> <li><em><strong>indole_deflected_states.txt</strong></em> ASCII file contains initial populations of rotational states of indole in the deflected beam.<br> The following data is arranged in columns: <em>m, J</em>, <em>id</em>, <em>energy</em>, <em>normalized population</em>. The <em>J</em> and <em>m</em> are rotational quantum numbers of the total angular momentum and its <em>Z</em>-projection, the <em>id</em> number refers to the state&#39;s index in file <em>indole_energies_j0_j20.txt</em> listing rotational states of indole.</li> <li><em><strong>indole_data.py</strong></em> Python module provides basic functions to extract information from HDF5 data files <em>monte_carlo_av_&lt;time&gt;.h5</em>. It can also be used to compute alignment and orientation.</li> <li><em><strong>example_cos.py</strong></em> and <em><strong>example_dens.py</strong></em> Python scripts that demonstrate how to use <em>indole_data.py</em> module for computing and plotting alignment traces and a 2D projection of the probability density distribution, respectively.</li> <li><em><strong>monte_carlo.py</strong></em> Python script that was used to compute through Monte-Carlo sampling probability density distributions for Cartesian positions of atoms in indole (<em>monte_carlo_av_&lt;time&gt;.h5</em> files) using probability density distribution functions in Euler angles (outputs of Richmol program <em>rotdens_av_&lt;time&gt;.gz</em>).</li> </ul> <p><strong>The <em>analysis_scripts.zip</em> contains the following files and folders:</strong></p> <ul> <li><strong><em>H_Plus</em></strong> folder contains codes relevant for the analysis of H<sup>+</sup>&nbsp;ion data. <ul> <li><strong><em>analyse_full_alignment_scans.m</em></strong>: subtracts background and combines delay scan data sets together, takes account of errors.</li> <li><em><strong>calculate_resamped_df.m</strong></em>:&nbsp;called by <em>analyse_full_alignment_scans.m</em> to calculate the frequency sampling.</li> <li><em><strong>unique_mean.m</strong></em>:&nbsp;called by <em>analyse_full_alignment_scans.m</em> when combining data sets. Combines non-unique data points into a single data point.</li> </ul> </li> <li><em><strong>C_Plus2</strong></em> folder contains codes for the analysis of C<sup>2+</sup>&nbsp;ion data.&nbsp;The file descriptions are identical to those in the <em>H_Plus</em> directory.</li> <li><em><strong>intensity/calculate_intensity.m</strong></em>:&nbsp;calculates peak intensity of the laser pulse from measured parameters as well as statistical error.</li> <li><em><strong>intensity/compare_exp_sim.m</strong></em>:&nbsp;fits experimental and theoretical tomography and delay-dependent 2D projection values.</li> <li><em><strong>intensity/nir2hdf5_kHz.py</strong></em>:&nbsp;converts&nbsp;raw data files (from&nbsp;<em>rawdata_part.z*</em> archive<em>)</em>&nbsp;into hdf5 files.</li> </ul>

ShareScore

40/100

Overall dataset sharing score

Score breakdown

These five areas show where the dataset supports — or may limit — practical reuse.

Stewardship
4
Harmonization
4
Access
20
Reuse readiness
8
Engagement
4