5.2
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
Corrigendum
Current Issue
Editorial
Erratum
Full Length Article
Full lenth article
Letter to Editor
Original Article
Research article
Retraction
Retraction notice
Review
Review Article
SPECIAL ISSUE: ENVIRONMENTAL CHEMISTRY
5.3
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
Corrigendum
Current Issue
Editorial
Erratum
Full Length Article
Full lenth article
Letter to Editor
Original Article
Research article
Retraction
Retraction notice
Review
Review Article
SPECIAL ISSUE: ENVIRONMENTAL CHEMISTRY
View/Download PDF

Translate this page into:

Original article
12 (
8
); 4612-4626
doi:
10.1016/j.arabjc.2016.06.016

Shedding light on molecular structure, spectroscopic, nonlinear optical and dielectric properties of bis(thiourea) silver(I) nitrate single crystal: A dual approach

Advanced Functional Materials & Optoelectronic Laboratory (AFMOL), Department of Physics, College of Science, King Khalid University, P.O. Box 9004, Abha 61413, Saudi Arabia
Research Center for Advanced Materials Science (RCAMS), King Khalid University, P.O. Box 9004, Abha 61413, Saudi Arabia
Department of Chemistry, Faculty of Science, King Khalid University, P.O. Box 9004, Abha 61413, Saudi Arabia

⁎Corresponding authors at: Department of Physics, College of Science, King Khalid University, Abha 61413, Saudi Arabia. Fax: +966 72418319. shkirphysics@gmail.com (Mohd. Shkir), shkirphysics@kku.edu.sa (Mohd. Shkir), shabbir193rb@gmail.com (Shabbir Muhammad)

Disclaimer:
This article was originally published by Elsevier and was migrated to Scientific Scholar after the change of Publisher.

Peer review under responsibility of King Saud University.

Abstract

The current work is to spotlight the key structure-property features of a novel thiourea complex bis(thiourea) silver(I) nitrate (BTSN) for future applications in nonlinear optical devices through dual approach involving experimental and theoretical techniques. The synthesis and growth of good quality single crystals of a relatively new material BTSN have been done. The crystal structure, quantitative and qualitative, vibrational, optical and dielectric analysis of the grown single crystals was carried out. The optical transparency of the grown crystals is found to be more than 80% confirms its colorless nature and applications in optoelectronic devices. Further, the state-of-art computational methods have been applied to obtain the ground state molecular geometry at B3LYP/6-31G∗, MP2/6-31G∗ and M06/6-31G∗ levels of theory. Various key electro-optical properties (complementary to experimental results) such as IR, Raman, and polarizabilities have been determined at the same levels of theory. The static and dynamic polarizability and first hyperpolarizability both were calculated to comprehend the potential applications of BTSN in nonlinear optics. In addition to the above, numerous novel molecular level insights have been achieved in the form of total and partial density of states, HOMO-LUMO gap and molecular electrostatic potential map. The calculated values of static and frequency dependent dynamic first hyperpolarizability are found to be 8.46 × 10−30 and 2.30 × 10−30 esu which are about 23 times and 13 times greater than those of prototype urea molecule, respectively, calculated at the same B3LYP/6-31G∗ level of theory. From the analyzed results it is clear that the titled compound possesses excellent electro-optic properties that make it a decent contestant for photonic devices.

Keywords

Crystal growth
Optical properties
Nonlinear optical material
Density functional theory
1

1 Introduction

Over the past few decades, various kinds of materials have been functionalized to match the demand of functional materials with unique optical and nonlinear optical (NLO) properties (Marder et al., 1991; Papadopoulos et al., 2006; Muhammad et al., 2013). Since the invention of photon as the fastest carrier of information, the field of nonlinear optics has received an enormous consideration by the researchers and photonic industry. To satisfy day to day technological requirements the scientists are continuously growing the single crystals of new as well as existing nonlinear optical (NLO) materials with exceptional optoelectronic properties which can be used in various applications such as microelectronics, semiconductors, photonics, sensors, laser technology, and optical communication (Saleh, 1991; Abbas et al. 2019; Penn et al., 1991; Shkir and Abbas, 2014; Shakir et al., 2009c; Zyss et al., 1984; Fujiwara et al., 2006; Muhammad et al., 2015; Shkir et al., 2015b). As per the current literature, the research and development is attentive on new kind of NLO materials known as semiorganic class of materials that combine the properties of organic and inorganic materials (Muhammad et al., 2009, 2013; Shkir et al., 2015a). Thiourea is a well-known material in the history belongs to centrosymmetric material category, and has large dipole moment and due to this it has the ability to modify the inorganic as well organic matrix and forms various applicable metal complexes with noncentrosymmetric structure (Allen, 2002; Petrova et al., 1997; Nardelli et al., 1965; Crow, 1986; Burrows et al., 1999; Campo et al., 2004; Jiang and Fang, 1999; Zachariadis et al., 2003).

Recently, in addition to the above reported thiourea complexes, a new thiourea complex Bis(thiourea) silver(I) nitrate (BTSN) has been reported by Sivakumar et al. (2014) in which they have described the synthesis, crystal growth, structure, spectroscopic, nonlinear optical and mechanical properties. The above literature shows that there has been a focus on preliminary experimental studies of BTSN molecule.

However, the first time use of a twofold approach including experimental and computational techniques will emphasize the numerous key features of BTSN such as its molecular geometry, pattern of frontier molecular orbitals (FMOs), linear as well as static and dynamic nonlinear optical properties, which are unknown to scientific community so far. For instance, the main aim of the present investigation was to deliver the key insight at molecular level to sightsee the prospective of BTSN as an advanced functional material. Moreover, the synthesis, crystal growth, surface morphology/elemental composition, optical and dielectric analysis will be performed.

Therefore, the main goal of the present investigation is many folds such as the proper synthesis and its single crystal growth by slow evaporation technique and subjected to powder X-ray diffraction and Fourier transform – Raman spectroscopy study for structural and vibrational analysis. The quality check or crystalline perfection of the grown single crystals was performed by UV–vis.-NIR, diffuse reflectance spectroscopy and dielectric measurements. The dielectric studies were studied in the high frequency range at ambient temperature. In addition to the above reported characterizations, the stable molecular geometry optimization is performed using diverse innovative computational methods. The Infrared, Raman, ultraviolet–visible spectra, polarizability, static and dynamic first hyperpolarizability are calculated using density functional theory (DFT), time dependent DFT (TD-DFT), and finite field (FF) methodologies. The computational and experimentally attained results are discussed and likened with each other as well as compared with the reported data wherever possible.

2

2 Experimental and theoretical details

2.1

2.1 Synthesis

