In silico geometric and energetic data of all possible simple rotamers made of non-metal elements.

This article presents theoretical data on geometric and energetic features of halogenated rotamers of the following backbone structures, C-C, N-N, P-P, O-O, S-S, N-P, O-S, C-N, C-P, C-O, C-S, N-O, N-S, P-O and P-S. The data is considered to be comprehensive combinations of non-metal elements in the form abcx-ydef whereby a,b,c,d,e,f are halogen (fluorine to iodine), hydrogen or a lone pair and x,y are carbon, nitrogen, phosphorus, oxygen and sulfur. Data were obtained from ab initio geometry optimization and frequency calculations at HF, B3LYP, MP2 and CCSD levels of theory on 6-311++G(d,p) basis set. In total, 8535 non-enantiomeric structures were produced by custom-made codes in Mathematica and Q-Chem quantum chemical package. Extracted geometric and energetic data as well as raw output files, codes and scripts associated with the data production are presented in the data repository.

Specifications table

Value of the data

The origin of energetic preference for staggered structure in pan class="Chemical">ethane [1], [2], [3], [4], [5] and gauche structure in class="Chemical">pan class="Chemical">1,2-difluoroethane [6], [7], [8], [9], [10], [11] has long been debated and sometimes controversial. The comprehensive data set presented in this article fills in the gap in the literature and can be used for further analysis and discussion in relevant topics such as gauche effect [7,8,12] and bent bond [7,13].

Similar to cis effect where the cis or (Z) isomer is more stable than trans or (E) isomer [14] and relative stability of positional isomers of substituted pan class="Chemical">benzenes [15], gauche effect is demonstrated in this data set by many examclass="Chemical">ples where steric hindrance alone fails to account for the observed relative stability trend.

For reference purpose, 15665 rotamers are identified with internal numbering, SMILES and pan class="Chemical">PubChem CID. (Out of 15665 rotamers, 1713 rotamers (11%) are identified with CID, of which only 631 are unique.) These can be used in future theoretical or exclass="Chemical">perimental work involving two-center non-class="Chemical">pan class="Chemical">metal rotamers.

Source codes and raw data are available for reproduction of the work and further analysis. For example, molecular dipole moment and vibrational spectrum can be extracted from the raw output. Source codes can be used to generate molecules of related classes for further calculation.

Data description

There are 15 folders for C—C, N—N, class="Chemical">P—class="Chemical">pan class="Chemical">P, O—O, S—S, N—P, O—S, C—N, C—P, C—O, C—S, N—O, N—S, P—O and P—S. In each folder, there are four subfolders for four different methodologies, HF, B3LYP, MP2 and CCSD. In addition to raw output files (.out) and geometry in Z-matrix and Cartesian coordinate format (.xyz), the following summary table files (.csv) are provided in each subfolder:

A single csv file in xyz subfolder containing geometric data of 7 bond lengths in Å, 12 bond angles and 9 torsional angles in degree (If lone pair(s) are involved, there will be less numbers of geometric parameters and ‘de’ is shown in place of a numerical value.)

Energetic data, in separate csv files, include electronic energy (Eelec) in a.u. (Hartree), thermal correction to enthalpy (Hcorr) in kcal mol−1, zero-point vibrational energy (EZpan class="Chemical">PE) in kcal mol−1 and entroclass="Chemical">py (S) in cal mol−1 K−1.

An example of these data is shown in Fig. 1. Names for compounds exist in two different formats and due to symmetry, there are up to six ways to write these out regardless of the format. Therefore, a rotamer name may not exactly match a file name in many instances. Source codes, scripts and examples are provided in a separate folder.

Fig. 1

An example of data for BrClFCCBrICl rotamer calculated at B3LYP/6-311++G(d,p) level of theory.

An example of data for BrClFCCBrICl rotamer calculated at B3LYpan class="Chemical">P/6-311++G(d,class="Chemical">p) level of theory.

Experimental design, materials, and methods

