Next Article in Journal
New Ethylated Derivatives of Sulfur- and Nitrogen-Containing Artifacts from Tenodera sinensis Egg Pod and Their Anti-Renal Fibrosis
Next Article in Special Issue
Why Do Ionic Surfactants Significantly Alter the Chemiluminogenic Properties of Acridinium Salt?
Previous Article in Journal
Effect of Stepwise Exposure to High-Level Erythromycin on Anaerobic Digestion
Previous Article in Special Issue
Application of Strain Engineering in Solar Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Criteria to Provide a Locality/Normality Degree in Molecules and Their Relevance in Physical Chemistry

by
Eduardo Suárez
,
Oscar Guzmán-Juárez
and
Renato Lemus
*
Instituto de Ciencias Nucleares, Universidad Nacional Autónoma de México, A.P. 70-543, Circuito Exterior, C. U., Mexico City 04510, Mexico
*
Author to whom correspondence should be addressed.
Molecules 2024, 29(15), 3490; https://doi.org/10.3390/molecules29153490
Submission received: 19 June 2024 / Revised: 18 July 2024 / Accepted: 19 July 2024 / Published: 25 July 2024
(This article belongs to the Special Issue Exclusive Feature Papers in Physical Chemistry, 2nd Edition)

Abstract

:
In contrast to the traditional analysis of molecules using local mode behavior, where the degree of locality is given through a function in terms of Morse potential parameters, new criteria for locality/normality (LN) suitable for application to any molecular system are proposed. The approach is based on analysis of the connection between the algebraic normal and local mode representations. It is shown that both descriptions are equivalent as long as the polyad (total number of quanta) in the local representation is not conserved. The constraint of a local polyad conservation naturally provides a criterion for assigning an LN degree in quantitative form, without an analogue in configuration space. The correlation between the different parameters reveals the physical properties of molecules. A clear connection between the LN degree (based on the fundamentals) and spectroscopic properties is also presented, suggesting a promising approach for identifying mixtures of isotopologues.

Graphical Abstract

1. Introduction

Vibrational degrees of freedom can be identified either with a local or a normal mode behavior. Until the seventies, the point of view of description in terms of normal modes dominated, due in part to the success in describing spectra through the inclusion of resonances in the Hamiltonian [1,2]. Examples of this success include the descriptions of H2O and CO2 [3,4]. This situation changed during the eighties with the advent of modern spectroscopy techniques based on lasers [5,6,7], unveiling doublets in the energy spectra of molecules involving bonds with large mass differences, although evidence of such patterns had been identified many years before [8,9,10,11]. This kind of spectra, although difficult to describe in terms of normal modes, were relatively easy to interpret in terms of interacting local oscillators [12,13,14,15,16]. Indeed, the doublets signal the presence of a local mode behavior and they were explained by the simple model of anharmonic oscillators harmonically coupled (AOHC) [17]. Regarding this model, the general features of the spectra are explained in terms of both the anharmonicity and strength of the interaction between the oscillators, as a consequence of the close relation between anharmonicity and locality [18]. It was realized that in every molecule with bonds involving large mass differences, like for instance H2O, CH4, and AsH3, the stretching degrees of freedom cn be well described in terms of interacting Morse oscillators. In contrast, the bending degrees of freedom are still treated in terms of normal modes, due to their collective features [19]. Within a set of molecules presenting a local mode behavior, it is possible to assign a locality degree using the parameter ξ defined by [12,13]
ξ = 2 π arctan λ ω x ,
where ω x stands for the Morse anharmonicity and λ corresponds to the interaction strength of the oscillators. In the local limit ( ω x large and λ small) ξ 0 , while in the normal limit ( ω x small and λ large) ξ 1 . This parameter has been calculated for several molecules following the expected behavior, in accordance with the mass ratio ligand/central-atom dictated by the structure of the molecule. In general, the criterion established by the parameter ξ is satisfied for a great variety of molecules. However, when in a series of molecules the anharmonicities are similar and the strengths do not increase in accordance with the mass difference, the local-normal mode behavior leads to unexpected results. This is the case for the series of pnictogen pyramidal hydrides XH3, with X = N, P, As, Sb, and Bi. In order to elucidate this type of behavior a new perspective on local-normal behavior was considered [20].
Even for molecules with local mode behavior, a local-normal transition may appear regarding different states of a given multiplet. This situation has been analyzed from a dynamical point of view using methods of non-linear classical mechanics [21,22,23,24,25,26,27,28,29]. Recently, the local to normal mode transition has been studied from a quite different point of view [20,30,31,32]. The basic feature of the proposal is to focus on the problem from a polyad breaking perspective. A molecule with a local mode behavior is characterized by a set of interacting levels associated with a local polyad, defined in terms of local quantum numbers. As long as this polyad is conserved, the molecule maintains a local character. This viewpoint is partial, in the sense that the polyad-conservation depends on the energy. At sufficiently high energy, the local polyad stops being conserved. However, a local behavior may still be present in a wide energy range, as long as the local polyad is conserved. From this perspective, the energy range becomes very narrow for molecules with evident normal mode behavior. This is indeed the case of the CO2 molecule, where the central carbon atom is lighter than the ending oxygen atoms. In this situation, it is not possible to define a local polyad. In fact, the concept of a polyad as a pseudo quantum number that embraces the set of interacting states is well defined in a normal mode scheme, and only when this polyad is suitable to be translated into a local scheme is the molecule said to have a local character, otherwise the concept of a local polyad is lost. The question that arises is how to measure the LN degree from this perspective. The traditional parameter (1) is not useful, because it is constrained to molecules with local mode behavior and consequently a new criterion is needed, in order to involve any possible case covering the different molecular systems. On the other hand, this new criterion is not expected to be based on a system of interacting Morse oscillators, because, in principle, molecules with strong normal mode behavior cannot be described starting with local oscillators without breaking the polyad. In this work, we propose a new general different approach based on the analysis of the limit of the normal to the local description. Following this route, it is possible to derive more than one LN criteria. Surprisingly, these criteria include only the analysis of normal modes associated with the fundamentals. In this contribution, we present in detail the analysis of the new LN criteria for two, three, four, and six oscillators. The case of two oscillators has been partially presented in the context of the isotopologues of CO2 [33], as well as the three equivalent oscillators in pyramidal molecules [20], but here both are included in detail in the context of an unified treatment that includes any molecular system. We shall show that the proposed parameters are sensitive enough to distinguish any molecule in either the extreme normal or local regimes. The four-oscillator system is included for tetrahedral molecules, to study the situation where spurious states are present, while octahedral molecules involving six oscillators are studied, because several types of interactions are involved. Embracing this set of systems, we consider every situation to test the proposed parameters. We stress that, for the first time, LN parameters are proposed with incidence in molecular properties, such as non-rigidity and resonances, but also in spectroscopy where vibrational excitations are involved.
This paper is organized as follows: The general theory to establish the different LN criteria is presented in Section 2. In Section 3, an analysis involving two oscillators is presented in detail. Section 4 is devoted to the case of three equivalent oscillators. In Section 5, tetra-atomic molecules are analyzed, while in Section 6, the stretching modes of octahedral molecules are studied. Section 7 is devoted to discussing the relation between the parameters proposed and their physical properties. Finally, the conclusions are drawn in Section 8.

2. LN Mode Criteria: General Formalism