The synthesis of bis(thiourea) silver(I) nitrate (BTSN) was successfully done by chemical reaction between commercially available thiourea and silver nitrate purchased from Sigma Aldrich by taking their equimolar ratio (2:1). Primarily the calculated amount of thiourea (CH4N2S) and Silver nitrate (AgNO3) was dissolved separately in double distilled water and filtered, then was mixed slowly in another beaker and again stirred well for a long time using a temperature controlled magnetic stirrer above the room temperature for proper chemical reaction and to yield a homogeneous mixture of solution. The reaction mechanism involved in the synthesis of BTSN is as follows: 2 ( CH 4 N 2 S ) Thiourea + AgNO 3 Silver Nitrate Ag [ CH 4 N 2 S ] 2 . NO 3 Bis ( thiourea ) silver nitrate

Then the prepared solution was evaporated at room temperature to yield the crystalline powder salt of Bis(thiourea) silver nitrate (BTSN).

2.2

2.2 Crystal growth

By using the synthesized materials the good quality single crystal of BTSN was grown by slow evaporation solution technique from the aqueous solutions. Saturated solution of BTSN was prepared at high temperature ∼35 °C. The continuous stirring was done for more than 72 h to get transparent and homogeneous solution of the titled compound. The prepared solution was then filtered twice in another flat surface beaker using good quality filter paper and covered with a perforated lid and was housed in a constant temperature bath with temperature stability ±0.01 °C at constant temperature (35 °C). After one day we have started to reduce the temperature up to room temperature and left for crystal growth. Good quality and optically transparent single crystals were harvested after a span of 15 days from the beaker as shown in Fig. 1.

As grown single crystals of BTSN.
Figure 1 As grown single crystals of BTSN.

2.3

2.3 Characterization

The homogeneous powder of the grown crystals of BTSN was prepared for performing the powder X-ray diffraction analysis on Shimadzu X-600 Japan powder X-ray diffractometer (PXRD) having Cu Kα radiation (40 kV, 30 mA, λ  = 0.1543 nm). The specimen was scanned at the scan rate of 0.002°/s in the angular range of 2θ (5–90°) to find the crystal structure and lattice parameters of BTSN.

The surface topography of the samples was investigated by using scanning electron microscope (JSM 6360 LA, Japan) equipped with energy dispersive analytical X-ray unit (EDAX). EDAX was used to find out the elemental composition of the samples. Before SEM examination, the samples were sputter coated with 10 nm of platinum to eliminate charge buildup.

FT-Raman spectrum was recorded for BTSN single crystal to study the vibrational modes. To record the spectrum the laser beam was made to incident normally on the most prominent surfaces of the grown crystal and scattered intensity was collected at 300 K using THERMO SCIENTIFIC, DXR FT-RAMAN spectrometer coupled with microscope using 532 nm laser of power up to 10 mW, full range grating (3500–5 cm−1), at estimated resolution 5.1–8.3 cm−1 and the size of the aperture pinhole was kept 50 μm.

For optical testing BTSN single crystal of thickness ∼1.2 mm and its powder in a holder of thickness 1.5 cm were subjected to JASCO V-570 UV–vis.-NIR and Shimadzu UV–VIS-NIR spectrophotometer model (UV-3600) spectrophotometer to measure its optical absorbance and diffuse reflectance (DR) in wavelength range of 190–2300 nm at 300 K, respectively.

The dielectric studies were carried out in the higher frequency range i.e. 1 kHz–10 MHz at constant temperature under the dark and UV light using a KEITHLEY 4200-SCS and indigenously designed and developed UV light system.

2.4

2.4 Computational details

All the calculations have been performed using Gaussian suit of programs (Frisch et al., 2009). The equilibrium geometry of BTSN was successfully achieved corresponding to the true minimum on the potential energy surface (PES) by solving self-consistent field equation. The density functional theory (DFT) method has been used to optimize the molecular geometry of BTSN at B3LYP and hybrid Gen basis set. The hybrid Gen basis consists of explicitly two types of basis sets i.e. 6-31G∗ and LANL2DZ. For instance, 6-31G∗ has been implemented on carbon, nitrogen, oxygen and hydrogen atoms while LANL2DZ basis set has been applied for silver metal atom. The LANL2DZ basis set LANL (Los Alamos National Laboratory) provides effective core potential (ECP) for heavy atoms (Ag) and the double zeta (DZ) type basis set is utilized as valence basis set (Wadt and Hay, 1985). The use of ECP for the core electrons simplifies the problem by eliminating the need for the core basis functions, which requires large computation, and at the same time the accuracy is also maintained. This basis set takes into account possible relativistic effects due to the presence of the metal atoms. In addition to B3LYP, highly correlated method MP2 has also used to optimize the initial molecular geometry for comparison of results with B3LYP method. The optimized structural parameters of BTSN were used to characterize all stationary points as minima using IR and Raman frequency calculations. The IR and Raman spectra of BTSN were calculated by taking the second derivative of the energy that is computed analytically. Similarly, all other electronic properties including dipole moment, polarizability, first hyperpolarizability of BTSN were calculated at the B3LYP/6-31G∗ level of theory. The time dependent density functional theory (TD-DFT) has been used to calculate excitation energies with relatively larger basis set of 6-31+G∗ in place of 6-31G∗ basis set. The static first hyperpolarizability (βtot) and its components for BTSN were calculated by finite field (FF) approach at B3LYP/6-31G∗ level of theory. The FF method has been broadly applied to investigate the NLO properties of organic materials because this methodology can be used in concert with the electronic structure method to calculate β values. Recently, several βtot amplitudes calculated by this method are found to be in a semi quantitative agreement with experimental structure property relationship (Muhammad et al., 2010, 2015). In FF method, a molecule is subjected to a static electric field (F), and the energy (E) of the molecule is expressed by the following equation:

(1)
E = E ( 0 ) - μ 1 F 1 - 1 2 α ij F i F j - 1 6 β ijk F i F j F k - 1 24 γ ijkl F i F j F k F l - where E(0) is the energy of molecule in the absence of an electric field, μ is the components of the dipole moment vector, α is the linear polarizability, β and γ are the first and second hyperpolarizabilities, respectively, while i, j and k label the x, y and z components respectively. It is clear from the above equation that the values of μ, α, β, and γ can be obtained by differentiating E with respect to F.

3

3 Results and discussion

3.1

3.1 Structural analysis

3.1.1

3.1.1 Powder X-ray diffraction analysis