Exhaustive listing of all rotamers can be done in many different approaches. We completed our comprehensive lists of all rotamers by extending the approach we have used for substituted class="Chemical">benzenes [15]. Rotamers as viewed by Newman class="Chemical">projection can be equivalent to substituted class="Chemical">pan class="Chemical">benzenes with two additional conditions. First, the list of substituent elements must include a lone pair of electrons. Second, rotamers are less symmetric compared to benzenes with regards to rotation and flipping. Q-Chem 5.2.1 [16], IQmol 2.13 [17] and Wolfram Mathematica 12.0 [18] were used in the same way as described previously [15]. In addition to compounds in Table 1, Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Preliminary calculations were also completed for all combination of single atom of C, N, P, O, S and hydrogen/halogen atoms from F to I.

Table 1

List of 2675 possible C—C rotamers in 210 formulas (1505 non-enantiomeric rotamers).

				Rotamers per formula		Number of rotamers
Number of elements	Empirical formula (Distribution of elements)	Number of empirical formulas	Rotamer structure	Inclusive of enantiomers	Exclusive of enantiomers	Inclusive of enantiomers	Exclusive of enantiomers
1	C2α6 (6)	5	α3C−Cα3	1	1	5	5
2	C2α5β (1–5)	20	α2βC−Cα3	1	1	20	20
	C2α2β4 (2–4)	20	αβ2C−Cαβ2	3	2	60	40
			α2βC−Cβ3	1	1	20	20
	C2α3β3 (3–3)	10	αβ2C−Cα2β	3	2	30	20
			α3C−Cβ3	1	1	10	10
3	C2αβγ4 (1–1–4)	30	αβγC*−Cγ3	2	1	60	30
			αγ2C−Cβγ2	3	2	90	60
	C2αβ2γ3 (1–2–3)	60	αβγC*−Cβγ2	6	3	360	180
			αγ2C−Cβ2γ	3	2	180	120
			αβ2C−Cγ3	1	1	60	60
	C2α2β2γ2 (2–2–2)	10	αβγC*−C*αβγ (meso compound)	9	5	90	50
			αβ2C−Cαγ2 α2βC−Cβγ2 α2γC−Cβ2γ	9	6	90	60
4	C2αβγδ3 (1–1–1–3)	20	αβδC*−Cγδ2 αγδC*−Cβδ2 βγδC*−Cαδ2	18	9	360	180
			αβγC*−Cδ3	2	1	40	20
	C2αβγ2δ2 (1–1–2–2)	30	αγδC*−C*βγδ	12	6	360	180
			αβγC*−Cγδ2 αβδC*−Cγ2δ	12	6	360	180
			αγ2C−Cβδ2 αδ2C−Cβγ2	6	4	180	120
5	C2αβγδε2 (1–1–1–1–2)	5	αβεC*−C*γδε αγεC*−C*βδε αδεC*−C*βγε	36	18	180	90
			αβγC*−Cδε2 αβδC*−Cγε2 αγδC*−Cβε2 βγδC*−Cαε2	24	12	120	60

Table 2

List of 975 possible N—N (or P—P) rotamers in 70 formulas (500 non-enantiomeric rotamers).

				Rotamers per formula		Number of rotamers
Number of elements	Empirical formula (Distribution of elements)	Number of empirical formulas	Rotamer structure	Inclusive of enantiomers	Exclusive of enantiomers	Inclusive of enantiomers	Exclusive of enantiomers
1	N2α4 (4)	5	α2N−Nα2	3	2	15	10
2	N2αβ3 (1–3)	20	αβN*−Nβ2	6	3	120	60
	N2α2β2 (2–2)	10	αβN*−N*αβ (meso compound)	9	5	90	50
			α2N−Nβ2	3	2	30	20
3	N2αβγ2 (1–1–2)	30	αγN*−N*βγ	12	6	360	180
			αβN*−Nγ2	6	3	180	90