In this section, we present the general ideas that lead to establishing several independent parameters to measure the LN degree of molecules. We start by considering the situation for an arbitrary number of equivalent oscillators, albeit for the particular case in which the normal coordinates coincide with symmetry-adapted coordinates. This is a common situation in small and medium size molecules, where the energy of the bends are quite different from the stretches. The generalization to normal coordinates involving several coordinates of the same symmetry involves GF formalism and will be discussed in the next section.
The simplest model to describe vibrational degrees of freedom consists in considering the system as a set of independent harmonic oscillators associated with normal coordinates. When normal coordinates span irreducible representations Γ of the symmetry group G, the Hamiltonian in terms of bosonic operators takes the form [34]
H ^ N = Γ ћ Ω Γ 2 γ ( A ^ Γ γ A ^ Γ γ + A ^ Γ γ A ^ Γ γ ) ,
where the sum over Γ runs over all irreducible representations contained in the reducible representation spanned by the equivalent local oscillators. The subscript index N in the Hamiltonian emphasizes the normal mode representation. From the fundamental energies and the explicit expression of the frequencies Ω Γ = G Γ Γ ( 0 ) F Γ Γ , we are able to estimate the force constants for both normal and local mode schemes
{ F Γ Γ } { f q i q j ( N ) } .
Here, G Γ Γ ( 0 ) and F Γ Γ correspond to the elements of the Wilson matrix and force constants, respectively, in the normal basis. The local force constants f q i q j ( N ) associated with the local coordinates { q i } are obtained from the force constants F Γ Γ  using the chain rule. The superscript index (N) emphasizes the precedence of the local force constants. The connection between the normal bosonic operators A ^ Γ γ and the local bosonic operators a ^ i ( a ^ i ) corresponding to the internal coordinates is given by the Bogoliubov transformation
A ^ Γ , γ = δ + Γ i m i , Γ γ a ^ i + δ Γ i m i , Γ γ a ^ i ,
where
δ ± Γ = 1 4 F Γ Γ G Γ Γ ( 0 ) f q q g q q ( 0 ) g q q ( 0 ) G Γ Γ ( 0 ) + f q q F Γ Γ ± 2 f q q g q q ( 0 ) F Γ Γ G Γ Γ ( 0 ) .
The term g q q ( 0 ) stands for the Wilson matrix elements in the local scheme. Given our consideration of equivalent diagonal operators, we have f q i q i = f q q ; g q i q i ( 0 ) = g q q ( 0 ) i . The elements | | m i , Γ γ | | stand for the coefficients connecting the normal and local coordinates through
Q Γ γ = i m i , Γ γ q i ,
which are obtained by symmetry projection [34]. The relation (4) is obtained from the definition of the operators A ^ Γ , γ in terms of the normal coordinates (6) and assuming bosonic operators for the local coordinates with frequencies ω i = f q i q i g q i q i ( 0 ) . Notice that, in Equation (5), there is no superscript index in the force constants f q q . The reason for this is that the same expression will be used to estimate the force constants in both local and normal limits. The coefficients δ ± Γ satisfy the relation δ + Γ δ Γ = 1 , as a consequence of the commutator [ A ^ Γ γ , A ^ Γ γ ] = δ Γ Γ δ γ γ . In addition, we can see that in the pure local limit δ + Γ = 1 and δ Γ = 0 , and consequently it is useful to introduce the average
δ = 1 N Γ Γ δ Γ
as an LN degree, with  N Γ standing for the number of irreducible representations contained in the subspace of N oscillators. The substitution of (4) into the Hamiltonian (2) leads to an algebraic representation of the Hamiltonian in the local scheme
H ^ L = i λ 0 ( i ) ( a ^ i a ^ i + a ^ i a ^ i ) + i j λ 1 ( i j ) ( a ^ i a ^ j + a ^ i a ^ j ) + i j λ 2 ( i j ) ( a ^ i a ^ j + a ^ i a ^ j ) ,
with coefficients
λ 0 ( i ) = Γ ћ Ω Γ 2 ( δ + Γ + δ Γ ) γ | m i , Γ γ | 2 ,
λ 1 ( i j ) = Γ ћ Ω Γ 2 ( δ + Γ + δ Γ ) γ m i , Γ γ m j , Γ γ ,
λ 2 ( i j ) = Γ ћ Ω Γ ( δ + Γ δ Γ ) γ m i , Γ γ m j , Γ γ .
These expressions provide a deep insight into the problem, since it implies that, in the limit δ 0 , the coefficients λ 2 ( i j ) vanish. Consequently, only when δ is negligible is the total number of local quanta P L = i a ^ i a ^ i conserved. We will shortly return to this point. In practice, it is more convenient to carry out the explicit substitution of (5) to express these parameters in the following form
λ 0 ( i ) = ћ 4 Γ F Γ Γ g q q ( 0 ) f q q + G Γ Γ ( 0 ) f q q g q q ( 0 ) γ | m i , Γ γ | 2 ,
λ 1 ( i j ) = ћ 4 Γ F Γ Γ g q q ( 0 ) f q q + G Γ Γ ( 0 ) f q q g q q ( 0 ) γ m i , Γ γ m j , Γ γ ,
λ 2 ( i j ) = ћ 4 Γ F Γ Γ g q q ( 0 ) f q q G Γ Γ ( 0 ) f q q g q q ( 0 ) γ m i , Γ γ m j , Γ γ .
For equivalent oscillators the coefficients λ 0 ( i ) are independent of the oscillator, and are reduced to λ 0 ( i ) = ћ ω 0 / 2 with ω 0 = ( f q q g q q ( 0 ) ) 1 / 2 for any oscillator. The parameters λ 1 ( i j ) and λ 2 ( i j ) can be simplified to the form
λ 1 ( i j ) = ћ ω 0 2 ( x f ( i j ) + x g ( i j ) ) ; λ 2 ( i j ) = ћ ω 0 2 ( x f ( i j ) x g ( i j ) ) ,
with definitions
x f ( i j ) = f q i q j f q i q i ; x g ( i j ) = g q i q j ( 0 ) g q i q i ( 0 ) .
Not all of these parameters are different. Some of them are expected to be equal, depending on the symmetry dictated by the geometrical conformation. This information is contained in the matrix | | m i , Γ γ | | . To simplify the notation, we shall use x f ( x g ) when only one type of interaction is present and x f ( x g ) for a second type of interaction. It is important to notice that, as long as spurious modes are not present, the same Hamiltonian (8) is obtained starting from the quadratic Hamiltonian in local coordinates and momenta, which explains the subscript index L. The presence of redundancies imposes constraints on the local representation (8), as will be shown in Section 5.
The Hamiltonians (2) and (8) are equivalent. Their difference lies in the representation. Both Hamiltonians conserve the normal total number of quanta P ^ N = Γ γ A ^ Γ γ A ^ Γ γ , but from the point of view of (8), P ^ L = i a ^ i a ^ i is not conserved. From now on, we assign the name polyad to the total number of quanta, which is justified by that fact that we shall be dealing with Hamiltonians with interactions up to second order (normal modes). In order to obtain a local polyad-conserving Hamiltonian, we have two alternatives. One possibility consists in just neglecting the interactions associated with λ 2 ( i j ) in the Hamiltonian (8) to obtain
H ^ L ( P L ) = i λ 0 ( i ) ( a ^ i a ^ i + a ^ i a ^ i ) + i > j λ 1 ( i j ) ( a ^ i a ^ j + a ^ i a ^ j ) ,
where the coefficients are kept to be identified by (10). This Hamiltonian (13) can be used to estimate the force constants
{ f q i q j ( L ) } ,
by choosing parameters to fit the fundamental energies and using the matrix representation in the local basis | n 1 n 2 = Π i | n i , with n i being the number of local quanta for the i-th oscillator. It is clear that this route is feasible as long as an evident local mode behavior is present. In Table 1, analytical expressions for the force constants obtained from (3) and (11) are displayed for the systems we shall discuss.
As a second alternative to obtain a Hamiltonian with the property [ H ^ , P ^ L ] = 0 , we may take δ Γ = 0 in (4) and consequently also in (9). In this case, the renormalization δ + Γ = 1 must be imposed to satisfy [ A ^ Γ γ , A ^ Γ γ ] = δ Γ Γ δ γ γ valid. In this case, the canonical transformation takes the form
A ^ Γ , γ = i m i , Γ γ c ^ i .
The substitution of the operators (15) into the Hamiltonian (2) leads to the polyad-conserving Hamiltonian
H ^ L ( P L ) = ω n o r i ( c ^ i c ^ i + c ^ i c ^ i ) + i j λ n o r ( i j ) ( c ^ i c ^ j + c ^ i c ^ j ) ,
with coefficients
ω n o r = Γ ћ Ω Γ 2 γ | m i , Γ γ | 2 ,
λ n o r ( i j ) = Γ ћ Ω Γ 2 γ m i , Γ γ m j , Γ γ ,
The Hamiltonians (16) and (8) look similar when λ 2 ( i j ) = 0 , but they are not the same. Two features make them different: (i) the bosonic operators c ^ i ( c ^ i ) are not strictly local, which explains the subscript index L ; and (ii) the relation between the spectroscopic parameters and the force constants is different, a fact that suggests an additional LN criterion, as we next discuss.
The sets of parameters { ω n o r , λ n o r ( i j ) } in (16) and { λ 0 ( i ) , λ 1 ( i j ) } in (13) are functions of both the force and the structure constants, although with different functional form. To establish the connection between both sets, it is convenient to recall the definitions (12), since they are expected to be small, actually vanishing at the pure local limit. The set { λ 0 ( i ) , λ 1 ( i j ) } is expected to be recovered from { ω n o r , λ n o r ( i j ) } near the local limit. The latter set can be considered a function of the variables (12), which can be put together in vector form x = { x f ( i j ) , , x g ( i j ) , } . The connection between the parameters is given by the Taylor series expansion
ω n o r = 2 λ 0 ( i ) ћ 1 + 1 2 ! x H ω ( 0 ) x ˜ + | | x | | 2 E 2 ω ( 0 , x ) ,
λ n o r ( i j ) = λ 1 ( i j ) + ω l o c 1 2 ! x H λ ( i j ) ( 0 ) x ˜ + | | x | | 2 E 2 λ ( i j ) ( 0 , x ) ,
where H ( 0 ) is the Hessian matrix evaluated at x = 0 and | | x | | 2 E 2 ( 0 , x ) denotes the error involved up to second order. It is clear that the second-order terms measure the deviation of the parameters from the local mode description, but they also indicate the degree of local polyad conservation. We thus propose the set
γ ( ω ) = | 1 2 ! x H ω ( 0 ) x ˜ | ; γ ( i j ) = | 1 2 ! x H λ ( i j ) ( 0 ) x ˜ | ,
as new parameters to provide an LN degree. In addition, the spectroscopic parameters (11) from the local scheme and the parameters involved in (16) provide different force constants, both connected through (18). This fact suggests the introduction of the following parameters to estimate the LN degree
ϵ 1 = 1 f q q ( L ) f q q ( N ) ; ϵ ( i j ) = 1 f q i q j ( L ) f q i q j ( N ) ,
since lim x 0 ϵ 1 = ϵ ( i j ) = 0 .
On the other hand, since the energy splitting of a set of degenerate equivalent oscillators is expected to be proportional to the their interaction strength, the following LN parameter may also be proposed [30,31]
ζ = | 2 π arctan E d E u ( E d + E u ) / 2 |
where E d and E u correspond to the lowest and highest energy of the multiplet.
Summarizing, we have identified four parameters, namely δ Γ , γ , { ϵ i , ϵ ( i j ) } , and ζ , which vanish in the local limit, and consequently provide a way to measure the locality(normality) of a molecule. It is convenient to emphasize that these criteria do not assume a model of interacting Morse oscillators, but only harmonic oscillators. They arise from bosonic operators and are determined by the fundamental energies. We now describe their characteristics:
δ Γ : This parameter measures the suitability of applying the polyad-conserving canonical transformation (15) for each symmetry (normal mode). For convenience, we also introduced the average (7), which takes into account the contribution of all normal modes. It is worth stressing that this parameter δ can be calculated for any molecular system, and from its definition (5) it measures the degree of locality/normality from the point of view of the normal mode scheme. In contrast, the polyad-conserving transformation (15) has been assumed in every study of the local-to-normal mode transition involving stretching degrees of freedom of molecules with a clear local mode behavior. Indeed, this assumption leads to the x-K relations [18,35,36,37,38,39,40,41,42], which stop being valid when the molecules move to a normal mode behavior.
γ s : These parameters correspond to the Hessians in (18) and they provide the quadratic approximation for the sets { ω n o r , λ n o r ( i j ) } starting from { λ 0 ( i ) , λ 1 ( i j ) } when the local parameters become small. These parameters establish LN criteria from the point of view of the local mode scheme and are expected to be correlated with δ , in some cases in a perfect linear trend when the higher order terms in the expansion (18) are neglected compared with the quadratic terms, as we shall discuss.
ϵ s : This set of parameters are introduced to see the impact of the LN degree on the estimation of the force constants and takes into account both the normal and local mode schemes, since both estimations (3) and (14) for the force constants are involved. However, in order to obtain reliable results, the force constants should not be too small, since the errors may hide the criterion.
ζ : This parameter takes into account the correlation between the strength of the interaction and the LN degree, and it is a natural parameter based on the correlation between the splitting of the fundamentals and the strength of the interaction. However, we shall prove that the parameter ζ , although intuitive, is not appropriate for establishing a LN degree.
Here, we have assumed that the development of this approach lies in the formulation of the model in terms of internal coordinates. This route can be quite elaborate from a theoretical point of view because of the calculation of the Wilson matrix and the identification of redundancies. However, since these parameters involve the concept of normal modes, the calculations are quite fast.
A comment regarding the correlation between the polyad breaking and the LN degree deserves special attention. The canonical transformation (15) conserves the local polyad, but even when the local polyad is not conserved, the Hamiltonian satisfies [ P ^ N , H L ] = 0 . Substitution of (4) into the definition of a normal polyad leads to the relation
P ^ N = ζ 0 + i β 0 ( i ) ( a ^ i a ^ i ) + i j β 1 ( i j ) ( a ^ i a ^ j + a ^ i a ^ + j ) + i β 2 ( i ) ( a ^ i 2 + a ^ i 2 ) + i j β 3 ( i j ) ( a ^ i a ^ j + a ^ i a ^ j ) ,
with coefficients given by
ζ 0 = N Γ Γ δ Γ = N Γ 2 δ ,
β 0 ( i ) = Γ γ ( δ + Γ + δ Γ ) | m i , Γ γ | 2 ,
β 1 ( i j ) = Γ γ ( δ + Γ + δ Γ ) m i , Γ γ m j , Γ γ ,
β 2 ( i ) = Γ γ δ + Γ δ Γ | m i , Γ γ | 2 ,
β 3 ( i j ) = Γ γ δ + Γ δ Γ m i , Γ γ m j , Γ γ .
The transformation (15) assumes δ Γ = 0 and δ + Γ = 1 , leading to the values ζ 0 = β 1 ( i j ) = β 2 ( i j ) = β 3 ( i j ) = 0 and β 0 ( i j ) = 1 , with P ^ N = P ^ L . The explicit behavior of these parameters will be studied later on in the context of two equivalent oscillators. At the moment, we just reinforce the argument that local polyad breaking is strongly correlated with the LN degree.
In the next sections, we present an analysis for different numbers of oscillators. The aim is to show that the correlation between different parameters unveils physical properties. Due to its importance, the case of two oscillators will be studied in detail, including the case of non-equivalent oscillators. The latter illustrates the way our approach is modified when the normal mode coordinates do not coincide with symmetry-adapted coordinates.

3. LN Degree in Triatomic Molecules

In this section, we revisit the stretching degrees of freedom of triatomic molecules. Since this system has already been discussed in part [30,31], we present only its salient features, although in modified form in accordance with the general framework presented in the previous section. First, we consider the case of equivalent oscillators.

3.1. Equivalent Oscillators

For two equivalent oscillators, the irreducible representations are two, which we name as Γ = g , u . The expressions for the force constants (3) are [30,31]
F g g = f r r + f r r ; G g g ( 0 ) = g r r ( 0 ) + g r r ( 0 ) ;
F u u = f r r f r r , G u u ( 0 ) = g r r ( 0 ) g r r ( 0 ) ,
from which we obtain the force constants f q i q j ( N ) provided by Table 1, with μ 2 = u . The parameters (5) can be simplified to
δ ± Γ = 1 4 ( r Γ ± 1 ) 2 r Γ ; r Γ = 1 + ( ) σ Γ x f 1 + ( ) σ Γ x g ,
where σ g = 0 and σ u = 1 . The matrix | | m i , Γ γ | | corresponding to (6) takes the form
| | m i , Γ γ | | = 1 2 1 1 1 1 .
The substitution of (4) into (2) leads to an algebraic representation in terms of local bosonic operators a ^ i ( a ^ i ) :
H ^ L = λ 0 i = 1 2 ( a ^ i a ^ i + a ^ i a ^ i ) + λ 1 ( 1 ) ( a ^ 1 a ^ 2 + a ^ 1 a ^ 2 ) + λ 2 ( 1 ) ( a ^ 1 a ^ 2 + a ^ 1 a ^ 2 ) ,
with coefficients given by (11). Here, the superscript index (1) in λ ’s indicates that only one type of interaction is present. The basis for constructing the matrix representation of the Hamiltonian (27) is given by the direct product | n 1 n 2   =   | n 1     | n 2 . When the interaction strength λ 2 ( 1 ) is large, the whole space of states | n 1 n 2 is mixed, leading to a time-consuming diagonalization. In contrast, when λ 2 ( 1 ) is negligible, the Hamiltonian commutes with the operator P ^ L . Hence, if we make the approximation λ 2 ( 1 ) 0 and fit the λ 0 and λ 1 ( 1 ) to reproduce the fundamentals, we can estimate the force constants (14) from (11). The analytical expressions are given in Table 1.
A more fruitful viewpoint for obtaining a local polyad-conserving Hamiltonian is provided by realizing that λ 2 ( 1 ) also vanishes when the contribution of annihilation operators in (4) is null. To follow this route, we apply the canonical transformation (15), which when substituted into (2), leads to the Hamiltonian
H ^ L ( P L ) = ω n o r i = 1 2 ( c ^ i c ^ i + c ^ i c ^ i ) + λ n o r ( 1 ) ( c ^ 1 c ^ 2 + c ^ 1 c ^ 2 ) ,
where
ω n o r = ћ ω 0 2 1 2 ( 1 + x f ) ( 1 + x g ) + 1 2 ( 1 x f ) ( 1 x g )
λ n o r ( 1 ) = ћ ω 0 2 ( 1 + x f ) ( 1 + x g ) ( 1 x f ) ( 1 x g ) .
The relation between these coefficients and (29) is obtained through the Equation (18), which takes the form [30,31]
ω n o r = ћ ω 0 2 1 1 8 ( x g x f ) 2 + O ( x 3 ) ,
λ n o r ( 1 ) = ћ ω 0 2 x g + x f + O ( x 3 ) .
with the identification
γ ( ω ) = 1 8 ( x f x g ) 2 .
This parameter is related to λ 2 ( 1 ) appearing in the Hamiltonian (27) by γ ( ω ) = 1 2 λ 2 ( 1 ) ћ ω 0 2 , which is consistent with the fact that both of them are associated with local polyad breaking. From the expressions for the force constants obtained from Table 1, we can calculate the parameters ϵ i , i = 1 , 2 defined in (20). Finally, we should recall that we have the additional parameter ζ defined in (21), which was used in Refs. [30,31] for the case of two oscillators.
We now proceed to show the relation between these parameters for different molecules. We start by considering the series of symmetric triatomic molecules analyzed in Refs. [30,31]. In Table 2, the fundamentals as well as the force constants are displayed, while in Table 3 all the parameters suitable for measuring the LN degree are listed. In order to appreciate the behavior of the parameters, it is convenient to display the results in graphical form. In Refs. [30,31], the plot ζ vs. γ was presented. Although a general local-to-normal trend was identified, a clear correlation between the parameters was not manifested. This fact suggests that the splitting of the interacting oscillators is not necessarily a quantitative criterion for assigning an LN degree. In Figure 1, a plot of δ vs. γ ( ω ) for all the molecules included in Table 3 is depicted. Now, we can see a linear trend for the molecules near the local limit. The molecules with normal mode behavior are shifted to the upper part of the line. The reason for this is that the parameter γ ( ω ) only takes into account the first term of the expansion (30). A discussion of this behavior will be provided in Section 3.3. The importance of the correlation between δ and γ ( ω ) is two fold: first both parameters represent a consistent LN criterion and second this correlation allows unveiling molecular properties, as we shall see later on.
We now turn our attention to the correlation between the parameters ϵ 1 vs. δ and ϵ 1 vs. γ ( ω ) displayed in Figure 2. A clear linear trend is obtained, which confirms the validity of the parameter ϵ 1 for measuring the LN degree.
Regarding the correlation of ϵ 2 vs. δ , a somewhat unclear behavior is obtained for the molecules with strong local behavior. This is explained as, in the region close to locality, a weak interaction between oscillators is present, leading to small values of f r r . This fact makes it difficult to assign an LN degree; the error is of the same order of magnitude as the parameters themselves.

3.2. Non-Equivalent Bonds