Powder X-ray diffraction technique is one of the best techniques to analyze the crystalline nature as well as crystal structure of any material. For better understanding and comparison between thiourea and BTSN we have subjected both for powder X-ray diffraction analysis (PXRD). The recorded PXRD pattern of pure thiourea as well as BTSN is shown in Fig. 2(a). From the PXRD patterns (Fig. 2(a)), it is clearly visible that the pattern of BTSN is totally different than that of thiourea, which confirms its formation. Moreover, the XRD pattern of current work with previously reported (Sivakumar et al., 2014) is also shown in Fig. 2(b). It is clear from Fig. 2(b) that the PXRD pattern is similar to the reported one (CCDC# 913226) which confirms the similar structure of the currently grown crystals. The obtained X-ray data were used to confirm the crystal structure and to calculate the lattice parameters using ‘POWDERX’ and ‘CHECKCELL’ software’s. The crystal system, lattice parameters, space group, etc., are presented in Table 1 and were found in close covenant with the earlier report (Sivakumar et al., 2014). From figure it is clear that the peak intensities are very high which suggest good crystallinity of the grown crystals and its usefulness.

(a) Powder X-ray diffraction patterns of BTSN with thiourea and (b) BTSN with earlier reports (CCDC#913226).
Figure 2 (a) Powder X-ray diffraction patterns of BTSN with thiourea and (b) BTSN with earlier reports (CCDC#913226).
Table 1 Lattice parameters of BTSN crystal.
Lattice parameters Reported work (Sivakumar et al., 2014) Current work (powder XRD)
Single crystal XRD Powder XRD POWDERX CHECK CELL
a (Å) 33.3455(6) 33.3440 33.3547 33.3459
b (Å) 45.2957(7) 45.2916 45.2536 45.2858
c (Å) 20.3209(5) 20.3064 20.3565 20.3249
V (Å)3 30692.8(10) 30,666 30729.038 30692.545
α (°) 90 90 90 90
β (°) 90 90 90 90
γ (°) 90 90 90 90
System Orthorhombic Orthorhombic Orthorhombic Orthorhombic
S.G. C2221 C2221 C2221 C2221

3.2

3.2 Scanning Electron Microscope (SEM) and energy dispersive X-ray spectroscopy (EDXS) analysis

SEM/EDXS is found to be a very effective technique to analyze the qualitative (surface morphology/inclusions) as well as quantitative (the presence of elements in any crystalline matrix) of single crystals (Shkir et al., 2014a, 2014b).

Different crystals of BTSN and its various places were chosen to record the SEM micrographs and EDXS patterns. The captured SEM micrographs and recorded EDXS patterns of BTSN crystal along prominent plane are shown in Fig. 3. The recorded EDXS pattern reveals the elemental or atomic percentage of C (9.50%), N (43.06%), O (26.98%), S (13.64%), Ag (6.24%) and Pt (0.57% due to Pt coating). The energy (keV), mass %, error % and atomic % of all the elements present in BTSN are presented in Table 2 with theoretical values. As clear from the EDX and theoretically calculated data available in Table 2 that there is a difference in both the values of each element composition. The difference in experimental and theoretical elemental composition is may be due to inability of EDX to detect lighter atoms such as H in the system (Srinivasan and Narvekar, 2016). However, the experimental data obtained from EDX study are very useful and supportive as it is a good technique to identify the elemental composition in any complex (Kalaivani et al., 2016). SEM micrographs were captured at different resolution 27×, 30× and 650× and show that the surface of grown crystal is of good quality.

EDX spectra and SEM images of BTSN single crystals.
Figure 3 EDX spectra and SEM images of BTSN single crystals.
Table 2 The percentage composition of elements in BTSN crystal.
Elements E (keV) Mass (%) Error (%) At (%) Theoretical mass (%)
C 0.277 4.81 0.39 9.50 C-7.99
N 0.392 25.43 1.48 43.06 H-2.68
O 0.525 18.20 1.16 26.98 N-22.70
S 2.307 18.44 0.19 13.64 O-13.96
Ag 2.983 28.39 0.70 6.24 S-21.27
Pt 2.048 4.72 0.80 0.57 Ag-31.39
Total 100.00 100.00 99.99

3.3

3.3 Molecular geometry

In the present investigation, we have adopted a tetrahedral coordination of thiourea with silver atom as in it has been suggested in some previous studies including single crystal analysis of BTSN (Sivakumar et al., 2014). The initial geometry of monomeric molecule BTSN has been adopted from its experimental crystallographic coordinates (see Fig. S1 of supporting information) and fully optimized to get the lowest energy global minimum structure. The optimized geometrical structure of monomeric BTSN has been illustrated in Fig. 4 as calculated at B3LYP/6-31G∗ level of theory. The geometrical parameters including bond lengths and bond angles are calculated at two different levels of theory as collected in Table 3. Table 3 shows that the B3LYP and MP2 methods in conjugation with hybrid Gen basis sets (as explained in computational details section) are used to study the geometrical parameters of BTSN. A comparison of results from Table 3 shows that the Ag–S bond lengths as at B3LYP method are in reasonable agreement with highly correlated MP2 method and also available experimental parameters. The slight variations among all four Ag–S bond lengths from average Ag–S bond length (2.751 Å) are perhaps due to the steric orientations of thiourea ligands. The nitrate ion has been stabilized through hydrogen bonding with two of the thiourea ligands. The hydrogen bond distances between thiourea and nitrate ion might be considered in medium (1.5–2.2 Å) range of hydrogen bonds (Scheiner, 1997). It is important to mention here that due to the gas phase optimization of BTSN monomeric structure and its structure relaxation from polymer to monomer, its nitrate ion configuration (hydrogen bond angles) has been changed or became irrelevant (for comparison) from its original experimental positions as shown in Fig. S1 of supporting information, while the other details of different bond angles have been remain relevant and shown in Table 3 for comparison purpose.

Optimized molecular geometry of monomeric BTSN at B3LYP/6-31G∗ level of theory.
Figure 4 Optimized molecular geometry of monomeric BTSN at B3LYP/6-31G∗ level of theory.
Table 3 Important bond lengths [Å] and bond angles [°] of BTSN optimized molecule.
Bond lengths (Å) Methods Bond angle (°) Methods
B3LYP MP2 Exp. B3LYP MP2 Exp.
Ag1–S18 2.816 2.743 2.700 S2–Ag1–S18 102.29 102.47 105.16
Ag1–S10 2.671 2.669 2.700 S2–Ag1–S26 98.82 101.54 93.87
Ag1–S2 2.631 2.631 2.484 S10–Ag1–S26 97.48 98.16 93.73
Ag1–S26 2.886 2.806 2.700 S18–Ag1–S26 107.70 102.52 108.39
S2–C3 1.753 1.736 1.738 C3–S2–Ag1 93.55 90.53 108.32
S10–C11 1.719 1.703 1.719 C19–S18–Ag1 95.72 92.65 114.13
S18–C19 1.722 1.708 1.713 C11–S10–Ag1 105.34 103.13 105.90
S26–C27 1.709 1.694 1.709 C27–S26–Ag1 105.47 105.45 106.454
H6⋯O35 1.806 1.836 2.110 N21–H25–O37 158.86 159.09
H25⋯O37 1.747 1.758 2.216 N4–H6–O35 173.20 171.76

3.4

3.4 Experimental and theoretical IR and FT-Raman spectroscopic analysis

To characterize the organic, inorganic and organic-inorganic hybrid materials infrared (IR) and Raman spectroscopy plays a vital role. These techniques help us to identify the presence of functional groups as well as to understand different molecular conformations and reaction mechanism. Raman Spectra of BTSN have been recorded in the wavenumber range of 3500–50 cm−1 as shown in Fig. 5(a). There is no theoretical report available so far on the titled compound; therefore, the IR and Raman spectra were calculated theoretically for the first time by DFT method as explained in methodology section and shown in Figs. 5(b) and 6. To compare the theoretically calculated IR and Raman spectra with experimental we have used a scaling factor 0.9777 (Andrade et al., 2008). It is considerable that several chemical groups can be formed in titled compound due to the coordination of metal with thiourea and that may be possible either through bonding to the nitrogen or sulfur of the thiourea (Selvakumar et al., 2007; Kannan et al., 2004). Some previous reports have already explained the bond formation patterns in the titled compound (Sivakumar et al., 2014). The vibrations of thiourea and silver ion coordinated tetrahedral bonds show its Raman peaks in the lower frequency region i.e. 900–50 cm−1. It is clear from the molecular geometry of BTSN molecule (Fig. 4) that a tetrahedral arrangement forms as the silver atom coordinates to sulfur atom from thiourea and to oxygen atoms from nitrate group. For pure thiourea generally NH2 symmetric stretching vibration occurs in the range of 3400–3100 cm−1 (Selvakumar et al., 2007; Angeli Mary and Dhanuskodi, 2001). Experimentally, the symmetric and asymmetric stretching vibrations of NH2 occur at 3220 and 3343 cm−1, while theoretically these bands are observed at 3062, 3106, 3176, 3194 and 3238 respectively. The asymmetric and symmetric stretching vibrations of C–N in BTSN are observed at 1510 and 1055 cm−1, while theoretically these bands are observed at 1504 and 1082 cm−1 respectively.

(a) Experimental and (b) theoretically calculated Raman spectra of BTSN.
Figure 5 (a) Experimental and (b) theoretically calculated Raman spectra of BTSN.
Calculated IR spectrum of BTSN at B3LYP/6-31G∗ level of theory.
Figure 6 Calculated IR spectrum of BTSN at B3LYP/6-31G∗ level of theory.

The asymmetric and symmetric vibrations of C–S stretching vibrations are observed at 1410 and 722 cm−1 experimentally, while theoretically at 1381 and 712 cm−1 due to the formation of BTSN complex as thiourea is bonded to silver (Ag+). The N–C–N vibration band is observed at 477 cm−1 in experimental Raman spectrum, while theoretically at 466 cm−1 in metal complex BTSN crystals. The stretching vibration band of AgS occurs at 177 cm−1 (Bowmaker et al., 2009; Sivakumar et al., 2014) while theoretically at 193 cm−1. Spectral features associated with N–C–S, and S–M vibrations are observed due to the delocalization of lone pairs of S atoms (Jayalakshmi et al., 2006). The experimental Raman vibration bands are in agreement with earlier report (Sivakumar et al., 2014). In Fig. 5(a) and (b) the band occurs at 1639 cm−1 in experimental Raman spectrum is due to NH2 bending vibrations while this vibration occurs at 1627 cm−1 in theoretically calculated Raman spectrum. The band at 1636 cm−1 in IR spectrum (Fig. 6) was also observed due to above said vibrations. The low intensity vibration band at ∼1480 cm−1 in experimental Raman and at 1487 cm−1 in theoretically calculated IR spectrum occurs due to C–N asymmetric stretching vibration. The two bands appeared in the lower energy region at ∼1082 and 712 cm−1 in theoretical spectrum may be assigned to C–N symmetric and C–S stretching vibrations, respectively (Kumari et al., 2009; Sivakumar et al., 2014). The asymmetric bending vibrations of N–C–S, N–C–N are observed at 624 and 475 cm−1 respectively which are found to be shifted in comparison of pure thiourea (Kumari et al., 2009), indicating the formation of thiourea metal complex.

3.5

3.5 Optical studies

The study of optical light absorbance of the grown crystal is key parameter to understand the suitability of the grown crystal of titled material for optical applications. The UV–Vis–NIR spectra of BTSN crystal were recorded as shown in Fig. 6 as we can get crucial information about its structure because the absorption of ultra violet and visible light involves the promotion of electron in σ and π orbital from the ground to higher energy states. A nonlinear optical material can be of great practical use only if it has wide transparency window; therefore, it is essential to study the ultraviolet spectra of any NLO crystal. The grown single crystal has UV cutoff wavelength around 315 nm as clear from Fig. 7(a), which also shows that the grown crystal has very low absorbance or high transmittance i.e. ∼80% in visible region which indicates its use in optical window for SHG laser radiation applications in blue region (Rao, 1974). The optical transparency of the grown crystal is higher than many other materials.

(a) Absorbance spectra and (b) optical band gap plot of BTSN.
Figure 7 (a) Absorbance spectra and (b) optical band gap plot of BTSN.

We have calculated the optical band gap of the grown crystal from absorbance data as it will give very profound and good insight about its optical performance and application in optoelectronic devices. Optical band gap (Egopt) of the BTSN crystal was calculated according to the following steps:

  1. Initially we have calculated optical absorption coefficient (α) from absorbance data using the following relation:

(2)
α = Absorbance d
where d is the thickness of the measured crystal (= 0.12 cm).

  1. Then Egopt was evaluated from the following relation:

(3)
( α h ν ) 1 r = A ( h υ - E g )
where r is an index that differentiates the process of optical absorption in the crystal and its value corresponds to r = 1/2 for direct allowed transitions, r = 3/2 for direct forbidden transitions, r = 2 for indirect allowed transitions, r = 3 for indirect forbidden transitions and A, h, ν are constants having their usual meanings. Egopt of BTSN crystal was calculated at all the transition values of r to be ensured that the grown crystal comes in direct band gap material category. The straight line was extrapolated in the graph for all r values to the x-axis () as shown in Fig. 7(b). The best fitting of the curve was observed at n = 1/2 which exemplifies to direct allowed transition in BTSN crystal and the value of band gap is found to be ∼3.63 eV. The high band gap of BTSN makes it a promising candidate for optoelectronic device applications (Shakir et al., 2009b; Shkir et al.; Shkir et al., 2012a, 2012b, 2014a, 2014b, 2014c, 2015c; Shkir, 2011).

Diffuse reflectance (DR) is a destructive technique for surface analysis using a mirror like reflection from the sample using an integrating sphere attachment. The ratio of the light scattered from an infinitely thick layer and from an ideal non-absorbing reference sample is measured as a function of the wavelength in diffused reflectance spectrum measurement. Illumination of the powdered samples by incident radiation leads to diffuse illumination of the samples and the incident light is partially absorbed and scattered. The scattered radiation, coming from the sample is collected in an integration sphere and detected (Weckhuysen and Schoonheydt, 1999). In general the Kubelka-Munk theory is used for diffuse reflectance spectra analysis from the weakly absorbing samples. DR of the titled crystal has been measured for the first time. For diffuse reflectance measurements we have crushed the grown crystals in homogeneous powder form and then loaded into the deep sample holder. The recorded DR spectra of BTSN crystal are shown in the Fig. 8(a). The measure data were used to evaluate the band gap.

(a) Diffuse reflectance and (b) optical band gap from DR data.
Figure 8 (a) Diffuse reflectance and (b) optical band gap from DR data.

At any wavelength the Kubelka-Munk equation is given by the following:

(4)
F ( R ) = ( 1 - R ) 2 2 R where F(R) is called Kubelka-Munk function and R is the absolute reflectance of the sample.

Eq. (2) can be written in terms of F(R) as follows:

(5)
α = absorbance t = F ( R ) t where the symbols are having their usual meanings and the value of t of the loaded circular sample is 1.5 mm.

For calculating the optical band gap from DR data the Eq. (2) can be written as follows:

(6)
( α h ν ) 1 r = F ( R ) h υ t 1 r = A ( h υ - E g ) where the symbols are having their usual meanings. The band gap of BTSN was calculated for all the transition values of r by same procedure as above and the graph between vs [F(R)/t]1/r is shown in Fig. 8(b). Best-fitted curve was observed at n = 1/2 which indicates that the grown BTSN crystal has direct allowed transition and the value of band gap is ∼3.6 eV. The band gap calculated from diffuse reflectance was found to be almost similar as calculated from UV–vis-NIR data.

Time dependent density functional theory (TD-DFT) is considered to be a reliable quantum chemical technique for studying the UV–Visible absorption spectra (corresponds to vertical electronic transitions) of different organic compounds (Jacquemin et al., 2004, 2005; Cossi and Barone, 2001). Furthermore, the comparison of calculated and experimental absorption wavelengths has also been made as per suggestion of respected reviewer. For instance, different TD-DFT methods have been applied to study the absorption spectrum of BTSN on its ground state optimized geometry in the framework of B3LYP/Gen level of theory. The experimental electronic absorption spectrum of title BTSN molecule is showing maximum absorption at ∼315 nm. The absorption spectra of BTSN at five different TD-DFT methods have been shown in Fig. 9 and its experimental spectrum is also illustrated through the inset of Fig. 9. It can be seen from Fig. 9 that all the TD-DFT methods (except TD-CAM-B3LYP) have reasonably reproduced the experimental absorption spectrum of BTSN molecule.

The calculated absorption spectra of BTSN at different levels of theory while inset is its experimental absorption spectrum and shape of grown crystal.
Figure 9 The calculated absorption spectra of BTSN at different levels of theory while inset is its experimental absorption spectrum and shape of grown crystal.

A careful analysis of Fig. 9 shows that the absorption energies calculated with TD-B3LYP and TD-M06 are closer to experimental wavelength, which represents their better accuracy in case of BTSN molecule. Calculations of optimized molecular orbital geometry of BTSN show that the visible absorption maxima are of the electron transition between frontier orbitals such as transition from HOMO to LUMO which will be explained in next Section 3.6. The maximum absorption band observed in UV region 315 nm (Exp.) and 300 nm (Theoretical) could be attributed to the redistribution of intramolecular electronic charge involving HOMO and LUMO orbitals.

3.6

3.6 Frontier molecular orbital (FMO) analysis

The frontier molecular orbitals (FMOs) play a very crucial role in the reactivity of any molecule. Among FMOs, the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) are very important. These FMOs of a molecule determine the way that it interacts with other species. The spatial distributions of the molecular orbitals for BTSN molecule have been shown in Fig. 10 as calculated at B3LYP/6-31G∗ level of theory. It can be seen from Figure that the HOMO and LUMO orbitals are localized on different components of BTSN molecule. The HOMO-LUMO transition indicates the redistribution of electronic density from nitrate ion fragment to main silver thiourea part. As the transition does not involve any significant intramolecular charge transfer within the BTSN molecule so the transition energy is relatively larger and lie in UV region, which give us an advantage of significant transparency. Based on energies of FMOs, we have also calculated the global reactivity parameters as given in Table 4. The global reactivity parameters play a very crucial role in the reactivity of any molecule and also determine the way that it interacts with other species.

FMOs involve in critical transition in BTSN molecule calculated at B3LYP/6-31G∗ level of theory.
Figure 10 FMOs involve in critical transition in BTSN molecule calculated at B3LYP/6-31G∗ level of theory.
Table 4 The calculated energy values of frontier molecular orbitals (FMOs) and energy difference at B3LYP/6-31G∗ level of theory for BTSN.
Orbitals a.u. eV
EHOMO −0.202 −5.497
EHOMO-1 −0.204 −5.551
EHOMO-2 −0.216 −5.878
ELUMO −0.031 −0.844
ELUMO+1 −0.018 −0.489
ΔEHOMO–LUMO 0.171 4.653
EHOMO-1–ELUMO 0.173 4.707
EHOMO-2–ELUMO 0.185 5.034
ΔEHOMO–1−LUMO+1 0.186 5.062
Hardness (η) 0.085 2.3265
Potential (μ) −0.116 −3.1515
Softness (S) 0.043 1.1814
Electronegativity (χ) 0.116 3.1515
Electrophilic index (ω) 0.078 2.1345

3.7

3.7 Density of states (DOS)

The total density of states (TDOS) and partial density of states (PDOS) are investigated for BTSN molecule using AOMix wave function analysis program (Gorelsky and Lever, 2001). The total density of states (TDOS) and partial density of states (PDOS) are calculated to understand the role of individual molecules into the bonding properties of BTSN molecule. We generate the total and partial density-of-states plots by dividing BTSN molecule into three fragments including Ag, nitrate and thiourea fragments as shown in Fig. 11. The TDOS and PDOS plots show population analysis per orbital and demonstrate the modest view of the makeup of the molecular orbitals in a certain energy range while PDOS plot shows percentage contribution of each group to each molecular orbital in the final molecule. A careful analysis of Fig. 11 shows that thiourea contributes more to the total number of states per interval of energy as compared with other fragments. The most of the high energy occupied states and lower energy unoccupied states around the band gap are composed of thiourea fragment. Similarly the contribution of Ag is significant over energy range of −9 to −8 eV and a further contribution is at lower conduction bands around −1 to 6 eV. Additionally, it can be seen that the contributions of nitrate are between −14 to −12 eV and −8 to −5 eV in valance bands; and 0 to −1 eV and 5 to 7 eV in conduction bands.

Total density of state (TDOS) and partial density of state (PDOS) plots for BTSN molecule.
Figure 11 Total density of state (TDOS) and partial density of state (PDOS) plots for BTSN molecule.

3.8

3.8 Static polarizability and first hyperpolarizability (β)

In our present investigation, we have calculated the electronic dipole moment, molecular polarizability and first hyperpolarizability. For a molecule, its dipole moment (μ) is defined as follows:

(7)
μ = μ x 2 + μ y 2 + μ z 2 1 / 2

The average polarizability ( α 0 ) can be calculated by following equations:

(8)
α o = 1 3 ( α xx + α yy + α zz )

For anisotropy of polarizability (Δα)

(9)
Δ α = 1 2 [ ( α xx - α yy ) 2 + ( α yy - α zz ) 2 + ( α zz - α xx ) 2 + 6 α xz 2 ]

Similarly, the amplitude of first static hyperpolarizability (β0) is defined as

(10)
β 0 = ( β x + β y + β z ) 1 / 2 where
(11)
β i = 3 5 ( β iii + β ijj + β ikk ) i , j , k = x , y , z ,

Similarly, the total first hyperpolarizability (βtot) can also be calculated using following equation:

(12)
β tot = ( β xxx + β xyy + β xzz ) 2 + ( β yyy + β yzz + β yxx ) 2 + ( β zzz + β zxx + β zyy ) 2 1 / 2

The second-order polarizability (β) is a third rank tensor that can be described by a 3 × 3 × 3 matrix. According to Kleinman symmetry (βxyy = βyxy = βyyx, βyyz = βyzy = βzyy,…likewise other permutations also take same value), the 27 components of the 3D matrix can be reduced to 10 components (Kleinman, 1962). The details of 10 components have been shown in Table 5. It is well acknowledged fact that the importance of polarizability and hyperpolarizabity of a molecular system is dependent on the electronic communication of two different parts of a molecule. The calculated values of dipole moment (a.u.), average polarizability (α0), anisotropy of polarizability (Δα) and hyperpolarizability (β) are given in Table 4.

Table 5 The calculated values of polarizability (α), hyperpolarizability (β) and dipole moment (μ) along their individual tensor components for BTSN.
Components a.u. esu (×10−24) Component a.u. esu (×10−30)
αxx 250 37.02 βxxx −437 −3.77
αxy −11 −1.69 βxxy −28 −0.24
αyy 232 34.45 βxyy −405 −3.50
αxz −1 −0.16 βyyy 83 0.72
αyz 10 1.53 βxxz −37 −0.32
αzz 211 31.27 βxyz 4 0.04
α0 231 34.24 βyyz −375 −3.24
Δα 613 90.94 βxzz −119 −1.03
μx −5.468 −13.89D βyzz −54 −0.47
μy 0.885 2.24D βzzz 218 1.89
μz 0.747 1.89D βtot 980 (580)b 8.47 (5.01)
μtot 5.589 14.20D β0 588 (348) 5.08 (3.00)
μtot (urea) 1.66 4.24D (4.56D)a βtot (urea) 43 0.37

For α, 1 a.u. = 0.1482 × 10−24 esu, for β, 1 a.u. = 0.008629 × 10−30 esu.

Values in parenthesis are with MP2 method.

Table 5 shows the calculated value of total dipole moment (14.20D) for BTSN molecule. The highest component of dipole moment is μx having a value of −13.89 D, where the negative sign implies that it is almost entirely directed along negative x-axis as shown in Fig. 12. Fig. 12 clearly indicates that the dipole moment of BTSN molecule is entirely directed from nitrate ion toward tetrakis(thiourea) silver(I) along the negative x-axis as shown in Fig. 10. As the electronic dipole moment is usually directed from negative to positive point charges, the Milliken population analysis indicates a significantly negatively charged Nitrate ion and positively charged silver(I) present along x-axis. A similarly but somewhat different analysis can also be performed for polarizability values i.e. average polarizability (α0), anisotropy of polarizability (Δα) and hyperpolarizability (βtot) of BTSN molecule, which have nonzero amplitudes of 34.24 × 10−24, 90.94 × 10−24 and 5.01 × 10−30 esu, respectively. The nonzero value of βtot shows that the titled molecule possesses microscopic first static hyperpolarizability. The first hyperpolarizability value of BTSN molecule is 27 times larger than that of urea as calculated in the present investigation at the same B3LYP/6-31G∗ level of theory, which shows that BTSN can be considered as a potential candidate for NLO applications having a further advantage of its good transparency in visible region as discussed in UV–Visible spectra of BTSN in next section.

The representation of partial charges (relative color format) and dipole moment vector calculated using Milliken population analysis at B3LYP/6-31G∗ level of theory for BTSN.
Figure 12 The representation of partial charges (relative color format) and dipole moment vector calculated using Milliken population analysis at B3LYP/6-31G∗ level of theory for BTSN.

3.9

3.9 Frequency dependent polarizability and first hyperpolarizability

In addition to static first hyperpolarizability values, we have also calculated the dynamic frequency dependent first hyperpolarizabilities that are corresponding to electric field induced second harmonic generation (EFISHG) β ω SHG ( - 2 ω ; ω , ω ) and the electro-optical Pockels effect (EOPE) β ω OP ( - ω ; 0 , ω ) . The EFISHG and EOPE are usually approximated as complement to the experimental frequency dependent first hyperpolarizability values. The frequency dependent coupled-perturbed Kohn-Sham (CPKS) method with the B3LYP functional has been applied to calculate the dynamic first hyperpolarizability values. In case of EFISHG experiments, the measurement provides information on the projection of the vector part of β on the dipole moment vectors as given by following equation

(13)
β ω ( - 2 ω ; ω , ω ) = β ω = 3 5 ζ x , y , z μ ζ β ζ μ where μ is the norm of dipole moment vector and μζ and βζ are the components of μ and β vectors. In similar way, the EOPE β ω OP ( - ω ; 0 , ω ) has been also calculated with different frequency input orientations. All the βω values have been given in electrostatic units (10−30 esu) within T-convention of reference (Willetts et al., 1992). In CPKS method, the matrices of CPKS equation are expanded in Taylor series of external dynamic electric field and are solved analytically order by order. According to experimental setup for EFISHG first hyperpolarizability (βω) measurement, the frequency dependent calculations are carried out using different optical wavelengths. Furthermore, we have also calculated the dynamic frequency dependent values of polarizability including isotropic polarizability and anisotropic polarizability for BTSN. In the present investigation, a number of optical wavelengths have been used to determine the effect of incident laser wavelengths on polarizabilities (αiso, and αiso), βω(−2ω; ω, ω) and βω(−ω; 0, ω) values. Three different optical wavelengths have been used in the present investigation including 487.8, 799.4 and 1064 nm. From Table 6, it can be seen that for frequency dependent dynamic polarizability both its isotropic and anisotropic values show a slight gradual increase with decrease in the optical wavelengths of laser. A similar trend can also be seen for the frequency dependent βω (−ω; 0, ω) and βω (−2ω; ω, ω) values and there is slight increase with increasing the frequency wavelength. The frequency dependent βω (−2ω; ω, ω) and βω (−ω; 0, ω) values are also found to be slightly larger than static first hyperpolarizability (β0), which represent that the dispersion effects on SHG and OP values are not very large for BTSN molecule. The frequency dependent βω (−2ω; ω, ω) and βω (−ω; 0, ω) amplitudes are ∼13 times larger than those of their respective amplitudes in urea molecule as calculated at the same level of theory.
Table 6 The calculated values of frequency dependent polarizabilities including isotropic and anisotropic polarizability along with dynamic EFISHG hyperpolarizability (βωSHG, in ×10−30 esu) and dynamic EOPE hyperpolarizability (βωOP in ×10−30 esu) for BTSN.
Frequency (nm) αiso (×10−24 esu) αaniso (×10−24 esu) βωSHG (×10−30 esu) βωOP (×10−30 esu)
487.8 39.35 6.70 3.76 4.15
799.4 37.25 6.24 3.57 4.22
1064 36.78 6.14 3.50 4.41
1064 (urea) 4.05 4.51 0.274 (0.45 ± 0.12)a 0.304
Experimental SHG value of Urea as calculated at ω = 1064 nm in water (Ledoux and Zyss, 1982).

3.10

3.10 Molecular electrostatic potential

To have molecular level understanding, we have calculated 3-D plots of molecular electrostatic potential of BTSN molecule as shown in Fig. 13. The MEP is the measurement of electrostatic potential on constant electron density surface. The 3-D plots of MEP surface overlap on the top of total energy density. The MEP is helpful property to investigate the reactivity of molecular species by predicting that the approaching nucleophile is attracted to a positive region of molecule. In MEP plot, the maximum positive region that is preferred site for nucleophilic attack is indicated as blue color, while a maximum negative region is preferred site for electrophilic attack is indicated as red surface. The MEP of BTSN has been drawn in Fig. 13 to get simultaneous information about its molecular size, shape along with its positive, negative and neutral electrostatic potential regions in terms of color grading. The red color on MEP surface shows highest negative potential while blue color shows highest positive potential. The green color represents the slight positive or neutral potential over the thiourea atoms indicating their non-reactive nature toward any substitutions. The MEP is a real physical property, which is an observable in experimental diffraction method.

Molecular electrostatic potential (MEP) plot of BTSN molecule.
Figure 13 Molecular electrostatic potential (MEP) plot of BTSN molecule.

3.11

3.11 Dielectric studies

As per the concern over storage and dissipation of electric and magnetic energy the investigation of dielectric properties of any material is a necessary task (Von Hippel, 1954; von Hippel and Morgan, 1955). It is well known that the dielectric properties play a key role in explaining the various phenomena in electronics, optics, and solid-state physics. There is clear evidence in the literature that the dielectric and electro-optical properties of a crystal (Boomadevi and Dhanasekaran, 2004; Shakir et al., 2010a, 2010b, 2010c, 2010d) are interlocked, especially for non-conducting materials. Generally, a larger voltage is required to the material which has high value of dielectric constant needs in order to polarize the dipoles, due to which it may suffer changes in its refractive index. But as per our study on the title material it has low value of dielectric constant as well as dielectric loss, which eliminates the need of poling during maintaining the refractive index.

The grown crystals were subjected to dielectric measurements in the high frequency range from 3 kHz to 10 MHz under dark and high intensity UV light and the plots of variation of dielectric constant ( ε r ), dielectric loss (tan δ) are shown in Fig. 14(a) and (b). As per the available literature the value of dielectric constant founds to be higher in the lower-frequency region and then it decreases and becomes stable at higher frequencies, the reason behind this behavior is well explained (Shakir et al., 2009a, 2010b, 2010c, 2010d; Xue and Kitamura, 2002). By keeping this reason in mind we have done our dielectric measurement in higher frequency range. As clear from Fig. 14(a) the values of ε r are found to be stable in the entire testing range and based on the fact that the electric field, the dipole does not follow the alternating field beyond a certain frequency (Batra et al., 2005; Hill, 1969). The value of dielectric constant was found to be almost doubled under UV light which suggests its use in UV region to get good out from the device. The value of dielectric loss (Fig. 14(b)) also follows the similar path as that of dielectric constant. The value of dielectric loss depends on several factors such as size of crystal, defects and can be controlled as well (Przeslawski et al., 1995; Shakir et al., 2010b). As the value of dielectric is found to be very low which suggests the grown crystals are free from any major defects and high crystalline perfection, therefore it may be useful in electro-optic modulation applications.

Plots of variation of (a) dielectric constant and (b) dielectric loss measured on BTSN crystal.
Figure 14 Plots of variation of (a) dielectric constant and (b) dielectric loss measured on BTSN crystal.

4

4 Conclusions

We have successfully applied a dual approach to study various key parameters of a thiourea metal complex i.e. BTSN. We have successfully grown the single crystals of bis(thiourea) silver nitrate (BTSN) by slow evaporation solution technique in an aqueous solution. The grown crystals were subjected to various characterizations to analyze its key opto-electrical properties. Crystal structure and lattice parameters were calculated and found in well agreement with earlier reported experimental findings. The surface morphology and elemental composition were determined by SEM/EDXS studies, which show that the grown crystals are of good quality. The vibrational analysis was carried out using experimental and calculated FT-Raman and FT-IR spectra. The grown crystal was cut and polished for UV–vis.-NIR and diffuse reflectance studies as well as various other optical parameters were also calculated. The grown crystals were found to be highly transparent in the entire testing range, which warrant their use in optoelectronic devices. The cutoff wavelength of BTSN was observed at ∼315 nm. Their optical band gaps were calculated from absorbance and diffuse reflectance data and found to be 3.63 eV and 3.6 eV respectively, which are very close to each other. The dielectric study also shows that the grown crystals are having good crystalline perfection and the value of dielectric constant was found to be almost double under the UV light. In addition to the above, we have used different computational methods such as B3LYP/6-31G∗, M06/6-31G∗ and MP2/6-31G∗ and the stable ground state molecular geometry of BTSN was optimized. Various theoretically calculated properties such as geometrical, spectroscopic, linear and nonlinear optical were compared with our as well as reported experimental results and found in strong correlation. The values of static (βtot) and dynamic first hyperpolarizability (βω) were found to be 23 and 13 times greater than a typical prototype urea molecule, respectively, calculated theoretically. Further, the MEP, TDOS and PDOS were also calculated and analyzed, which offer in-depth understanding about numerous key features of individual molecular components in BTSN complex. From above results it is clear that BTSN has good linear and nonlinear optical properties, which makes it very useful in various optoelectronic device applications.

Acknowledgement

The authors would like to express their gratitude to King Khalid University, Saudi Arabia for providing administrative and technical support.

References

  1. Abbas, H., Shkir, M., AlFaify, S., 2019. Arab. J. Chem. 12(8), 2336–2346.
  2. , . Acta Crystallogr. Sect. B: Struct. Sci.. 2002;58:380-388.
  3. , , , . J. Mol. Struct. (Thoechem.). 2008;864:20-25.
  4. , , . Cryst. Res. Technol.. 2001;36:1231-1237.
  5. , , , , , . Cryst. Res. Technol.. 2005;40:757-760.
  6. , , . J. Cryst. Growth. 2004;261:70-76.
  7. , , , . Inorg. Chem.. 2009;48:3185-3197.
  8. , , , . Polyhedron. 1999;18:2665-2671.
  9. , , , , , , , , , . J. Inorg. Biochem.. 2004;98:1307-1314.
  10. , , . J. Chem. Phys.. 2001;115:4708-4717.
  11. , . J. Chem. Soc., Faraday Trans. 1: Phys. Chem. Condens. Phases. 1986;82:3415-3430.
  12. Frisch, M.J., Trucks, G.W., Schlegel, H.B., Scuseria, G.E., Robb, M.A., Cheeseman, J.R., Scalmani, G., Barone, V., Mennucci, B., Petersson, G.A., Nakatsuji, H., Caricato, M., Li, X., Hratchian, H.P., Izmaylov, A.F., Bloino, J., Zheng, G., Sonnenberg, J.L., Hada, M., Ehara, M., Toyota, K., Fukuda, R., Hasegawa, J., Ishida, M., Nakajima, T., Honda, Y., Kitao, O., Nakai, H., Vreven, T., Montgomery, J.A., Peralta, J.E., Ogliaro, F., Bearpark, M., Heyd, J.J., Brothers, E., Kudin, K.N., Staroverov, V.N., Kobayashi, R., Normand, J., Raghavachari, K., Rendell, A., Burant, J.C., Iyengar, S.S., Tomasi, J., Cossi, M., Rega, N., Millam, J.M., Klene, M., Knox, J.E., Cross, J.B., Bakken, V., Adamo, C., Jaramillo, J., Gomperts, R., Stratmann, R.E., Yazyev, O., Austin, A.J., Cammi, R., Pomelli, C., Ochterski, J.W., Martin, R.L., Morokuma, K., Zakrzewski, V.G., Voth, G.A., Salvador, P., Dannenberg, J.J., Dapprich, S., Daniels, A.D., Farkas, Ö., Foresman, J.B., Ortiz, J.V., Cioslowski, J., Fox, D.J., 2009. Gaussian 09, Revision B.01. Wallingford CT.
  13. , , , , , , , , , , . Jpn. J. Appl. Phys.. 2006;45:8676.
  14. , , . J. Organomet. Chem.. 2001;635:187-196.
  15. , , . Spectrochim. Acta Part A Mol. Biomol. Spectrosc.. 2014;122:130-141.
  16. , . Dielectric Properties and Molecular Behaviour. Van Nostrand Reinhold; .
  17. , , , , , , . J. Chem. Phys.. 2004;121:1736-1743.
  18. , , , . Chem. Phys. Lett.. 2005;410:254-259.
  19. , , , , . Physica B. 2006;371:1-4.
  20. , , . Adv. Mater.. 1999;11:1147-1151.
  21. , , , . J. Cryst. Growth. 2016;440:113-114.
  22. , , , , . J. Cryst. Growth. 2004;269:565-569.
  23. , . Phys. Rev.. 1962;126:1977.
  24. , , , , , , . Spectrochim. Acta Part A Mol. Biomol. Spectrosc.. 2009;73:263-267.
  25. , , . J. Chem. Phys.. 1982;73:203-213.
  26. Marder, S.R., Sohn, J.E., Stucky, G.D., 1991. DTIC Document.
  27. , , , , , , , , , , . J. Comput. Chem.. 2015;36:118-128.
  28. , , , , , , . J. Mater. Chem.. 2013;1:5439-5449.
  29. , , , , , . PCCP. 2010;12:4791-4799.
  30. , , , , , . J. Am. Chem. Soc.. 2009;131:11833-11840.
  31. , , , , , , , . Dalton Trans.. 2013;42:15053-15062.
  32. , , , . Acta Crystallogr. A. 1965;18:618-623.
  33. , , , . Non-linear Optical Properties of Matter. Springer; .
  34. , , , , , . Prog. Cryst. Growth Charact. Mater.. 1991;22:19-51.
  35. , , , . Acta Crystallogr. Sect. C: Cryst. Struct. Commun.. 1997;53:565-568.
  36. , , , . Solid State Commun.. 1995;96:195-198.
  37. , . Transition-metal oxides. Solid state chem. 1974
    [Google Scholar]
  38. Saleh, B. E. A., 1991. MC Teich Fundamentals of photonics.
  39. , . Hydrogen Bonding: A Theoretical Perspective (Topics in Physical Chemistry). Oxford: Oxford University Press; .
  40. , , , , , , , . Mater. Chem. Phys.. 2007;103:153-157.
  41. , , , , , . Mater. Sci. Eng., B. 2010;172:9-14.
  42. , , , , , , . J. Appl. Crystallogr.. 2010;43:491-497.
  43. , , , , , . J. Cryst. Growth. 2009;311:3871-3875.
  44. , , , , , . Solid State Commun.. 2009;149:2047-2049.
  45. , , , , . Mater. Chem. Phys. 2010;120:566-570.
  46. , , , , , , . J. Cryst. Growth. 2010;312:3171-3177.
  47. , , , , . Appl. Phys. Lett.. 2009;95 (pp. 252902–252902-252903)
  48. , . Chalcog. Lett.. 2011;8:435-440.
  49. , , , , , , . Nanosci. Nanotechnol. Lett.. 2012;4:405-408.
  50. , , . Spectrochim. Acta Part A Mol. Biomol. Spectrosc.. 2014;118:172-176.
  51. , , , . J. Phys. Chem. Solids. 2012;73:1309-1313.
  52. , , , , , . J. Phys. Chem. Solids. 2014;75:959-965.
  53. , , , , , . J. Cryst. Growth. 2014;391:104-110.
  54. Shkir, M., Muhammad, S., AlFaify, S. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy.
  55. , , , . Spectrochim. Acta Part A Mol. Biomol. Spectrosc.. 2015;143:128-135.
  56. , , , , , , , , . RSC Adv.. 2015;5:87320-87332.
  57. , , , , , . Spectrochim. Acta Part A Mol. Biomol. Spectrosc.. 2015;137:432-441.
  58. , , , , , , , , . Opt. Mater.. 2014;36:675-681.
  59. , , , , , , . Spectrochim. Acta Part A Mol. Biomol. Spectrosc.. 2014;118:603-613.
  60. , , . J. Cryst. Growth. 2016;440:110-112.
  61. , . Dielectric Materials and Applications. Artech House on Demand; .
  62. , , . J. Electrochem. Soc.. 1955;102:68C.
  63. , , . J. Chem. Phys.. 1985;82:284-298.
  64. , , . Catal. Today. 1999;49:441-451.
  65. , , , , . J. Chem. Phys.. 1992;97:7590-7599.
  66. , , . Solid State Commun.. 2002;122:537-541.
  67. , , , , , , , . Inorg. Chim. Acta. 2003;343:361-365.
  68. , , , . J. Chem. Phys.. 1984;81:4160-4167.

Appendix A

Supplementary material

Supplementary data associated with this article can be found, in the online version, at http://dx.doi.org/10.1016/j.arabjc.2016.06.016.

Appendix A

Supplementary material

Supplementary data 1

Supplementary data 1

Show Sections