4	N2αβγδ (1–1–1–1)	5	αβN*−N*γδ αγN*−N*βδ αδN*−N*βγ	36	18	180	90

Table 3

List of 45 possible O—O (or S—S) rotamers in 15 formulas (30 non-enantiomeric rotamers).

				Rotamers per formula		Number of rotamers
Number of elements	Empirical formula (Distribution of elements)	Number of empirical formulas	Rotamer structure	Inclusive of enantiomers	Exclusive of enantiomers	Inclusive of enantiomers	Exclusive of enantiomers
1	O2α2 (2)	5	αO−Oα	3	2	15	10
2	O2αβ (1–1)	10	αO−Oβ	3	2	30	20

Table 4

List of 1875 possible N—P rotamers in 70 formulas (950 non-enantiomeric rotamers).

				Rotamers per formula		Number of rotamers
Number of elements	Empirical formula (Distribution of elements)	Number of empirical formulas	Rotamer structure	Inclusive of enantiomers	Exclusive of enantiomers	Inclusive of enantiomers	Exclusive of enantiomers
1	NPα4 (4)	5	α2N−Pα2	3	2	15	10
2	NPαβ3 (1–3)	20	αβN*−Pβ2 αβP*−Nβ2	12	6	240	120
	NPα2β2 (2–2)	10	αβN*−P*αβ	12	6	120	60
			α2N−Pβ2 α2P−Nβ2	6	4	60	40
3	NPαβγ2 (1–1–2)	30	αγN*−P*βγ αγP*−N*βγ	24	12	720	360
			αβN*−Pγ2 αβP*−Nγ2	12	6	360	180
4	NPαβγδ (1–1–1–1)	5	αβN*−P*γδ (42)= 6 structures	72	36	360	180

Table 5

List of 75 possible O—S rotamers in 15 formulas (50 non-enantiomeric rotamers).

				Rotamers per formula		Number of rotamers
Number of elements	Empirical formula (Distribution of elements)	Number of empirical formulas	Rotamer structure	Inclusive of enantiomers	Exclusive of enantiomers	Inclusive of enantiomers	Exclusive of enantiomers
1	OSα2 (2)	5	αO−Sα	3	2	15	10
2	OSαβ (1–1)	10	αO−Sβ αS−Oβ	6	4	60	40

Table 6

List of 3125 possible C—N (or C—P) rotamers in 126 formulas (1625 non-enantiomeric rotamers).

				Rotamers per formula		Number of rotamers
Number of elements	Empirical formula (Distribution of elements)	Number of empirical formulas	Rotamer structure	Inclusive of enantiomers	Exclusive of enantiomers	Inclusive of enantiomers	Exclusive of enantiomers
1	CNα5 (5)	5	α3C−Nα2	1	1	5	5
2	CNαβ4 (1–4)	20	αβ2C−Nβ2	3	2	60	40
			β3C−N*αβ	2	1	40	20
	CNα2β3 (2–3)	20	αβ2C−N*αβ	6	3	120	60
			α2βC−Nβ2	3	2	60	40
			β3C−Nα2	1	1	20	20
3	CNαβγ3 (1–1–3)	30	αβγC*−Nγ2	6	3	180	90
			αγ2C−N*βγ βγ2C−N*αγ	12	6	360	180
			γ3C−N*αβ	2	1	60	30
	CNαβ2γ2 (1–2–2)	30	αβγC*−N*βγ	12	6	360	180
			β2γC−N*αγ βγ2C−N*αβ	12	6	360	180
			αβ2C−Nγ2 αγ2C−Nβ2	6	4	180	120
4	CNαβγδ2 (1–1–1–2)	20	αβδC*−N*γδ αγδC*−N*βδ βγδC*−N*αδ	36	18	720	360
			αβγC*−Nδ2	6	3	120	60
			αδ2C−N*βγ βδ2C−N*αγ γδ2C−N*αβ	18	9	360	180