In the previous analysis, we considered two equivalent oscillators, a case where the principal isotopologues are embraced. However, we can incorporate the full variety of triatomic molecules through a symmetry reduction. In this case, it is more convenient to start with an algebraic representation of the Hamiltonian in terms of local operators.
H ^ L = ћ ω 1 2 ( a ^ 1 a ^ 1 + a ^ 1 a ^ 1 ) + α β ( a ^ 2 a ^ 2 + a ^ 2 a ^ 2 ) λ 1 ( 1 ) ( a ^ 1 a ^ 2 + a ^ 1 a ^ 2 ) + λ 2 ( 1 ) ( a ^ 1 a ^ 2 + a ^ 1 a ^ 2 ) ,
with the new definitions
ω 1 = f 11 g 11 ( 0 ) ; λ 1 ( 1 ) = α β x f + x g ; λ 2 ( 1 ) = α β x f x g
and
x f = f 12 f 11 f 22 ; x g = g 12 ( 0 ) g 11 ( 0 ) g 22 ( 0 ) ; α 2 = g 22 ( 0 ) g 11 ( 0 ) ; β 2 = f 22 f 11 .
The algebraic Hamiltonian (33) is obtained from the quadratic Hamiltonian for two interacting oscillators in configuration space. For the general case like FCN, α 1 and β 1 . Assuming the Born–Oppenheimer approximation, for asymmetric isotopologues of type x 1 B y A x 2 B , we have β = 1 , while for symmetric molecules, it is clear that α = β = 1 ; the case previously analyzed.
We proceed to obtain the normal representation of the Hamiltonian (33). To consider this scenario, we invoke the GF formalism to obtain the normal modes. The normal modes are defined in terms of the internal coordinates using the transformation L 1 G 0 F L = Λ , where the matrices G 0 and F are chosen in terms of the symmetry-adapted coordinates S = { S g = ( q 1 + q 2 ) / 2 , S u = ( q 1 q 2 ) / 2 } [33]. Here, the condition L G 0 1 L = 1 provides the normalization of L .
In the normal mode scheme, the algebraic Hamiltonian equivalent to (2) takes the general form
H ^ N = ι ћ Ω ι ( Γ ) 2 γ ( A ^ ι γ A ^ ι γ + A ^ ι γ A ^ ι γ ) ,
where we have taken into account that the ι -th normal mode carries the Γ -th irreducible representation with components γ . For the particular case of two non-equivalent interacting oscillators, the Hamiltonian (36) simplifies to
H ^ N = ι ћ Ω ι ( Γ ) 2 ( A ^ ι A ^ ι + A ^ ι A ^ ι ) ,
with Γ identified with the totally symmetric irreducible representation. It was found that the connection between the bosonic normal operators A ^ ι ( A ^ ι ) associated with the ι -th normal mode and the local operators a ^ j ( a ^ j ) is the following [33]
A ^ ι = j = 1 2 ( c + ( ι , j ) a ^ j + c ( ι , j ) a ^ j ) ,
with coefficients
c ± ( ι , j ) = 1 2 2 σ ( 1 ) ( σ + 1 ) ( 1 δ j 1 ) β ι ћ ω j μ j ( L 1 ) ι σ ± ( L ˜ ) ι σ ω j μ j ћ β ι 2 .
where β ι 2 = 1 ћ ( λ ι ) 1 / 2 , with λ ι = Λ ι ι given by
λ ± = g 11 f 11 2 1 + η 2 + 2 η x f x g ± 1 + η 2 ( η 2 2 ) + 4 η [ η ( x f 2 + x g 2 ) + ( 1 + η 2 ) x f x g ] ,
with η = α β . In Equation (39), σ denotes a sum over the symmetry-adapted coordinates σ = { 1 g , 2 u } . The matrix L is defined through S = L Q :
Q ι = α ( L 1 ) ι α S α ; α = 1 , 2 .
If we now into account that μ j = 1 / g q j q j ( 0 ) and ω j = f q j q j g q j q j ( 0 ) , Equation (38) allows the parameter δ to be defined
δ = 1 2 ι = 1 2 j = 1 2 | c ( ι , j ) | 2 ,
which measures the LN degree and also the feasibility of establishing a polyad preserving canonical transformation from (38). Here, we have to stress that, because of the non-equivalence between the oscillators, the coefficients c ± ( ι j ) cannot be factorized as happened in (4).
In order to identify the parameters (19), we recall the connection between (32) and λ 2 ( 1 ) in (27). Hence, from λ 2 ( 12 ) λ 2 ( 1 ) , we define
γ ( ω ) = 1 8 ( x f x g ) 2 ,
as the parameter connected with the Hessian in (18), albeit with the new definitions (35).
As we know, the substitution of (38) into the Hamiltonian (2) leads to the algebraic Hamiltonian (33) given in the local representation. But now, in order to obtain the equivalent expression to (28), we need to consider the polyad-conserving canonical transformations
A ^ 1 = α c ^ 1 + 1 α 2 c ^ 2
A ^ 2 = 1 α 2 c ^ 1 α c ^ 2
to ensure the commutation relations [ A ^ ι , A ^ ι ] = δ ι ι , where α is chosen to be α = c + ( 1 , 1 ) . Substitution of (44) into (2) gives rise to the Hamiltonian
H ^ L ( P L ) = ω n o r ( 1 ) ( c ^ 1 c ^ 1 + c ^ 1 c ^ 1 ) + ω n o r ( 2 ) ( c ^ 2 c ^ 2 + c ^ 2 c ^ 2 ) + λ n o r ( 12 ) ( c ^ 1 c ^ 2 + c ^ 1 c ^ 2 ) ,
where
ω n o r ( 1 ) = ћ 2 ( Ω 1 α 2 + Ω 2 ( 1 α 2 ) )
ω n o r ( 2 ) = ћ 2 ( Ω 1 ( 1 α 2 ) + Ω 2 α 2 )
λ n o r ( 12 ) = ћ α 1 α 2 ( Ω 1 Ω 2 ) .
Both the Hamiltonians (33) and (45) involve three force constants { f 11 , f 22 , f 12 } . This is indeed the case in a general situation like the FCN molecule, where the two fundamentals are not enough to determine them. In such situations, fundamentals for an isotopologue are needed to estimate the three force constants. For symmetric isotopologues f 11 = f 22 , and just the two fundamentals need to be determined.
The transformation (44) acquires a preponderant importance in describing molecules with clear normal mode behavior, like CO 2 , N 2 O, OCS, or FCN, for instance. In molecules with local mode behavior, the usual approach consists in describing the stretching internal coordinates in terms of interacting Morse oscillators, with the bends in terms of normal coordinates. The splitting of the stretching fundamental energies is small and consequently a local polyad-conserving Hamiltonian may be a good approximation. In contrast, for the bending degrees of freedom , this is not possible because the polyad is only well defined in the normal scheme. The same situation is present for the stretching degrees of freedom in molecules with strong normal behavior. In order to embrace both behaviors, we propose to apply the transformation (44) to a normal polyad Hamiltonian to obtain a local representation. The obtained Hamiltonian is later anharmonized by mapping the bosonic operators c ^ i ( c ^ i ) for stretches to S U ( 2 ) operators b ^ i ( b ^ i ) associated with Morse ladder operators. This approach allows us to take advantage of the Morse properties in systems where only the normal polyad is defined. This procedure has been tested on the whole series of isotopologues of carbon dioxide [33,53] and the FCN molecule [54].

3.3. Parameterization Local–Normal Mode Transition

Because of the close relation between the contribution of the annihilation operators in (38) and local polyad breaking P ^ L , it is convenient to recall the parameterization from H2O to CO 2 presented in Refs. [30,31], but now taking care of the linear path appearing in the diagram δ vs. γ . Both parameters depend on x g and x f . Hence, taking t x g , we have for x f ( t )
x f ( t ) = m f ( t x g ( N ) ) + x f ( N ) ; m f = x f ( N ) x f ( L ) x g ( N ) x g ( L ) ,
ω ( t ) = m ω ( t x g ( N ) ) + ω ( N ) ; m ω = ω ( N ) ω ( L ) x g ( N ) x g ( L ) ,
with t [ x g ( L ) , x g ( N ) ] , while ( x g ( L ) , x f ( L ) ) and ( x g ( N ) , x f ( N ) ) correspond to the values of the parameters for the molecules H 2 O and CO 2 , respectively. This parametrization can be used to study Equation (22) connecting the normal and local polyads. From Figure 3, we can see that, at the local limit, all parameters vanish, with the exception of β 0 going to unity and leading to P ^ N = P ^ L , as expected. We found that these parameters are basically correlated to the δ parameter and consequently they do not contribute to additional LN criteria.
Up to now, we have considered the correlation of δ vs. γ ( ω ) for molecules near a local mode behavior. Figure 1 suggests that the linear trend is a consequence of the local mode behavior, since the location of CO 2 and CS 2 lies outside of the linear correlation. This situation is clearly manifested in Figure 4, left side, where the series of molecules BeX 2 and MgX 2 are analyzed. We can see that by including molecules with strong normal mode behavior the function changes from linear to exponential form, making it clear that the linear correlation is just an approximation. To explain this behavior, we notice that from the point of view of the series (18), γ ( ω ) represents the first correction to the linear approximation. When higher-order terms are incorporated, the linear trend tends to be recovered, as seen in the same Figure 4 on the right-hand side. This behavior also explains the location of the molecules CO 2 and CS 2 above the line in Figure 1. On the other hand, the parameter δ is unique, while there are several γ s . This fact suggests that putting together the information of all the γ s , a closer approximation to δ is obtained. This guess is confirmed by the plot of δ vs. γ 10 ( ω ) + γ 9 ( i j ) depicted on the right of Figure 4. The subscript index in the γ ’s means the upper order is taken into account. In our analysis, the addition of higher-order terms to γ ( ω ) is not necessary to define a different parameter, because of the clear exponential form, which ensures a consistent behavior, and consequently both δ and γ ( ω ) provide LN parameters, even in extreme conditions of normality.

3.4. Isotopologues

The analysis of the series of isotopologues deserves special attention, because the plots δ vs. γ ( ω ) provide a perfect linear correlation, which is expected from Figure 4, since they represent a short segment of the exponential function. In addition, since each series of isotopologues is characterized using the same force constants, the slopes are expected to be correlated to x f .
The case of equivalent oscillators embraces the isotopologues of type x A y B 2 . Examples of molecules with sufficient experimental information include H 2 O , NO 2 , SO 2 , O 3 , and CO 2 . In Figure 5, the plots δ vs. γ are displayed for H 2 O , SO 2 , and O 3 . For the isotopologues of carbon dioxide, the corresponding plot is given in Figure 1 of Ref. [33]. For NO 2 , only two isotopologues are given, and no plot is necessary. In Figure 6, the force constants vs. the slope m are depicted. The molecules H 2 O , SO 2 , and CO 2 are in close agreement with the expected linear trend, although for SO 2 and CO 2 , we have corrected the slope, adding higher order terms to γ ( ω ) , as explained before. However, two molecules, NO 2 and O 3 , present a quite different behavior, which we believe is due to their resonant structure manifested by the different values of force constants.
In this analysis, our results are obtained either using force constants from ab initio calculations or from the fundamentals, but in either case the force constants must be the same for every isotopologue.

4. Pyramidal Molecules

In Ref. [20], the case of three equivalent oscillators was analyzed in the context of the pyramidal molecules for both stretching and bending modes, while in Ref. [61] a study of the stretching modes in the molecule BF 3 was considered. Here, we focus on pyramidal molecules, because of their abundance.
Pyramidal molecules present two stretching and two bending normal modes, both with symmetries A 1 E . The fundamental energies of the bending modes are far from the stretching frequencies, hence a good approximation consists in neglecting stretching–bending interactions up to second order, where the normal modes are defined. This means that the stretches and bends can be independently analyzed in terms of a Hamiltonian of type (2).
We start with the stretching degrees of freedom. The expressions for the force constants (3) are [20]
F A 1 A 1 = f r r + 2 f r r ; G A 1 A 1 ( 0 ) = g r r ( 0 ) + 2 g r r ( 0 ) ,
F E E = f r r f r r ; G E E ( 0 ) = g r r ( 0 ) g r r ( 0 ) ,
while for the matrix elements involved in (6)
| | m i , Γ γ | | = 1 3 2 6 0 1 3 1 6 1 2 1 3 1 6 1 2 ,
which was chosen to be associated with the group chain C 3 v C s a with C s a = { E , σ v a } , with the same notation for the symmetry elements used in Ref. [20]. In this case, the algebraic Hamiltonian (8) takes the form
H ^ L = λ 0 i = 1 3 ( a ^ i a ^ i + a ^ i a ^ i ) + λ 1 ( 1 ) i > j = 1 3 ( a ^ i a ^ j + a ^ i a ^ j ) + λ 2 ( 1 ) i > j = 1 3 ( a ^ i a ^ j + a ^ i a ^ j ) ,
with coefficients given by (11). This Hamiltonian does not conserve the local polyad. Again, if we demand λ 2 ( 1 ) = 0 , we are able to estimate the force constants using (11), as well as the matrix representation of the Hamiltonian in the local basis L = { | 100 , | 010 , | 001 } . The results were extracted from Table 1 with μ 2 = E , allowing the parameters ϵ s to be calculated.
To obtain a polyad-conserving Hamiltonian, we should consider the canonical transformation (15), which when substituted into the Hamiltonian (2) yields
H ^ L ( P L ) = ω n o r i = 1 3 ( c ^ i c ^ i + c ^ i c ^ i ) + λ n o r ( 1 ) i > j = 1 3 ( c ^ i c ^ j + c ^ i c ^ j ) ,
where the spectroscopic parameters given by
ω n o r = ћ ω 0 2 1 3 ( 1 + 2 x f ) ( 1 + 2 x g ) + 2 3 ( 1 x f ) ( 1 x g ) ,
λ n o r ( 1 ) = ћ ω 0 2 2 3 ( 1 + 2 x f ) ( 1 + 2 x g ) 2 3 ( 1 x f ) ( 1 x g ) .
are functions of x , in accordance with (18). Their expansion, equivalent to (18), leads to the identification
γ ( 12 ) γ ( 1 ) = γ ( ω ) = 1 4 ( x f x g ) 2 ,
which again turns out to be basically the Hessian of the Taylor expansions of (52). In Table 4, the fundamentals as well as the calculated force constants are given for several molecules. Because it is possible to establish an isomorphism between the stretches and the bends in such a way that both span the same irreducible representations, the general results are basically valid, with the proviso that the structure constants change for the bends [20]. In Table 5, the LN parameters for the stretches are presented, while in Table 6 the corresponding parameters for the bends are listed.
In Figure 7, a plot of δ vs. γ ( ω ) is displayed for the stretching modes for several molecules, including the pyramidal molecules analyzed in Ref. [20]. Again, a linear trend is manifested, because all the molecules are in the local mode region (small values of δ ). In Figure 8, the plot δ vs. γ ( ω ) for the bending modes is displayed. One point to stress is that the NH 3 molecule is the only one located out of the line, a fact we assumed to be due to the existence of the inversion mode [20]. We have thus identified two cases where the deviation of the linear behavior allows particular molecular properties to be identified: internal inversion and resonance.
Although we do not present the plot ϵ 1 vs. δ , it turned out to be similar for both stretches and bends, a behavior that can be extracted from Table 5 and Table 6. Regarding ϵ 2 , for the stretches, the expected correlation with δ is fulfilled, while this is not the case for the bends. From the same Tables, it is clear that the parameter ζ associated with the splitting of labels does not represent an LN parameter.

5. Tetrahedral Molecules

Here, we present a study of tetrahedral molecules for both stretching and bending degrees of freedom. This system is relevant because the bending modes present a spurious state. First, we consider the four stretching oscillators, which reduce to the A 1 F 2 irreducible representations. Again, the fundamental energies of the bending modes are far from the stretching frequencies, and consequently to obtain the normal modes, a good approximation consists in neglecting stretching–bending interactions up to second order.

5.1. Stretching Oscillators

The Hamiltonian in the normal scheme for the stretching modes takes the form (2), with
F A 1 A 1 = f r r + 3 f r r ; G A 1 A 1 ( 0 ) = g r r ( 0 ) + 3 g r r ( 0 ) ,
F F 2 F 2 = f r r f r r ; G F 2 F 2 ( 0 ) = g r r ( 0 ) g r r ( 0 ) .
From these expressions and the fundamentals, we obtained the estimation of the force constants (3) given in Table 1 with μ 2 = F 2 . Here, the components of F 2 are labeled in accordance with the canonical chain T d C 2 v . The connection between the normal bosonic operators and the local operators is given by (4), with δ ± Γ obtained through (5). The matrix elements | | m i , Γ γ | | define the symmetry adapted combinations for one quantum
| | m i , Γ γ | | = 1 2 1 2 1 2 0 1 2 1 2 0 1 2 1 2 1 2 1 2 0 1 2 1 2 0 1 2
with order { A 1 , ( F 2 , A 1 ) , ( F 2 , B 2 ) , ( F 2 , B 3 ) } .
The substitution of (4) into the Hamiltonian (2) leads to the algebraic representation of the Hamiltonian in the local scheme
H ^ L = λ 0 i = 1 4 ( a ^ i a ^ i + a ^ i a ^ i ) + λ 1 ( 1 ) i > j = 1 4 ( a ^ i a ^ j + a ^ i a ^ j ) + λ 2 ( 1 ) i > j = 1 4 ( a ^ i a ^ j + a ^ i a ^ j ) ,
with coefficients defined by (11) with the convention λ ( 12 ) λ ( 1 ) , since only one type of interaction is present. As expected, this Hamiltonian does not conserve the local polyad. Demanding λ 2 ( 1 ) = 0 in (56) and taking its diagonalization in the space of one local quanta L = { | 1000 , | 0100 , | 0010 , | 0001 } , we are able to estimate the force constants f r r ( L ) , f r r ( L ) from Table 1.
Let us now consider the equivalent expression (16). The normal operators neglecting the contribution of the annihilation local operators are given by (15). Their substitution into the Hamiltonian (2) leads to
H ^ L ( P L ) = ω n o r i = 1 4 ( c ^ i c ^ i + c ^ i c ^ i ) + λ n o r ( 1 ) i > j = 1 4 ( c ^ i c ^ j + c ^ i c ^ j ) ,
with coefficients
ω n o r = ћ ω 0 2 1 4 ( 1 + x f ) ( 1 + 3 x g ) + 3 4 ( 1 x f ) ( 1 x g ) ,
λ n o r ( 1 ) = ћ ω 0 2 1 2 ( 1 + 3 x f ) ( 1 + 3 x g ) 1 2 ( 1 x f ) ( 1 x g ) ,
which when expanded in terms of the x f and x g  allow the identification
γ ( ω ) = 3 8 ( x f x g ) 2 ;
γ ( 12 ) γ ( 1 ) = 1 2 ( x f x g ) 2 .
Since these expressions are proportional, we introduce the unique parameter
γ ( x f x g ) 2 .
The fits of the parameters involved in the Hamiltonians (56) and (57) to reproduce the fundamentals. Both fits provide the same numerical values for the sets { λ 0 , λ 1 ( 1 ) } and { ω n o r , λ n o r ( 1 ) } , but the connections to the force constants are different. From these results the parameters ϵ s  can be obtained.
In Table 7, the fundamentals and the force constants calculated from Table 1 are given, while in Table 8 the LN parameters can be found. Using these results in Figure 9 the plot of δ vs. γ is depicted, including several tetrahedral molecules for the stretching modes. On the right side of the figure, molecules with extreme normal mode behavior are displayed, while on the left, molecules with local mode behavior are included. The exponential trend in the former, behavior explained in Section 3.3.

5.2. Bending Oscillators

The set of bending oscillators in tetrahedral molecules is interesting because of the presence of a spurious state, a common situation in molecules with high symmetry in the framework of internal coordinates. The subspace of six oscillators spans the irreducible representations A 1 E F 2 , with A 1 identified as a spurious mode. Consequently, the Hamiltonian (2) only involves the E F 2 modes with
Ω E = ( f θ θ 2 f θ θ + f θ θ ) ( g θ θ ( 0 ) 2 g θ θ ( 0 ) + g θ θ ( 0 ) )
Ω F 2 = ( f θ θ f θ θ ) ( g θ θ ( 0 ) g θ θ ( 0 ) ) .
Here, we have three force constants and two fundamentals with the constraint
Ω A 1 = ( f θ θ + 4 f θ θ + f θ θ ) ( g θ θ ( 0 ) + 4 g θ θ ( 0 ) + g θ θ ( 0 ) ) = 0 .
The estimation of the force constants are given in Table 1, with μ 2 = E and μ 3 = F 2 , omitting from the sums the spurious state A 1 . The normal and the local operators are related through the transformation (4). The symmetry projection matrix takes the form
| | m i , Γ γ | | = 1 6 1 3 0 1 2 0 0 1 6 1 2 3 1 2 0 1 2 1 2 1 6 1 2 3 1 2 0 1 2 1 2 1 6 1 3 0 1 2 0 0 1 6 1 2 3 1 2 0 1 2 1 2 1 6 1 2 3 1 2 0 1 2 1 2
Here, the projection is associated with the group chain T d C 2 v with the following order for the irreducible representations { A 1 , ( E , A 1 ) , ( E , A 2 ) , ( F 2 , A 1 ) , ( F 2 , B 1 ) , ( F 2 , B 2 ) } .
Because of the presence of the spurious mode, this system will be analyzed following a different route. The substitution of the transformations (4) into (2), without eliminating the spurious state A 1 , leads to the local representation
H ^ L = λ 0 i a ^ i a ^ i + a ^ i a ^ i + λ 1 ( 1 ) j > i a ^ i a ^ j + a ^ i a ^ j + λ 1 ( 2 ) j > i a ^ i a ^ j a ^ i a ^ j + λ 2 ( 1 ) j > i a ^ i a ^ j + a ^ i a ^ j + λ 2 ( 2 ) j > i a ^ i a ^ j + a ^ i a ^ j ,
with coefficients
λ 0 = ћ Ω A 1 12 δ + , A 1 + δ , A 1 + Ω E 6 δ + , E + δ , E + Ω F 2 4 δ + , F 2 + δ , F 2 ;
λ 1 ( 1 ) = ћ Ω A 1 12 δ + , A 1 + δ , A 1 Ω E 12 δ + , E + δ , E ;
λ 1 ( 2 ) = ћ Ω A 1 12 δ + , A 1 + δ , A 1 + Ω E 12 δ + , E + δ , E Ω F 2 8 δ + , F 2 + δ , F 2 ;
λ 2 ( 1 ) = ћ Ω A 1 12 δ + , A 1 δ , A 1 Ω E 12 δ + , E δ , E ;
λ 2 ( 2 ) = ћ Ω A 1 12 δ + , A 1 δ , A 1 + 2 Ω E 12 δ + , E δ , E Ω F 2 4 δ + , F 2 δ , F 2 ,
where here we notice that two types of pairwise interactions appear. We decided to write down the coefficients in this form, because this makes evident that when δ Γ = 0 , we obtain the expected results (11). However, when the A 1 mode is eliminated, the coefficients become
λ 0 = ω 0 24 [ 10 + 4 ( x f + x g ) + ( x f + x g ) ]
λ 1 ( 1 ) = ω 0 24 [ 2 2 ( x f + x g ) + ( x f + x g ) ] ,
λ 1 ( 2 ) = ω 0 24 [ 2 + 4 ( x f + x g ) 5 ( x f + x g ) ] ,
λ 2 ( 1 ) = ω 0 24 [ 2 ( x f x g ) + ( x f x g ) ] ,
λ 2 ( 2 ) = ω 0 24 [ 4 ( x f x g ) 5 ( x f x g ) ] .
As noticed, the general form (11) has been broken. But we will not deal with the expressions (66) used to calculate the force constants. Instead, we return to (11), as we next explain. The simplest way to proceed consists in starting with the Hamiltonian (64), but neglecting the non-conserving polyad contributions:
H ^ L ( P L ) = λ 0 i ( a ^ i a ^ i + a ^ i a ^ i ) + λ 1 ( 1 ) i j ( a ^ i a ^ j + a ^ i a ^ j ) + λ 1 ( 2 ) i j ( a ^ i a ^ j + a ^ i a ^ j ) ,
with the identification (11). The diagonalization of this Hamiltonian in one quantum local basis L 6 = { | 100000 , | 010000 , | 001000 , | 000100 , | 000010 , | 000001 , } leads to three eigenvalues, which when they are identified with the fundamentals ћ Ω Γ lead to
ћ Ω A 1 = ћ ω 0 1 + 2 ( x f + x g ) + 1 2 ( x f + x g )
ћ Ω E = ћ ω 0 1 ( x f + x g ) + 1 2 ( x f + x g )
ћ Ω F 2 = ћ ω 0 1 1 2 ( x f + x g ) ,
When we impose the condition (62) with definitions (12), we are able to obtain the local force constants f q i q j ( L )  from the equations given in Table 1.
We now proceed to obtain the γ s parameters. Applying the transformation (15) to the Hamiltonian (2), preserving the spurious state, we obtain
H ^ L ( P L ) = ω n o r i ( c ^ i c ^ i + c ^ i c ^ i ) + λ n o r ( 1 ) i j ( c ^ i c ^ j + c ^ i c ^ j ) + λ n o r ( 2 ) i j ( c ^ i c ^ j + c ^ i c ^ j ) ,
where
ω n o r = ћ Ω A 1 12 + Ω E 6 + Ω F 2 4 ;
λ n o r ( 1 ) = ћ Ω A 1 12 Ω E 12 ;
λ n o r ( 2 ) = ћ Ω A 1 12 + Ω E 6 Ω F 2 4 ,
The Equations (70) are implicit functions of the variables x . Their expansion leads to
ω n o r = ћ ω 0 1 2 γ ( ω ) + O ( x 3 ) ;
λ n o r ( 1 ) = ћ ω 0 1 4 ( x f + x g ) γ ( 1 ) + O ( x 3 ) ;
λ n o r ( 2 ) = ћ ω 0 1 4 ( x f + x g ) γ ( 2 ) + O ( x 3 ) ,
with the identification
γ ( ω ) = 1 4 ( x f x g ) 2 + 1 16 ( x f x g ) 2 ;
γ ( 1 ) = 1 8 ( x f x g ) ( x f x g ) + ( x f x g ) ;
γ ( 2 ) = 1 4 ( x f x g ) 2 .
We stress that in these expressions the spurious state A 1 has been included. However, the redundancy takes into account the calculation of force constants f q i q j ( N ) through the use of (61) together with the constraint (62). In Table 9, the LN parameters are listed, while in Figure 10, the plot δ vs. γ ( ω ) is displayed. In addition in Figure 11, the corresponding plots involving the parameters γ ( 1 ) and γ ( 2 ) are presented. These results show the consistency of the parameters in establishing an LN degree, even in the presence of spurious states.

5.3. Isotopologues

Let us now consider the series of isotopologues α X β H 4 , X = Ge, Si, N, C; involving the bending degrees of freedom. The plots δ vs. γ ( ω ) are displayed in Figure 12. Since similar results were obtained for δ vs. γ ( 1 ) and δ vs. γ ( 2 ) , their corresponding plots are not included. The linear correlations are evident, as in the previous cases. This is clearly explained by the local mode behavior of the systems. However, this is not an obvious result, since we are dealing with bending degrees of freedom presenting a spurious state. We again stress that bending modes are not traditionally contemplated when assigning an LN degree, but with these results we confirm that the parameters we have introduced represent a valid measure of the LN degree.
From Figure 12, we may also consider the slopes to look for a correlation with the force constants x f . The result is presented in Figure 13. In order to create the plot, we had to take into account that several force constants are available. This fact is taken into account by the bars in the plot. We can see that with the exception of GeH 4 , a line can be assigned for the correlation. We believe that the odd behavior of GeH 4 with respect to the other compounds is due to the Germanium configuration [ A r ] 4 s 2 3 d 10 4 p 2 , in contrast to C, N, and Si, where no d orbitals are present. A similar situation appears in the series of pyramidal molecules, where the unexpected local-to-normal order is believed to have the same origin [20].

6. Octahedral Molecules