5	CNαβγδε (1–1–1–1–1)	1	αβγC*−N*δε (53)= 10 structures	120	60	120	60

Table 7

List of 625 possible C—O (or C—S) rotamers in 70 formulas (375 non-enantiomeric rotamers).

				Rotamers per formula		Number of rotamers
Number of elements	Empirical formula (Distribution of elements)	Number of empirical formulas	Rotamer structure	Inclusive of enantiomers	Exclusive of enantiomers	Inclusive of enantiomers	Exclusive of enantiomers
1	COα4 (4)	5	α3C−Oα	1	1	5	5
2	COαβ3 (1–3)	20	αβ2C−Oβ	3	2	60	40
			β3C−Oα	1	1	20	20
	COα2β2 (2–2)	10	αβ2C−Oα α2βC−Oβ	6	4	60	40
3	COαβγ2 (1–1–2)	30	αβγC*−Oγ	6	3	180	90
			αγ2C−Oβ βγ2C−Oα	6	4	180	120
4	COαβγδ (1–1–1–1)	5	αβγC*−Oδ αβδC*−Oγ αγδC*−Oβ βγδC*−Oα	24	12	120	60

Table 8

List of 375 possible N—O (or N—S, P—O, P—S) rotamers in 35 formulas (200 non-enantiomeric rotamers).

				Rotamers per formula		Number of rotamers
Number of elements	Empirical formula (Distribution of elements)	Number of empirical formulas	Rotamer structure	Inclusive of enantiomers	Exclusive of enantiomers	Inclusive of enantiomers	Exclusive of enantiomers
1	NOα3 (3)	5	α2N−Oα	3	2	15	10
2	NOαβ2 (1–2)	20	αβN*−Oβ	6	3	120	60
			β2N−Oα	3	2	60	40
3	NOαβγ (1–1–1)	10	αβN*−Oγ αγN*−Oβ βγN*−Oα	18	9	180	90

List of 975 possible N—N (or pan class="Chemical">P—class="Chemical">pan class="Chemical">P) rotamers in 70 formulas (500 non-enantiomeric rotamers).

List of 1875 possible N—pan class="Chemical">P rotamers in 70 formulas (950 non-enantiomeric rotamers).

List of 3125 possible C—N (or C—pan class="Chemical">P) rotamers in 126 formulas (1625 non-enantiomeric rotamers).

List of 375 possible N—O (or N—S, pan class="Chemical">P—O, class="Chemical">pan class="Chemical">P—S) rotamers in 35 formulas (200 non-enantiomeric rotamers).

Table 1, Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8 provide a comprehensive listing of all rotamers considered in this work. The listing is first arranged by the number of substituent elements and pattern of empirical formulas. An explanation on how to calculate the number of chemical empirical formulas in each table is given in Table 9. Rotamer structures are also listed for each pattern. Each rotamer structure can be rotated three times unless it is symmetric (cannot rotate) or has chiral center(s) (×2 for each center). Asterisks (*) shown in Tables 1, 2, 4, Table 6, Table 7, Table 8 indicate chiral centers. There are two special cases of meso compounds in Tables 1 and 2 which have a reduced number of rotamer structures. Similar symmetrical cases were also found in previous study of substituted pan class="Chemical">benzenes [15]. Since enantiomeric structures are identical in energy, only one of the two enantiomeric structures is considered for each class="Chemical">pair. Table 10 class="Chemical">provides an overview of all comclass="Chemical">putational jobs described in this class="Chemical">paclass="Chemical">per.

Table 9

Examples for number of empirical formula calculation.

Number of elements	Empirical formula (Distribution of elements)	Number of empirical formulasa
3	C2αβ2γ3 (1–2–3)	k = 3, n1 = 1, n2 = 1, n3 = 1 therefore (53)3!1!1!1! = 60
3	C2α2β2γ2 (2–2–2)	k = 3, n1 = 3 therefore (53)3!3! = 10
4	C2αβγ2δ2 (1–1–2–2)	k = 4, n1 = 2, n2 = 2 therefore (54)4!2!2! = 30
5	C2αβγδε2 (1–1–1–1–2)	k = 5, n1 = 4, n2 = 1 therefore (55)5!4!1! = 5