We now consider the stretching modes of octahedral molecules, with the O h symmetry group. This system is included in our study because of the presence of two types of interactions: contiguous and opposite ones, but also because, in the case of contiguous bonds, the parameter x g vanishes. A partial analysis of this system, including vibrational descriptions of the series of molecules SF 6 , WF 6 , and UF 6 , was discussed in Ref. [32].
The reduction of the internal coordinates corresponds to A 1 g E g T 1 u . In this case, the force and structure constants take the form
F A 1 g A 1 g = f r r + 4 f r r + f r r ; G A 1 g A 1 g ( 0 ) = g r r ( 0 ) + 4 g r r ( 0 ) + g r r ( 0 ) ;
F E g E g = f r r 2 f r r + f r r ; G E g E g ( 0 ) = g r r ( 0 ) 2 g r r ( 0 ) + g r r ( 0 ) ;
F T 1 u T 1 u = f r r f r r ; G T 1 u T 1 u ( 0 ) = g r r ( 0 ) g r r ( 0 ) .
from these expressions and the fundamental energies, the force constants can be obtained from Table 1, with μ 3 = T 1 u . The results are displayed in Table 10. The matrix elements | | m i , Γ γ | | defining the symmetry-adapted combinations for one quantum are given by
| | m i , Γ γ | | = 1 6 1 2 3 1 2 0 1 2 1 2 1 6 1 2 3 1 2 0 1 2 1 2 1 6 1 2 3 1 2 0 1 2 1 2 1 6 1 2 3 1 2 0 1 2 1 2 1 6 1 3 0 1 2 0 0 1 6 1 3 0 1 2 0 0 ,
in accordance with chain O h D 2 h . The substitution of (4) into the Hamiltonian (2) leads to the local representation
H ^ L = λ 0 i = 1 6 ( a ^ i a ^ i + a ^ i a ^ i ) + λ 1 ( 1 ) i > j ( a ^ i a ^ j + a ^ i a ^ j ) + λ 1 ( 2 ) i > j ( a ^ i a ^ j + a ^ i a ^ j ) + λ 2 ( 1 ) i > j ( a ^ i a ^ j a ^ i a ^ j ) λ 2 ( 2 ) i > j ( a ^ i a ^ j + a ^ i a ^ j ) ,
with coefficients provided by (11). A polyad-conserving Hamiltonian is obtained by setting λ 2 ( 1 ) = λ 2 ( 2 ) = 0 . The diagonalization of the Hamiltonian in the one-quantum basis L 6 = { | 100000 , | 010000 , | 001000 , | 000100 , | 000010 , | 000001 } allows us to express the spectroscopic parameters in terms of the fundamental energies, from which we estimate the force constants f q i q j ( L ) . The results are displayed in Table 1. We may now apply the polyad-conserving transformation (15) to the Hamiltonian (2) to obtain
H ^ L ( P L ) = ω n o r i = 1 6 ( c ^ i c ^ i + c ^ i c ^ i ) + λ n o r ( 1 ) i > j ( c ^ i c ^ j + c ^ i c ^ j ) + λ n o r ( 2 ) i > j ( c ^ i c ^ j + c ^ i c ^ j ) ,
with spectroscopic coefficients given by
ω n o r = ћ 2 1 6 Ω A 1 g + 1 3 Ω E g + 1 2 Ω T 1 u ,
λ n o r ( 1 ) = ћ 2 1 3 Ω A 1 g 1 3 Ω E g ,
λ n o r ( 2 ) = ћ 2 1 3 Ω A 1 g + 2 3 Ω E g Ω T 1 u ,
with
Ω A 1 g = ω 0 ( 1 + 4 x f + x f ) ( 1 + 4 x g + x g )
Ω E g = ω 0 ( 1 2 x f + x f ) ( 1 2 x g + x g )
Ω T 1 u = ω 0 ( 1 x f ) ( 1 x g ) .
In this case, the expansion in terms of x leads to the identification
γ ( ω ) = 1 2 ( x f x g ) 2 + 1 8 ( x f x g ) 2 ;
γ ( 1 ) = 1 2 ( x f x g ) ( x f x g ) + ( x f x g ) ;
γ ( 2 ) = ( x f x g ) 2 .
In Table 11, the corresponding parameters associated with the LN degree are presented. Based on these results, in Figure 14 the plot δ vs. γ ( ω ) is displayed, obtaining a clear linear correlation. In addition, in Figure 15, the plots involving ϵ 1 and ϵ 2 are included. As expected, the ϵ 1 shows a clear linear correlation, while for ϵ 2 , although a linear trend is obtained, it is not as clear as for ϵ 1 . The parameter ϵ 3 was not included, because no correlation appeared. Again, from Table 11, it is clear that the parameter ζ does not provide a general LN parameter.

7. Normal/Local Degree and Physical Properties

The importance of having a parameter measuring the local/normal degree is that it may be correlated with physical properties, depending on the vibrational degrees of freedom. First, we shall consider the effect on the partition function.

7.1. Partition Function

As a first case, we shall consider a molecular system in the framework of the Born–Oppenheimer approximation and rotor rigid approximation. If we focus on the stretching vibrational degrees of freedom and consider the Hamiltonian of two equivalent interacting oscillators reduced to the form (2), the partition function Z N takes the simple form
Z N = Γ = g , u e ћ Ω Γ 2 K T 1 e ћ Ω Γ K T ,
with Ω g = ω 0 ( 1 + x f ) ( 1 + x g ) and Ω u = ω 0 ( 1 x f ) ( 1 x g ) . Considering the parameterization (47), which is equivalent to moving along the parameter δ ( t ) and frequencies Ω Γ ( t ) , we are able to see the functional form Z N ( t ) , and consequently the properties depending on it. In Figure 16, left side, the change in the function Z N ( t ) along the LN parameter is displayed for different temperatures. We can see that in the local limit the partition function is close to zero and increases as the normal character becomes stronger. On the other hand, for two equivalent oscillators of frequency ω 0 without interaction, the partition function takes the form
Z L = e ћ ω 0 2 K T 1 e ћ ω 0 K T 2 ,
with ω 0 = f r r g r r ( 0 ) . It is clear that in the local limit both coincide: lim t 0 Z N = Z L . This behavior is shown on the right side of Figure 16. The importance of this analysis is that the dependence on δ of the partition function implies a a correlation with thermodynamic properties.

7.2. Spectroscopic Properties

We now study the correlation between the parameter δ and the wave functions obtained using an algebraic model based on s u ( 2 ) and s u ( 3 ) algebras [20,33,53,54]. First, we consider the series of pyramidal molecules PH 3 , AsH 3 , and SbH 3 , which present a local mode behavior. Thereafter, we analyze the isotopologues of carbon dioxide and cyanogen fluoride. The analysis of these systems permits proving that the parameters introduced are correlated with the spectroscopic properties of molecules.

7.2.1. Pyramidal Molecules

Let us start by considering the effect of the LN degree on the wave functions and consequently on the transition intensities. In Ref. [20], a spectroscopic description of the molecules PH 3 , AsH 3 , and SbH 3 was given. These molecules present a local mode behavior, which is manifested by the locality of the states. In Table 12, the maximum local components for several states characterized by having large local components are displayed. We have chosen states with experimental energies common to the three molecules. The first three states mostly have a stretching character, while the last one has a bending contribution. In order to appreciate the dependence on the LN parameter, the δ values have been included. It is interesting to notice the local to normal sequence for the states | 200000 and | 400000 is given by SbH 3   PH 3   AsH 3 , in accordance with the δ parameter. In contrast, for the state | 200100 , the sequence changes to SbH 3   AsH 3   PH 3 . This is explained by the fact that the former set is associated with stretches and the latter to the bends [20]. Hence, we have a correlation between the wave functions and the LN parameter δ , a result that may also exist in transition intensities.

7.2.2. Isotopologues of CO 2

In Refs. [33,53], vibrational analyses were carried out of the series of isotopologues of carbon dioxide using a S U 1 ( 2 ) × S U ( 3 ) × S U 2 ( 2 ) algebraic model based on the anharmonization of local operators applied to normal operators. It is known that, in these isotopologues, the Fermi interactions dominates the spectrum, a fact manifested in the wave functions. In Figure 3 of Ref. [53], a plot of the Fermi interaction strength | α F | vs. δ is displayed, showing a clear correlation. This result suggests that the wave functions will also be correlated. This is indeed the case, as we show next.
In Raman spectroscopy, four of the most intense lines are due to the transitions listed in Table 13 [78]:
These transitions are basically determined by the coefficients a and a , which turn out to be dominant in the Raman spectrum. The question that arises is whether or not these coefficients are correlated with the parameter δ . In Table 14, the coefficients of the basis together with the parameter δ are displayed. From this table, a correlation is clearly evident and graphically shown in Figure 4 of Ref. [53]. The real importance of this correlation lies in the impact on the description of the Raman spectrum, which may be useful in the identification of isotopologues. Preliminary results show that the Raman transition intensities are indeed correlated with the Raman intensities.

7.2.3. Isotopologues of FCN

As a second example of the importance of the LN degree in vibrational spectroscopy, we consider the vibrational degrees of freedom of the FCN molecule. Given the presence of the resonances ω 1 2 ω 2 and 2 ω 1 ω 3 , the appropriate polyad is P N = 2 ν 1 + 2 ν 2 + 4 ν 3 . In Ref. [54], vibrational description of this molecule was carried out. Here, we consider the four sets of states corresponding to polyads P N = 2 , 3 , 4 , 5 associated with symmetries Σ + and Π ± :
Σ + ; P N = 2 : a | 10 0 0 + b | 02 0 0 ,
Σ + ; P N = 4 : a | 04 0 0 + b | 12 0 0 + c | 20 0 0 + d | 00 0 1 ,
Π ± ; P N = 3 : α | 03 1 0 + β | 11 1 0 ,
Π ± ; P N = 5 : α | 05 1 0 + β | 13 1 0 + γ | 21 1 0 + δ | 01 1 1 .
For an specific polyad and symmetry, there is a multiplet of states interacting by the two resonances. Hence, the states are characterized by pairs. In Table 15, the squares of the maximum components are displayed, indicating the states in resonance, together with the parameter δ . Again, the correlation between the components and the LN parameter is manifested: going from local to normal mode behavior with the decrement in components manifested. This correlation is explicitly presented in Ref. [54].

8. Conclusions

In this contribution, we have presented, for the first time, a consistent set of LN criteria that can be applied to any molecule. In contrast to the long established theory of local molecules, where the criterion of locality depends on a model of interacting Morse oscillators, our criteria are based on the analysis of normal modes. Choosing a selected set of molecules, we have shown that the proposed LN criteria can be applied to a great variety of situations, from local to normal extremes. Each parameter presents its own features. The parameter δ measures the degree of locality from the normal point of view, while the Hessian γ ( ω ) can be associated with a local perspective. In addition, ϵ 1 (and in some cases ϵ 2 too) offers a third parameter, and this is defined taking ingredients from both local and normal mode schemes through the force constants. Analyses of several representative systems were presented. First, the most simple system of two oscillators was presented, in order to include the case where the normal modes do not coincide with symmetry-adapted coordinates. The pyramidal molecules were included, in order to show that the proposed parameters can also be applied to the bending degrees of freedom. In addition, tetrahedral molecules were studied because of the presence of spurious modes in the bending modes. Finally, octahedral molecules were incorporated in our analysis because of the presence of two types of interactions involving vanishing contributions of the Wilson matrix. Although this set of studied systems may be considered relatively simple, it was chosen to include every possibility encountered in the framework of local coordinates, which allowed us to conclude the validity of our LN criteria.
The correlations between the different parameters in these systems were investigated to prove their consistency, but also to show that through such correlations it is possible to identify particular signatures of the molecules with just the knowledge of the fundamental energies:non-rigidity in NH 3 , resonance structures in O 3 and NO 2 , and change in electronic configuration in the pnictogen pyramidal hydrides and GeH 4 . The isotopologues displayed a perfect linear correlation for the plots δ vs. γ ( ω ) . This fact allowed us to establish a linear correlation between the slopes and the force constants.
An important result is that the LN degree is correlated with physical properties. This conclusion was obtained by considering the behavior of the partition function for two oscillators, taking advantage of the parameterization H 2 O → CO 2 . A similar situation appears in spectroscopy; a clear dependence of the wave functions with the LN parameters, together with their correlation with the interaction strengths allowed us to conclude the importance of the LN parameters in the analysis of Raman and infrared spectroscopy. This finding is particularly relevant for series of isotopologues, where preliminary results indicated the existence of a correlation between the LN parameters and the relative transition intensities in the Raman spectra of the isotopologues of CO 2 .
The present formalism is based on a harmonic oscillator model, and only the fundamentals were involved in both the determination of the LN parameters and their correlation. Consequently, the proposed criteria can be applied to any molecular system, without limitations. A remarkable result is that these parameters are correlated with spectroscopic properties, taking into account the full description of the systems, where all the relevant interactions are included: anharmonicities and resonances. On the other hand, we selected internal coordinates in our treatment because of the physical meaning of the force constants, but in practice Cartesian coordinates are more appropriate to generalize our approach. In the latter case, efficient programs to obtain the normal modes are available and a work in this direction is in progress.
Finally, we conclude from this work that the LN parameters, in particular the δ , is a descriptor of a molecule, which opens the possibility of also being used in machine learning algorithms.

Author Contributions

Conceptualization, R.L.; methodology, R.L. and E.S.; software, E.S. and O.G.-J.; validation, E.S. and O.G.-J.; formal analysis, E.S., O.G.-J. and R.L.; investigation, E.S. and O.G.-J.; data curation, E.S. and O.G.-J.; writing—original draft preparation, R.L.; writing—review and editing, E.S. and O.G.-J.; visualization, E.S. and O.G.-J.; supervision, R.L.; funding acquisition, R.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by DGAPA-UNAM, México, under project IN-212224.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Herzberg, G. Molecular Spectra and Molecular Structure II. Infrared and Raman Spectra; D. Van Nostrand Company: Princeton, NJ, USA, 1945. [Google Scholar]
  2. Wilson, E.B.; Decious, J.C.; Cross, P.C. Molecular Vibrations. The Theory of Infrared and Raman Spectra; Dover Publications: Mineola, NY, USA, 1955. [Google Scholar]
  3. Fermi, E. Über den Ramaneffekt des Kohlendioxyds. Z. FüR Phys. 1931, 71, 250–259. [Google Scholar] [CrossRef]
  4. Darling, B.T.; Dennison, D.M. The Water Vapor Molecule. Phys. Rev. 1940, 57, 128–139. [Google Scholar] [CrossRef]
  5. Hollas, J.M. High Resolution Spectroscopy; John Wiley & Sons: Hoboken, NJ, USA, 1982. [Google Scholar]
  6. Walther, H. Laser Spectroscopy of Atoms and Molecules; Springer: Berlin/Heidelberg, Germany, 1976. [Google Scholar]
  7. Quack, M. Spectra and Dynamics of Coupled Vibrations in Polyatomic Molecules. Annu. Rev. Phys. Chem. 1990, 41, 839–874. [Google Scholar] [CrossRef]
  8. Mecke, R. Valenz-und Deformationsschwingungen mehratomiger Moleküle. III. Z. FüR Phys. Chem. 1932, 17B, 1–20. [Google Scholar] [CrossRef]
  9. Mecke, R. Das Rotationsschwingungsspektrum des Wasserdampfes. I. Z. FüR Phys. 1933, 81, 313–331. [Google Scholar] [CrossRef]
  10. Timm, B.; Mecke, R. Quantitative Absorptionsmessungen an den CH-Oberschwingungen einfacher Kohlenwasserstoffe. Z. FüR Phys. 1935, 98, 363–381. [Google Scholar] [CrossRef]
  11. Rumpf, K.; Mecke, R. Über die Absorptionsspektren einiger Benzolderivate im nahen Ultraroten bei großen Schichtdicken. Z. FüR Phys. Chem. 1939, 44B, 299–312. [Google Scholar] [CrossRef]
  12. Halonen, L.; Child, M.S. Model stretching overtone eigenvalues for SF6, WF6, and UF6. J. Chem. Phys. 1983, 79, 559–570. [Google Scholar] [CrossRef]
  13. Child, M.S.; Halonen, L. Overtone Frequencies and Intensities in the Local Mode Picture; Wiley: Hoboken, NJ, USA, 1984; pp. 1–58. [Google Scholar] [CrossRef]
  14. Halonen, L. Local Mode Vibrations in Polyatomic Molecules; Wiley: Hoboken, NJ, USA, 1998; pp. 41–179. [Google Scholar] [CrossRef]
  15. Halonen, L. Highly Excited States and Local Modes; John Wiley & Sons: Hoboken, NJ, USA, 2000. [Google Scholar]
  16. Jensen, P. An introduction to the theory of local mode vibrations. Mol. Phys. 2000, 98, 1253–1285. [Google Scholar] [CrossRef]
  17. Child, M.S.; Lawton, R.T. Local and normal vibrational states: A harmonically coupled anharmonic-oscillator model. Faraday Discuss. Chem. Soc. 1981, 71, 273. [Google Scholar] [CrossRef]
  18. Mills, I.; Robiette, A. On the relationship of normal modes to local modes in molecular vibrations. Mol. Phys. 1985, 56, 743–765. [Google Scholar] [CrossRef]
  19. Halonen, L. Internal coordinate Hamiltonian model for Fermi resonances and local modes in methane. J. Chem. Phys. 1997, 106, 831–845. [Google Scholar] [CrossRef]
  20. Guzmán-Juárez, O.; Suárez, E.; Lemus, R. A study of the local to normal mode transition in pyramidal molecules and their vibrational description in terms of an algebraic model. J. Mol. Spectrosc. 2023, 393, 111775. [Google Scholar] [CrossRef]
  21. Sibert, E.L.; Hynes, J.T.; Reinhardt, W.P. Quantum mechanics of local mode ABA triatomic molecules. J. Chem. Phys. 1982, 77, 3595–3604. [Google Scholar] [CrossRef]
  22. Davis, M.J. Analysis of highly excited vibrational eigenstates. Int. Rev. Phys. Chem. 1995, 14, 15–66. [Google Scholar] [CrossRef]
  23. Heller, E.J.; Davis, M.J. Molecular overtone bandwidths from classical trajectories. J. Phys. Chem. 1980, 84, 1999–2001. [Google Scholar] [CrossRef]
  24. Kellman, M.E.; Tyng, V. The Dance of Molecules: New Dynamical Perspectives on Highly Excited Molecular Vibrations. Accounts Chem. Res. 2007, 40, 243–250. [Google Scholar] [CrossRef] [PubMed]
  25. Xiao, L.; Kellman, M.E. Unified semiclassical dynamics for molecular resonance spectra. J. Chem. Phys. 1989, 90, 6086–6098. [Google Scholar] [CrossRef]
  26. Schmid, G.; Coy, S.; Field, R.; Silbey, R. Duffing’s oscillator and the normal to local mode transition in AB2 triatomic molecules. Chem. Phys. Lett. 1994, 219, 331–338. [Google Scholar] [CrossRef]
  27. Kellman, M.E.; Lynch, E.D. Semiclassical phase space evolution of Fermi resonance spectra. J. Chem. Phys. 1988, 88, 2205–2215. [Google Scholar] [CrossRef]
  28. Sako, T.; Yamanouchi, K.; Iachello, F. Classical-limit analysis of the normal-to-local transition by the algebraic effective Hamiltonian approach. J. Chem. Phys. 2001, 114, 9441–9452. [Google Scholar] [CrossRef]
  29. Xiao, L.; Kellman, M.E. Catastrophe map classification of the generalized normal–local transition in Fermi resonance spectra. J. Chem. Phys. 1990, 93, 5805–5820. [Google Scholar] [CrossRef]
  30. Bermúdez-Montaña, M.; Lemus, R.; Castaños, O. Fidelity, entropy, and Poincaré sections as tools to study the polyad breaking phenomenon. EPL (Europhys. Lett.) 2016, 116, 13001. [Google Scholar] [CrossRef]
  31. Bermúdez-Montaña, M.; Lemus, R.; Castaños, O. Polyad breaking phenomenon associated with a local-to-normal mode transition and suitability to estimate force constants. Mol. Phys. 2017, 115, 3076–3103. [Google Scholar] [CrossRef]
  32. Coello, L.; Lemus, R. Vibrational description of the stretching modes of octahedral molecules: A local-to-normal mode transition criterion. Mol. Phys. 2021, 119, e1865582. [Google Scholar] [CrossRef]
  33. Bermúdez-Montaña, M.; Rodríguez-Arcos, M.; Carvajal, M.; Ostertag-Henning, C.; Lemus, R. A Spectroscopic Description of Asymmetric Isotopologues of CO2. J. Phys. Chem. A 2023, 127, 6357–6376. [Google Scholar] [CrossRef] [PubMed]
  34. Bunker, P.R.; Jensen, P. Molecular Synnetry and Spectroscopy; National Research Council of Canada: Ottawa, ON, Canada, 1998. [Google Scholar]
  35. Mills, I.; Mompean, F. Local modes and x, k relations in ethene and propadiene. Chem. Phys. Lett. 1986, 124, 425–431. [Google Scholar] [CrossRef]
  36. Valle, R.G.D. Local-mode to normal-mode hamiltonian transformation for X − H stretchings. Mol. Phys. 1988, 63, 611–621. [Google Scholar] [CrossRef]
  37. Lehmann, K.K. Beyond the xK relations. Mol. Phys. 1989, 66, 1129–1137. [Google Scholar] [CrossRef]
  38. Law, M.M.; Duncan, J.L. Anharmonic stretching vibrations expressed as local modes. Mol. Phys. 1998, 93, 809–819. [Google Scholar] [CrossRef]
  39. Law, M.M.; Duncan, J.L. Anharmonically-coupled local mode to normal mode Hamiltonian transformations: Beyond the x, K relations. Mol. Phys. 1998, 93, 821–830. [Google Scholar] [CrossRef]
  40. Gambi, A.; Valle, R.G.D. Local mode and normal mode models for molecules with two non-equivalent C–H bonds. Mol. Phys. 2007, 105, 1779–1787. [Google Scholar] [CrossRef]
  41. Lehmann, K.K. On the relation of Child and Lawton’s harmonically coupled anharmonic–oscillator model and Darling–Dennison coupling. J. Chem. Phys. 1983, 79, 1098. [Google Scholar] [CrossRef]
  42. Kellman, M.E. Algebraic resonance dynamics of the normal/local transition from experimental spectra of ABA triatomics. J. Chem. Phys. 1985, 83, 3843–3858. [Google Scholar] [CrossRef]
  43. Kjaergaard, H.G.; Henry, B.R.; Wei, H.; Lefebvre, S.; Carrington, T.; Mortensen, O.S.; Sage, M.L. Calculation of vibrational fundamental and overtone band intensities of H2O. J. Chem. Phys. 1994, 100, 6228–6239. [Google Scholar] [CrossRef]
  44. Naumenko, O.; Campargue, A. H232S: First Observation of the (70±,0) Local Mode Pair and Updated Global Effective Vibrational Hamiltonian. J. Mol. Spectrosc. 2001, 210, 224–232. [Google Scholar] [CrossRef]
  45. Flaud, J.; Camypeyret, C.; Burger, H.; Willner, H. High-Resolution Analysis of the ν2, 2ν2, ν1, and ν3 Bands of H 2 80 Se. J. Mol. Spectrosc. 1993, 161, 157–169. [Google Scholar] [CrossRef]
  46. Kauppi, E.; Halonen, L. A simple curvilinear internal coordinate model for vibrational energy levels. J. Chem. Phys. 1992, 96, 2933–2941. [Google Scholar] [CrossRef]
  47. Saarinen, T.; Kauppi, E.; Halonen, L. Potential energy surface for F2O. J. Mol. Spectrosc. 1990, 142, 175–181. [Google Scholar] [CrossRef]
  48. Aures, R.; Gericke, K.H.; Maul, C.; Trott-Kriegeskorte, G.; Kawasaki, M.; Nakano, Y. Excited state dynamics of Cl2O in the near ultraviolet. J. Chem. Phys. 2002, 117, 2141–2150. [Google Scholar] [CrossRef]
  49. Babikov, Y.L.; Mikhailenko, S.N.; Barbe, A.; Tyuterev, V.G. S&MPO—An information system for ozone spectroscopy on the WEB. J. Quant. Spectrosc. Radiat. Transf. 2014, 145, 169–196. [Google Scholar] [CrossRef]
  50. Delon, A.; Jost, R. Laser induced dispersed fluorescence spectra of jet cooled NO2: The complete set of vibrational levels up to 10000 cm−1 and the onset of the X ˜ 2 A1 A ˜ 2 B2 vibronic interaction. J. Chem. Phys. 1991, 95, 5686–5700. [Google Scholar] [CrossRef]
  51. Miller, C.E.; Brown, L.R. Near infrared spectroscopy of carbon dioxide I. 16O12C16O line positions. J. Mol. Spectrosc. 2004, 228, 329–354. [Google Scholar] [CrossRef]
  52. Suzuki, I. General Anharmonic Force Constants of Carbon Disulfide. Bull. Chem. Soc. Jpn. 1975, 48, 1685–1690. [Google Scholar] [CrossRef]
  53. Bermúdez-Montaña, M.; Rodríguez-Arcos, M.; Carvajal, M.; Ostertag-Henning, C.; Lemus, R. Algebraic vibrational description of the symmetric isotopologues of CO2: 13C16O2, 12C18O2, 12C17O2. Chem. Phys. 2022, 557, 111481. [Google Scholar] [CrossRef]
  54. Suárez, E.; Guzmán-Juárez, O.; Lemus, R. A general local algebraic approach for molecules with normal mode behavior: Application to FCN. Submitted to Vibrational Spectroscopy.
  55. Tennyson, J.; Bernath, P.F.; Brown, L.R.; Campargue, A.; Carleer, M.R.; Császár, A.G.; Gamache, R.R.; Hodges, J.T.; Jenouvrier, A.; Naumenko, O.V.; et al. IUPAC critical evaluation of the rotational–vibrational spectra of water vapor. Part I—Energy levels and transition wavenumbers for H 2 17 O and H 2 18 O. J. Quant. Spectrosc. Radiat. Transf. 2009, 110, 573–596. [Google Scholar] [CrossRef]
  56. Tennyson, J.; Bernath, P.F.; Brown, L.R.; Campargue, A.; Császár, A.G.; Daumont, L.; Gamache, R.R.; Hodges, J.T.; Naumenko, O.V.; Polyansky, O.L.; et al. IUPAC critical evaluation of the rotational–vibrational spectra of water vapor, Part III: Energy levels and transition wavenumbers for H 2 16 O. J. Quant. Spectrosc. Radiat. Transf. 2013, 117, 29–58. [Google Scholar] [CrossRef]
  57. Tennyson, J.; Bernath, P.F.; Brown, L.R.; Campargue, A.; Császár, A.G.; Daumont, L.; Gamache, R.R.; Hodges, J.T.; Naumenko, O.V.; Polyansky, O.L.; et al. IUPAC critical evaluation of the rotational–vibrational spectra of water vapor. Part IV. Energy levels and transition wavenumbers for D 2 16 O, D 2 17 O, and D 2 18 O. J. Quant. Spectrosc. Radiat. Transf. 2014, 142, 93–108. [Google Scholar] [CrossRef]
  58. Tóbiás, R.; Furtenbacher, T.; Császár, A.G.; Naumenko, O.V.; Tennyson, J.; Flaud, J.M.; Kumar, P.; Poirier, B. Critical evaluation of measured rotational–vibrational transitions of four sulphur isotopologues of S16O2. J. Quant. Spectrosc. Radiat. Transf. 2018, 208, 152–163. [Google Scholar] [CrossRef]
  59. Perrin, A.; Flaud, J.M.; Camy-Peyret, C.; Vasserot, A.M.; Guelachvili, G.; Goldman, A.; Murcray, F.; Blatherwick, R. The ν1, 2ν2, and ν3 interacting bands of 14N16O2: Line positions and intensities. J. Mol. Spectrosc. 1992, 154, 391–406. [Google Scholar] [CrossRef]
  60. Perrin, A.; Toon, G.; Orphal, J. Detection of atmospheric 15NO2 in the ν3 spectral region (6.3 µm). J. Quant. Spectrosc. Radiat. Transf. 2015, 154, 91–97. [Google Scholar] [CrossRef]
  61. Amezcua-Eccius, C.; Lemus, R. An approach for the description of vibrational excitations: Application to 11BF3. J. Mol. Spectrosc. 2010, 260, 36–49. [Google Scholar] [CrossRef]
  62. Rajamäki, T.; Miani, A.; Halonen, L. Six-dimensional ab initio potential energy surfaces for H3O+ and NH3: Approaching the subwave number accuracy for the inversion splittings. J. Chem. Phys. 2003, 118, 10929–10938. [Google Scholar] [CrossRef]
  63. Rajamäki, T.; Miani, A.; Halonen, L. Vibrational energy levels for symmetric and asymmetric isotopomers of ammonia with an exact kinetic energy operator and new potential energy surfaces. J. Chem. Phys. 2003, 118, 6358–6369. [Google Scholar] [CrossRef]
  64. Sawodny, W.; Härtner, H.; Minkwitz, R.; Bernstein, D. Vibrational spectra and force constants for the molecular series NH3−nCln. J. Mol. Struct. 1989, 213, 145–153. [Google Scholar] [CrossRef]
  65. Yurchenko, S.N.; Carvajal, M.; Thiel, W.; Jensen, P. Ab initio dipole moment and theoretical rovibrational intensities in the electronic ground state of PH3. J. Mol. Spectrosc. 2006, 239, 71–87. [Google Scholar] [CrossRef]
  66. Nakamoto, K. Infrared and Raman Spectra of Inorganic Compounds. Part A: Theory and Applications in Inorganic Chemistry, 6th ed.; John Wiley & Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
  67. Sánchez-Castellanos, M.; Álvarez Bajo, O.; Amezcua-Eccius, C.; Lemus, R. Vibrational excitations of arsine in the framework of a local unitary group approach. J. Mol. Spectrosc. 2006, 240, 81–92. [Google Scholar] [CrossRef]
  68. Bürger, H.; Ruland, H.; Demaison, J.; Dréan, P. The ν1 = 1 and ν3 = 1 states of AsF3 studied by high resolution IR and millimeterwave spectroscopy: The equilibrium structure of AsF3. J. Mol. Struct. 2000, 517–518, 105–117. [Google Scholar] [CrossRef]
  69. Fusina, L.; Lonardo, G.D. The Fundamental Bands in the Infrared Spectrum of Stibine (SbH3). J. Mol. Spectrosc. 2002, 216, 493–500. [Google Scholar] [CrossRef]
  70. Molnár, J.; Kolonits, M.; Hargittai, M. Molecular structure of SbF3 and BiF3: An electron diffraction study. J. Mol. Struct. 1997, 413–414, 441–446. [Google Scholar] [CrossRef]
  71. Jerzembeck, W.; Bürger, H.; Breidung, J.; Thiel, W. High resolution infrared spectra of the ν1ν4 bands of BiH3, and ab initio calculations of the spectroscopic parameters. J. Mol. Spectrosc. 2004, 226, 32–44. [Google Scholar] [CrossRef]
  72. Boudon, V.; Rey, M.; Loëte, M. The vibrational levels of methane obtained from analyses of high-resolution spectra. J. Quant. Spectrosc. Radiat. Transf. 2006, 98, 394–404. [Google Scholar] [CrossRef]
  73. Ulenikov, O.; Gromova, O.; Bekhtereva, E.; Raspopova, N.; Fomchenko, A.; Sydow, C.; Bauerecker, S. High resolution study of strongly interacting ν3(F2)/ν1(A1) bands of MSiH4 (M = 28, 29, 30). J. Quant. Spectrosc. Radiat. Transf. 2017, 201, 35–44. [Google Scholar] [CrossRef]
  74. Ulenikov, O.; Gromova, O.; Bekhtereva, E.; Raspopova, N.; Kashirina, N.; Fomchenko, A.; Sydow, C.; Bauerecker, S. High resolution study of MSiH4 (M = 28, 29, 30) in the Dyad Region: Analysis of line positions, intensities and half-widths. J. Quant. Spectrosc. Radiat. Transf. 2017, 203, 496–510. [Google Scholar] [CrossRef]
  75. Sydow, C.; Gromova, O.; Bekhtereva, E.; Raspopova, N.; Belova, A.; Bauerecker, S.; Ulenikov, O. First high–resolution analysis of the fundamental bands of 29SiD4 and 30SiD4: Line positions and strengths. J. Quant. Spectrosc. Radiat. Transf. 2019, 225, 125–155. [Google Scholar] [CrossRef]
  76. Richard, C.; Boudon, V.; Rizopoulos, A.; Auwera, J.V.; Tchana, F.K. Line positions and intensities for the ν2/ν4 bands of 5 isotopologues of germane near 11.5 μm. J. Quant. Spectrosc. Radiat. Transf. 2021, 260, 107474. [Google Scholar] [CrossRef]
  77. Boudon, V.; Grigoryan, T.; Philipot, F.; Richard, C.; Tchana, F.K.; Manceron, L.; Rizopoulos, A.; Auwera, J.V.; Encrenaz, T. Line positions and intensities for the ν3 band of 5 isotopologues of germane for planetary applications. J. Quant. Spectrosc. Radiat. Transf. 2018, 205, 174–183. [Google Scholar] [CrossRef]
  78. Bermúdez-Montaña, M.; Carvajal, M.; Pérez-Bernal, F.; Lemus, R. An algebraic alternative for the accurate simulation of CO2 Raman spectra. J. Raman Spectrosc. 2020, 51, 569–583. [Google Scholar] [CrossRef]