The number of empirical formulas is calculated by using the expression , where• 5 is the number of possible substituent elements (H, F, Cl, Br, I),• k is the actual number of substituent elements and• ∏n! is the product of the factorial of the number of substituent elements with the same subscript.

Table 10

Summary of 43450 computational jobs included in this paper (opt for geometry optimization and freq for frequency calculation).

Class of compound	Number of rotamers		HF		B3LYP		MP2		CCSD
	All	non-enantiomeric	opt	freq	opt	freq	opt	freq	opt	freq
Preliminary CNPOS	175	170	All	All	All	All	All	–	All	–
C—C	2675	1505	All	All	All	All	All	–	23	–
N—N	975	500	All	All	All	All	All	–	30	–
P—P	975	500	All	All	All	All	All	–	30	–
O—O	45	30	All	All	All	All	All	All	All	–
S—S	45	30	All	All	All	All	All	All	All	–
C—N	3125	1625	All	All	All	All	All	–	33	–
C—P	3125	1625	All	All	All	All	All	–	33	–
C—O	625	375	All	All	All	All	All	–	33	–
C—S	625	375	All	All	All	All	All	–	33	–
N—P	1875	950	All	All	All	All	All	–	34	–
N—O	375	200	All	All	All	All	All	–	34	–
N—S	375	200	All	All	All	All	All	–	34	–
P—O	375	200	All	All	All	All	All	–	34	–
P—S	375	200	All	All	All	All	All	–	34	–
O—S	75	50	All	All	All	All	All	All	All	–
Total	15840	8535	8535	8535	8535	8535	8535	110	665	–

All optimization jobs converged in Q-Chem before geometry information were extracted. The converged results are not necessarily the same conformer as the input. For example, some anti forms are turned into gauche forms during the geometry optimization. Almost all of the converged rotamers were confirmed to be local minima by frequency calculation (16,729 out of 16,840 jobs). However, there were 111 frequency jobs with exactly one imaginary frequency ranging from 4.9i cm−1 to 513.1i cm−1. The list of files is provided in data folder. Most of them (74) are HF calculations and the rest (37) are B3LYclass="Chemical">P. Uclass="Chemical">pon closer insclass="Chemical">pection, all rotamers of O—O, S—S and O—S with an imaginary frequency (55) are in an anti form. Similarly, 28 of 34 rotamers of N—N, class="Chemical">pan class="Chemical">P—P and N—P with an imaginary frequency are in an anti form from the perspective of the two lone pairs. As these observations suggest that the anti form is not stable, gauche effect is evident in these classes of compounds.

Competing Interests

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Subject	Chemistry
Specific subject area	Physical and theoretical chemistry/spectroscopy
Type of data	Tables and Q-Chem output files
How data were acquired	Quantum chemical computation on Q-Chem 5.2.1, developer version
Data format	Raw and analyzed
Parameters for data collection	Hartree-Fock (HF)/6-311++G(d,p),
	Becke, 3-parameter, Lee–Yang–Parr (B3LYP)/6-311++G(d,p),
	Second order Møller–Plesset perturbation theory (MP2)/6-311++G(d,p), Coupled Cluster Singles and Doubles (CCSD)/6-311++G(d,p)
Description of data collection	Data were obtained from ab initio geometry optimization and frequency calculations. In total, 8535 non-enantiomeric structures were produced and processed by custom-made codes.
Data source location	Mahidol University, Salaya, ThailandLatitude and longitude: 13.792790, 100.325707
Data accessibility	Repository name: mendeley.comData identification number: DOI: 10.17632/m2h6yg9nzpDirect URL to data: https://data.mendeley.com/datasets/m2h6yg9nzp