Figure 1. Plot of the parameters δ vs. γ for the series of molecules included in Table 3.
Figure 1. Plot of the parameters δ vs. γ for the series of molecules included in Table 3.
Molecules 29 03490 g001
Figure 2. Plots of ϵ 1 vs. δ and ϵ 1 vs. γ ( ω ) , for the same series of molecules considered in Figure 1.
Figure 2. Plots of ϵ 1 vs. δ and ϵ 1 vs. γ ( ω ) , for the same series of molecules considered in Figure 1.
Molecules 29 03490 g002
Figure 3. Coefficients involved in the relation (22) connecting the normal and local polyads.
Figure 3. Coefficients involved in the relation (22) connecting the normal and local polyads.
Molecules 29 03490 g003
Figure 4. At the left, the plot δ vs. γ ( ω ) is displayed, including the series of molecules BeX 2 and MgX 2 , which manifest a strong normal mode behavior. On the right, the equivalent plot using the parameterization (47), but adding to γ ( ω ) additional contributions up to n-th order indicated with γ n ( ω ) . The curve labeled γ 10 ( ω ) + γ 9 ( i j ) means that both gammas were summed up the indicated order.
Figure 4. At the left, the plot δ vs. γ ( ω ) is displayed, including the series of molecules BeX 2 and MgX 2 , which manifest a strong normal mode behavior. On the right, the equivalent plot using the parameterization (47), but adding to γ ( ω ) additional contributions up to n-th order indicated with γ n ( ω ) . The curve labeled γ 10 ( ω ) + γ 9 ( i j ) means that both gammas were summed up the indicated order.
Molecules 29 03490 g004
Figure 5. Plots of δ vs. γ ( ω ) for the series of isotopologues of H 2 O, SO 2 , and O 3 . A clear linear trend was obtained in all systems. The fundamentals needed to obtain the parameters were obtained from the references indicated: water [55,56,57], sulfur dioxide [58], and ozone [49]. The force constants were obtained from the principal isotopologues.
Figure 5. Plots of δ vs. γ ( ω ) for the series of isotopologues of H 2 O, SO 2 , and O 3 . A clear linear trend was obtained in all systems. The fundamentals needed to obtain the parameters were obtained from the references indicated: water [55,56,57], sulfur dioxide [58], and ozone [49]. The force constants were obtained from the principal isotopologues.
Molecules 29 03490 g005
Figure 6. Slope m vs. x f for the isotopologues associated with the molecules H 2 O, SO 2 , CO 2 , O 3 , and NO 2 . The fundamentals for nitrogen dioxide were obtained from Refs. [59,60].
Figure 6. Slope m vs. x f for the isotopologues associated with the molecules H 2 O, SO 2 , CO 2 , O 3 , and NO 2 . The fundamentals for nitrogen dioxide were obtained from Refs. [59,60].
Molecules 29 03490 g006
Figure 7. Plot of δ vs. γ ( ω ) for the stretching modes of pyramidal molecules.
Figure 7. Plot of δ vs. γ ( ω ) for the stretching modes of pyramidal molecules.
Molecules 29 03490 g007
Figure 8. Plots of δ vs. γ ( ω ) for the bending modes of pyramidal molecules.
Figure 8. Plots of δ vs. γ ( ω ) for the bending modes of pyramidal molecules.
Molecules 29 03490 g008
Figure 9. Plot of δ vs. γ for the stretching degrees of freedom of tetrahedral molecules. On the right, the exponential form is manifested due to the high degree of normality of the molecules.
Figure 9. Plot of δ vs. γ for the stretching degrees of freedom of tetrahedral molecules. On the right, the exponential form is manifested due to the high degree of normality of the molecules.
Molecules 29 03490 g009
Figure 10. Plot of δ vs. γ ( ω ) for tetrahedral molecules involving the bending degrees of freedom.
Figure 10. Plot of δ vs. γ ( ω ) for tetrahedral molecules involving the bending degrees of freedom.
Molecules 29 03490 g010
Figure 11. Plot of δ vs. γ ( 1 ) and δ vs. γ ( 2 ) for tetrahedral molecules involving the bending degrees of freedom.
Figure 11. Plot of δ vs. γ ( 1 ) and δ vs. γ ( 2 ) for tetrahedral molecules involving the bending degrees of freedom.
Molecules 29 03490 g011
Figure 12. Plot of δ vs. γ ( ω ) for the different series of isotopologues.
Figure 12. Plot of δ vs. γ ( ω ) for the different series of isotopologues.
Molecules 29 03490 g012
Figure 13. Plot of m ( ω ) vs. x f associated with the isotopologues of Figure 12. [NH4 ] + , CH4, and SiH4 lie along the line, while GeH4 is manifested outside of it.
Figure 13. Plot of m ( ω ) vs. x f associated with the isotopologues of Figure 12. [NH4 ] + , CH4, and SiH4 lie along the line, while GeH4 is manifested outside of it.
Molecules 29 03490 g013
Figure 14. Plot of δ vs. γ ( ω ) for different octahedral molecules.
Figure 14. Plot of δ vs. γ ( ω ) for different octahedral molecules.
Molecules 29 03490 g014
Figure 15. Plots of δ vs. ϵ 1 and δ vs. ϵ 2 for different octahedral molecules.
Figure 15. Plots of δ vs. ϵ 1 and δ vs. ϵ 2 for different octahedral molecules.
Molecules 29 03490 g015
Figure 16. Plots of Z N ( t ) at different temperatures, showing the dependence of the LN degree. On the right, the partition function Z L ( t ) is shown together with Z N ( t ) . Both coincide in the local limit, as expected.
Figure 16. Plots of Z N ( t ) at different temperatures, showing the dependence of the LN degree. On the right, the partition function Z L ( t ) is shown together with Z N ( t ) . Both coincide in the local limit, as expected.
Molecules 29 03490 g016
Table 1. Calculated force constants f q i q j ( N ) using expressions (3) and f q i q j ( L ) (11) for the different systems we have included in our analyses. Here, the ћ Ω Γ corresponds to the fundamental energy of the Γ -th mode and N stands for the number of oscillators. G Γ Γ ( 0 ) and g q q ( 0 ) are the elements of the Wilson matrix in the normal and local coordinate basis, respectively. The meanings of μ 2 and μ 3 are irreducible representations and will be specified for each case along the text.
Table 1. Calculated force constants f q i q j ( N ) using expressions (3) and f q i q j ( L ) (11) for the different systems we have included in our analyses. Here, the ћ Ω Γ corresponds to the fundamental energy of the Γ -th mode and N stands for the number of oscillators. G Γ Γ ( 0 ) and g q q ( 0 ) are the elements of the Wilson matrix in the normal and local coordinate basis, respectively. The meanings of μ 2 and μ 3 are irreducible representations and will be specified for each case along the text.
Force Constants from (3)Force Constants from (11)
f q q ( N ) = 1 N Γ n Γ Ω Γ 2 G Γ Γ ( 0 ) f q q ( L ) = 1 g q q ( 0 ) 1 N Γ n Γ Ω Γ 2
f q q ( N ) = 1 N Γ ( ) δ Γ , μ 2 Ω Γ 2 G Γ Γ ( 0 ) f q q ( L ) = f q q ( L ) 2 Γ ( ) δ Γ , μ 2 Ω Γ Γ n Γ Ω Γ x g
f q q ( N ) = 1 N Γ ( ) δ Γ , μ 3 n Γ Ω Γ 2 G Γ Γ ( 0 ) f q q ( L ) = f q q ( L ) 2 Γ ( ) δ Γ , μ 3 n Γ Ω Γ Γ n Γ Ω Γ x g
Table 2. Fundamentals and force constants calculated using the expressions in Table 1. References from which the fundamental energies were taken are indicated.
Table 2. Fundamentals and force constants calculated using the expressions in Table 1. References from which the fundamental energies were taken are indicated.
Molecule ν 1 ν 3 f r r ( N ) f r r ( N )
H 2 O [43]3657.0533755.9307.6756−0.0910
H 2 S [44]2614.4082628.4553.9569−0.0168
H 2 Se [45]2344.362357.653.2413−0.0177
SO 2  [46]1151.711362.0610.0004−0.0233
F 2 O [47]928.078314.08640.9397
Cl 2 O [48]641.97686.543.00320.5452
O 3  [49]1103.141042.085.79411.6114
NO 2  [50]1319.791616.8510.28891.8984
CO 2  [51]1285.412349.1414.88710.6891
CS 2  [52]658.001535.357.60210.5767
Table 3. LN parameters for the molecules given in Table 2. Parenthesis ( x ) means × 10 x .
Table 3. LN parameters for the molecules given in Table 2. Parenthesis ( x ) means × 10 x .
Molecule δ γ ( ω ) ζ ϵ 1 ϵ 2
H 2 O5.47 (−7)1.09 (−6)0.01702.19 (−6)−2.73 (−6)
H 2 S6.06 (−7)1.21 (−6)0.00342.42 (−6)−6.42 (−7)
H 2 Se1.73 (−6)3.46 (−6)0.00366.92 (−6)−2.51 (−7)
SO 2 1.68 (−3)3.26 (−3)0.10566.57 (−3)−0.4635
F 2 O8.05 (−3)0.01560.07000.03120.0167
Cl 2 O0.01180.02280.04260.04570.0618
O 3 0.01670.03170.03580.06340.0514
NO 2 0.02110.03860.12710.07790.1567
CO 2 0.03310.04770.33710.10331.2741
CS 2 0.07360.08070.42960.18551.7789
Table 4. Fundamentals (in cm 1 ) and calculated force constants (in aJ Å 2 ) for certain pyramidal molecules. References from which the fundamental energies were taken are indicated.
Table 4. Fundamentals (in cm 1 ) and calculated force constants (in aJ Å 2 ) for certain pyramidal molecules. References from which the fundamental energies were taken are indicated.
Molecule ν 1 ν 2 ν 3 ν 4 f rr ( N ) f rr ( N ) f θ θ ( N ) f θ θ ( N )
14 NH 3  [62]3336.08932.433443.681626.286.3871−0.01170.5252−0.0815
14 NT 3  [63]2014.1656.372184.76996.286.5989−0.01380.5402−0.0645
15 NCl 3  [64]541.7364.8632.3262.82.42750.38370.23210.0571
14 NCl 3  [64]554.2365.26442632.41890.40730.22020.0533
31 PH 3  [65]2321.12992.132326.871118.313.08530.00060.3303−0.0203
PCl 3  [66]515258.35041862.56950.31040.19630.0222
PBr 3  [66]390159.9384.4112.82.05860.35020.09590.0072
75 AsH 3  [67]2115.16906.752126.42999.232.6198−0.00800.2727−0.0139
AsF 3  [68]7403387032634.58240.25790.36440.0489
AsCl 3  [66]416.5192.5319150.22.28600.21800.18290.0211
AsBr 3  [66]272128287991.83260.11050.12360.0079
121 SbH 3  [69]1890.50782.251894.50827.862.0944−0.00250.1929−0.0053
SbF 3  [70]6662506342134.03050.18280.24360.0223
SbCl 3  [66]380.7150.8358.9121.82.17940.15290.13640.0156
BiH 3  [71]1733.25726.71734.47751.241.7632−0.00080.1619−0.0035
BiCl 3  [66]3421233221071.93510.11790.11120.0116
BiBr 3  [66]22077214631.59750.09980.07500.0076
BiI 3  [66]16259.7163.5471.24530.07470.05750.0056
Table 5. Parameters δ , γ ( ω ) , ζ , ϵ i associated with the stretching degrees of freedom for several pyramidal molecules. Parenthesis ( x ) means × 10 x .
Table 5. Parameters δ , γ ( ω ) , ζ , ϵ i associated with the stretching degrees of freedom for several pyramidal molecules. Parenthesis ( x ) means × 10 x .
Molecule δ γ ( ω ) ζ ϵ 1 ϵ 2
14 NT 3 4.00 (−4)5.91 (−4)0.05160.0012−0.3214
SbF 3 4.84 (−4)8.09 (−4)0.03130.00160.0180
AsF 3 8.79 (−4)0.00150.03260.00290.0268
BiBr 3 0.00160.00270.01760.00530.0463
BiI 3 0.00240.00370.00590.00750.0703
AsBr 3 0.00370.00560.03410.01140.1118
PCl 3 0.00720.01170.01370.02310.1138
PBr 3 0.01500.02370.00920.04680.1755
15 NCl 3 0.02640.03480.09750.07260.3282
14 NCl 3 0.02840.03750.09470.07820.3339
Table 6. LN parameters δ , γ ( ω ) , ζ , ϵ i for the bending degrees of freedom of pyramidal molecules. Parenthesis ( x ) means × 10 x .
Table 6. LN parameters δ , γ ( ω ) , ζ , ϵ i for the bending degrees of freedom of pyramidal molecules. Parenthesis ( x ) means × 10 x .
Molecule δ γ ( ω ) ζ ϵ 1 ϵ 2
121 SbH 3 5.09 (−5)7.70 (−5)0.03601.57 (−4)−0.0029
BiH 3 7.12 (−5)1.10 (−4)0.02112.23 (−4)−0.0052
14 NH 3 1.33 (−4)1.09 (−4)0.31642.81 (−4)−0.0012
75 AsH 3 2.60 (−4)3.78 (−4)0.06160.0008−0.0078
31 PH 3 3.64 (−4)5.16 (−4)0.07580.0011−0.0090
PCl 3 1.29 (−6)2.64 (−6)0.20034.85 (−6)1.64 (−5)
AsCl 3 4.40 (−4)8.56 (−4)0.15410.00160.0061
SbCl 3 8.42 (−4)0.00160.13350.00300.0123
BiCl 3 0.00170.00310.08800.00590.0286
14 NCl 3 0.00770.01530.20030.02810.0586
Table 7. Fundamental energies for several tetrahedral molecules, together with the force constants calculated in accordance with Table 1. References from which the fundamental energies were taken are indicated.
Table 7. Fundamental energies for several tetrahedral molecules, together with the force constants calculated in accordance with Table 1. References from which the fundamental energies were taken are indicated.
Molecule ν 1 ν 2 ν 3 ν 4 f rr ( N ) f rr ( N ) f θ θ ( N ) f θ θ ( N ) f θ θ ( N )
12 CH 4  [72]2916.4821533.3363019.4931310.7624.91440.04550.3424−0.0776−0.0321
13 CH 4  [72]-1533.493-1302.781--0.3413−0.0776−0.0310
12 CD 4  [72]-1091.652-997.871--0.3576−0.0786−0.0433
13 CD 4  [72]-1091.801-989.250--0.3565−0.0786−0.0422
[ 14 NH 4 ] +  [66]30401680314514005.39140.03220.4024−0.0931−0.0300
[ 15 NH 4 ] +  [66]-1646-1399--0.3957−0.0894−0.0381
[ND 4 ] +  [66]-1215-1065--0.4281−0.0973−0.0387
[NT 4 ] +  [66]-976-913--0.4290−0.0940−0.0528
28 SiH 4  [73,74]2186.8723970.934452189.1895913.468792.74660.03110.1576−0.0311−0.0332
29 SiH 4  [74]-970.94842-912.18312--0.1575−0.0311−0.0331
30 SiH 4  [74]-970.96148-910.97961--0.1574−0.0311−0.0329
29 SiD 4  [75]-689.88679-672.93384--0.1628−0.0314−0.0373
30 SiD 4  [75]-689.89950-671.43227--0.1626−0.0314−0.0371
70 GeH 4  [76]-929.90124-821.54462--0.1360−0.0285−0.0219
72 GeH 4  [76]-929.90513-821.11703--0.1360−0.0285−0.0218
73 GeH 4  [76]-929.90728-820.91126--0.1359−0.0285−0.0218
74 GeH 4  [76,77]2110.70051929.909102111.14192820.711652.61090.01150.1359−0.0285−0.0218
76 GeH 4  [76]-929.91308-820.32666--0.1359−0.0285−0.0217
CF 4  [66]908.4-1283.0-6.75450.8276---
CCl 4  [66]460-792.765-3.10000.4401---
CBr 4  [66]267-672-2.45460.3006---
CI 4  [66]178-555-1.73720.2106---
HfCl 4  [66]-101.5-112--0.0672−0.0120−0.0193
[PH 4 ] +  [66]-1086-974--0.1867−0.0389−0.0311
[AsH 4 ] +  [66]-1024-941--0.1729−0.0346−0.0345
Table 8. LN parameters associated with the stretching degrees of freedom of tetrahedral molecules.
Table 8. LN parameters associated with the stretching degrees of freedom of tetrahedral molecules.
Molecule δ γ ζ ϵ 1 ϵ 2
74 GeH 4 2.48 (−5)7.92 (−5)0.00015.94 (−5)0.0090
28 SiH 4 1.64 (−4)5.25 (−4)0.00073.94 (−4)0.0236
[ 14 NH 4 ] + 2.63 (−4)8.03 (−4)0.02166.13 (−4)0.0708
12 CH 4 4.03 (−4)1.23 (−3)0.02219.40 (−4)0.0702
CF 4 0.05450.10680.20970.09200.6536
CCl 4 0.10380.15280.31090.14140.9123
CBr 4 0.19760.16990.45310.18351.4336
CI 4 0.28280.18130.50900.20991.6816
Table 9. LN parameters associated with the bending degrees of freedom of tetrahedral molecules.
Table 9. LN parameters associated with the bending degrees of freedom of tetrahedral molecules.
Molecule δ γ ( ω ) γ ( 1 ) γ ( 2 ) ζ ϵ 1 ϵ 2 ϵ 3
13 CH 4 4.96 (−7)6.05 (−7)1.82 (−7)1.21 (−7)0.1027−0.06559.40 (−4)−0.7303
12 CH 4 3.00 (−6)3.64 (−6)1.09 (−6)7.28 (−7)0.0988−0.06116.87 (−4)−0.6576
[ 14 NH 4 ] + 1.10 (−5)1.34 (−5)4.03 (−6)2.69 (−6)0.1145−0.06887.02 (−4)−0.9332
[ 15 NH 4 ] + 1.29 (−5)1.57 (−5)4.71 (−6)3.14 (−6)0.1024−0.06671.43 (−3)−0.7056
13 CD 4 1.87 (−4)2.19 (−4)6.56 (−5)4.38 (−5)0.0625−0.01685.00 (−4)−0.1460
12 CD 4 2.56 (−4)2.98 (−4)8.94 (−5)5.96 (−5)0.0570−0.00949.76 (−4)−0.0849
[ND 4 ] + 3.25 (−4)3.87 (−4)1.16 (−4)7.74 (−5)0.0833−0.02238.30 (−4)−0.2548
HfCl 4 4.65 (−4)4.81 (−4)1.44 (−4)9.63 (−5)0.06247.69 (−3)2.46 (−3)0.0206
[NT 4 ] + 6.85 (−4)7.80 (−4)2.34 (−4)1.56 (−4)0.04240.01644.42 (−3)0.1018
[PH 4 ] + 1.01 (−3)1.22 (−3)3.65 (−4)2.44 (−4)0.0690−0.09860.0153−0.6693
70 GeH 4 1.30 (−3)1.58 (−3)4.73 (−4)3.16 (−4)0.0784−0.11270.0198−0.8034
72 GeH 4 1.30 (−3)1.58 (−3)4.75 (−4)3.16 (−4)0.0787−0.11300.0199−0.8073
73 GeH 4 1.30 (−3)1.58 (−3)4.75 (−4)3.17 (−4)0.0789−0.11320.0199−0.8093
74 GeH 4 1.31 (−3)1.59 (−3)4.76 (−4)3.17 (−4)0.0790−0.11330.0199−0.8111
76 GeH 4 1.31 (−3)1.59 (−3)4.77 (−4)3.18 (−4)0.0793−0.11360.0200−0.8147
29 GeD 4 1.41 (−3)1.63 (−3)4.88 (−4)3.25 (−4)0.0158−0.07010.0162−0.3603
30 GeD 4 1.44 (−3)1.66 (−3)4.98 (−4)3.32 (−4)0.0173−0.07170.0166−0.3705
28 GeH 4 1.80 (−3)2.13 (−3)6.38 (−4)4.25 (−4)0.0388−0.09450.0226−0.5329
29 GeH 4 1.81 (−3)2.15 (−3)6.44 (−4)4.29 (−4)0.0397−0.09540.0229−0.5403
30 GeH 4 1.83 (−3)2.16 (−3)6.49 (−4)4.33 (−4)0.0405−0.09630.0231−0.5474
[AsH 4 ] + 2.16 (−3)2.59 (−3)7.76 (−4)5.17 (−4)0.0537−0.11210.0283−0.6754
Table 10. Fundamentals (in cm 1 ) from Ref. [66] and force constants (in aJ Å 2 ) for octahedral molecules.
Table 10. Fundamentals (in cm 1 ) from Ref. [66] and force constants (in aJ Å 2 ) for octahedral molecules.
Molecule ν 1 ν 2 ν 3 f rr ( N ) f rr ( N ) f rr ( N )
PuF66255196123.54760.2262−0.0801
NpF66465256183.64910.2643−0.0352
[PtF6 ] 2 6115765713.46170.07750.4071
[PbCl6 ] 2 2812092621.11320.12280.0449
[HfCl6 ] 2 3252572751.39280.13780.2623
[GaF6 ] 3 5353984811.96310.23850.2869
[SbBr6 ] 1921692391.31890.06510.1560
[SiF6 ] 2 6634777412.97510.39560.3629
[AlF6 ] 3 5414005681.89280.24750.3932
[TiCl6 ] 2 3202713161.28810.10080.4475
Table 11. LN parameters associated with the stretching degrees of freedom of octahedral molecules.
Table 11. LN parameters associated with the stretching degrees of freedom of octahedral molecules.
Molecule δ γ ( ω ) γ ( 1 ) γ ( 2 ) ζ ϵ 1 ϵ 2 ϵ 3
PuF60.00190.00230.00360.00410.11670.00450.0541−0.1762
NpF60.00240.00310.00500.00520.12970.00590.0645−0.5237
[PtF6 ] 2 0.00310.00560.00260.00050.04300.01100.11050.0117
[PbCl6 ] 2 0.00710.01040.01640.01220.18200.01850.13740.2959
[HfCl6 ] 2 0.01200.02060.02240.00980.14610.03770.20490.0673
[GaF6 ] 3 0.01420.02360.02930.01480.18180.04220.22120.1258
[SbBr6 ] 0.02300.03430.01400.00240.21040.06910.30480.2421
[SiF6 ] 2 0.02950.04340.04380.01770.26040.08070.33030.3378
[AlF6 ] 3 0.03680.05670.04920.01710.21270.10620.36700.2481
[TiCl6 ] 2 0.04810.07770.03330.00610.10460.15220.41350.1936
Table 12. Maximum components of the wave function compared with δ for the pyramidal molecules PH 3 , AsH 3 , and SbH 3 , taken from Ref. [20].
Table 12. Maximum components of the wave function compared with δ for the pyramidal molecules PH 3 , AsH 3 , and SbH 3 , taken from Ref. [20].
δ × 10 6 0.1430.6851.042
State SbH3PH3AsH3
| 200000 * 0.9890.9680.951
| 400000 * 0.9810.6270.555
| 200100 0.4670.7560.722
Table 13. Transition involved in the most intense lines in the Raman spectrum. Labels Σ + ( Π ± ) and Σ g + ( Π u ± ) for the asymmetric and symmetric isotopologues, respectively. Given the normalization, the parameters a and a determine the eigenstates and the transition strength.
Table 13. Transition involved in the most intense lines in the Raman spectrum. Labels Σ + ( Π ± ) and Σ g + ( Π u ± ) for the asymmetric and symmetric isotopologues, respectively. Given the normalization, the parameters a and a determine the eigenstates and the transition strength.
Symmetry | ν | ν
Σ + , Σ g + | 00 0 0 a | 10 0 0 + b | 02 0 0
Σ + , Σ g + | 00 0 0 b | 10 0 0 a | 02 0 0
Π ± , Π u ± | 01 ± 1 0 a | 11 ± 1 0 + b | 03 ± 1 0
Π ± , Π u ± | 01 ± 1 0 b | 11 ± 1 0 a | 03 ± 1 0
Table 14. Components of the wave functions, as well as the parameter δ . The notation for the isotopologues is 1 x O 1 y C 1 z O ( x y z ) .
Table 14. Components of the wave functions, as well as the parameter δ . The notation for the isotopologues is 1 x O 1 y C 1 z O ( x y z ) .
(636)(637)(638)(626)(738)(627)(628)(727)(728)(828)
a0.55870.6048590.6453120.69830.7477150.750710.7943560.78260.8302450.8488
a 0.60760.6362410.6562850.70450.7220740.7437650.7784240.75860.8057480.8150
δ 0.033750.034870.035940.36840.037170.0380350.039150.039330.040510.04181
Table 15. Square of the maximum components for the wave functions involved in the multiplets (84). The components were taken from Ref. [54].
Table 15. Square of the maximum components for the wave functions involved in the multiplets (84). The components were taken from Ref. [54].
δ 0.026260.028940.03035
PolyadStates (934)(924) (925)
2 | 02 0 0 | 10 0 0 0.9610.9440.937
4 | 04 0 0 | 12 0 0 0.8750.8340.816
4 | 20 0 0 | 00 0 1 0.8440.8270.811
3 | 03 1 0 | 11 1 0 0.9300.9040.890
5 | 05 1 0 | 13 1 0 0.8200.7810.750
5 | 21 1 0 | 01 1 1 0.7560.7350.712
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Suárez, E.; Guzmán-Juárez, O.; Lemus, R. Novel Criteria to Provide a Locality/Normality Degree in Molecules and Their Relevance in Physical Chemistry. Molecules 2024, 29, 3490. https://doi.org/10.3390/molecules29153490

AMA Style

Suárez E, Guzmán-Juárez O, Lemus R. Novel Criteria to Provide a Locality/Normality Degree in Molecules and Their Relevance in Physical Chemistry. Molecules. 2024; 29(15):3490. https://doi.org/10.3390/molecules29153490

Chicago/Turabian Style

Suárez, Eduardo, Oscar Guzmán-Juárez, and Renato Lemus. 2024. "Novel Criteria to Provide a Locality/Normality Degree in Molecules and Their Relevance in Physical Chemistry" Molecules 29, no. 15: 3490. https://doi.org/10.3390/molecules29153490

APA Style

Suárez, E., Guzmán-Juárez, O., & Lemus, R. (2024). Novel Criteria to Provide a Locality/Normality Degree in Molecules and Their Relevance in Physical Chemistry. Molecules, 29(15), 3490. https://doi.org/10.3390/molecules29153490

Article Metrics

Back to TopTop