Ancient Turkey Red Dye Alizarin Emerges as Safe, Sustainable Colorant for Toys

Alizarin, the brilliant orange-red dye that once colored the coveted Turkey Red textiles of the Ottoman Empire and fueled one of the most consequential breakthroughs in industrial chemistry, is now being reexamined through the lens of twenty-first-century quantum computing. A team of researchers from the Centre of Toy Science at Children’s Research University in Gandhinagar, Gujarat Technological University, Saurashtra University, and collaborating institutions has carried out a comprehensive computational characterization of the molecule, asking a deceptively simple question: could this ancient natural dye replace the synthetic colorants that currently saturate the toy industry? Their answer, published in Discover Chemistry, is a cautiously enthusiastic yes, backed by a dense lattice of quantum chemical calculations, simulated spectra, and toxicological predictions.

The stakes are higher than they might appear. Children are not simply small adults when it comes to chemical exposure. Studies cited by the research team indicate that children can ingest up to 24 milligrams of toy material per kilogram of body weight per day through mouthing behavior, making any leachable dye a genuine safety concern. Conventional synthetic colorants, particularly azo dyes, have long been scrutinized because they can release carcinogenic aromatic amines under reducing conditions. Regulatory frameworks such as the European EN 71-9/10/11 standards and the American ASTM F963 impose strict limits on toxic colorants in playthings, and the industry is actively hunting for alternatives that deliver vivid color without the toxicological baggage. Alizarin, with its centuries-long history of human use in cosmetics and consumer products, its natural origin in the root of the madder plant Rubia tinctorum, and its biodegradability, is an obvious candidate. What has been missing is a rigorous, quantitative description of the molecule’s electronic and optical behavior that could anchor such a substitution in hard science.

The research team, led by N. H. Vaosya and Rajveersinh Zala with supervision from Bharat Kataria and J. H. Markna, turned to density functional theory and its time-dependent extension, the workhorse methods of modern computational chemistry. Using the ORCA 6.1.1 software package with the B3LYP functional and the def2-SVP basis set developed by Weigend and Ahlrichs, they optimized the geometry of the alizarin molecule, which contains fourteen carbon atoms, eight hydrogens, and four oxygens arranged across three fused aromatic rings. The choice of basis set was deliberate: def2-SVP offers geometries within 0.01 angstroms of experimental structures while keeping the calculation tractable, requiring roughly twelve hours on a standard workstation for the 204 basis functions needed to describe the molecule. The optimization converged within eighteen cycles, and the resulting structure matched experimental X-ray crystallographic data with a root-mean-square deviation of just 0.023 angstroms, confirming a planar anthraquinone core with near-symmetric geometry.

That planarity is not a cosmetic detail. The flat, extended pi-conjugated system is precisely what makes alizarin a dye at all. Delocalized pi-electrons flowing across the three-ring skeleton create a small energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital, and the calculations placed that HOMO-LUMO gap at 1.79 electron volts, with the HOMO sitting at minus 5.53 eV and the LUMO at minus 3.74 eV. This is squarely in the range needed for visible-light absorption, which typically requires gaps between 1.5 and 3.0 eV. The spatial separation of the two orbitals is equally telling: the HOMO concentrates on the electron-rich dihydroxy-substituted ring, while the LUMO spreads over the electron-poor carbonyl regions. This arrangement creates a classic push-pull architecture, with hydroxyl groups donating electron density and carbonyl groups accepting it, setting up a strong intramolecular charge transfer when the molecule absorbs a photon.

Time-dependent DFT calculations of the first thirty singlet-singlet transitions revealed exactly how that charge transfer plays out optically. The most intense visible absorption appears at 483.8 nanometers with an oscillator strength of 0.130, corresponding to a pi-to-pi-star transition from HOMO-1 to LUMO+1. Electron density difference maps showed roughly 0.38 electrons of charge shifting from the dihydroxy ring toward the carbonyl region upon excitation, a hallmark of charge-transfer character that explains both the molecule’s saturated color and its known solvatochromism, the way its absorption shifts by 15 to 25 nanometers in more polar solvents. The computed absorption maximum lands within three nanometers of experimental values measured in ethanol, and the derived molar extinction coefficient of approximately 2.95 times ten to the third liters per mole per centimeter agrees with measured anthraquinone data within eight percent. Translating the absorption into perceived color, the team calculated CIE chromaticity coordinates of 0.62 and 0.35, placing the dye firmly in the orange-red region that toy designers prize for its visual appeal.

The simulated spectra went further, dissecting the molecule’s vibrational and magnetic behavior. The calculated infrared spectrum reproduced the signature carbonyl stretch at 1668 wavenumbers, notably redshifted by 15 to 30 wavenumbers relative to isolated quinone carbonyls because of intramolecular hydrogen bonding between the ortho-hydroxyl groups and adjacent carbonyl oxygens. Those hydrogen bonds, with oxygen-oxygen distances of 2.56 to 2.59 angstroms and bond energies estimated at 12 to 15 kilojoules per mole, rigidify the molecule and reduce the HOMO-LUMO gap by about 0.15 eV, deepening the red color. Overall, the computed IR frequencies matched literature values with a correlation coefficient of 0.998 and a root-mean-square deviation of 9.8 wavenumbers. Gauge-including atomic orbital NMR calculations predicted carbonyl carbon signals separated by 8 to 10 parts per million, a diagnostic fingerprint that distinguishes alizarin from its isomer quinizarin, and a linear relationship between calculated carbon charges and chemical shifts emerged with a correlation coefficient of 0.92, consistent with Karplus-Pople shielding theory.

Mulliken population analysis added a quantitative picture of the molecule’s chemical personality. Oxygen atoms carry charges between minus 0.161 and minus 0.170 elementary charge units, while the phenolic hydrogens sit at plus 0.174 to plus 0.180, reflecting their acidity with a calculated pKa of 6.2. This polarization, producing a molecular dipole moment of 4.09 Debye, underpins two properties with direct relevance to interactive toys: alizarin’s ability to chelate metal cations through its catechol-like 1,2-dihydroxyanthraquinone motif, and its dramatic halochromic response. When deprotonated in basic conditions, the absorption band shifts by 40 to 60 nanometers toward 520 to 540 nanometers, flipping the color from orange-red to purple-red. The authors suggest this pH sensitivity and metal-binding behavior could enable educational and color-changing toys, and the same catechol binding motif has already made alizarin a candidate sensitizer for titanium dioxide photoanodes in dye-sensitized solar cells.

On the safety front, the team employed the U.S. Environmental Protection Agency’s Toxicity Estimation Software Tool, version 5.1.2, using consensus quantitative structure-activity relationship modeling. The predictions are encouraging: an acute oral LD50 exceeding 5000 milligrams per kilogram, no predicted mutagenicity, no skin sensitization potential, and existing regulatory approval by the FDA and EU authorities. Aquatic toxicity modeling predicted a fathead minnow 96-hour LC50 of 5.19 on the negative logarithmic molar scale, indicating moderate rather than acute hazard, with the model showing improved reliability for structurally similar compounds. The authors are candid about the limits of these in silico results, however. They emphasize that computational predictions are screening tools, not regulatory approvals, and that experimental validation including cytotoxicity assays, dermal irritation studies, and migration testing under EN 71-10/11 protocols remains essential before any toy bearing alizarin color reaches a shelf.

The study also flags genuine open questions. The team did not model how strongly alizarin binds to polymer matrices such as polyethylene or polypropylene, a critical gap given that poor dye-substrate affinity drives leaching and wastewater contamination. Photostability is another concern: experimental work has shown that alizarin degrades under ultraviolet light and chlorination into products including phthalic anhydride and benzoic acid derivatives, and transformation products can sometimes rival or exceed their parent compounds in toxicity. The researchers call for molecular dynamics simulations of dye-polymer interactions, systematic photodegradation studies identified through mass spectrometry, and full toxicological testing of degradation products. Even with those caveats, the convergence of a 1.79 eV electronic gap, a strong visible absorption at 483.8 nanometers, exceptional structural stability from intramolecular hydrogen bonding, and a favorable predicted safety profile positions alizarin as a scientifically grounded, environmentally sustainable alternative to synthetic azo dyes. A molecule first synthesized by Carl Graebe and Carl Liebermann in 1868, which helped launch the synthetic dye industry, may now help unwind some of its chemical excesses, one children’s toy at a time.

Subject of Research: DFT and TD-DFT computational analysis of the natural dye alizarin's photophysical properties and toxicological safety for toy applications

Article Title: Computational insights into Alizarin for toy applications using DFT calculations of photophysical properties and spectroscopic safety

Article References: Vaosya, N. H., Zala, R., Jadav, G., Katariya, D. K., Khunt, R., Kataria, B., & Markna, J. H. (2026). Computational insights into Alizarin for toy applications using DFT calculations of photophysical properties and spectroscopic safety. Discover Chemistry, 3(1), Article 482. https://doi.org/10.1007/s44371-026-00939-4

Image Credits: AI Generated

DOI: 10.1007/s44371-026-00939-4

Keywords: alizarin, DFT, TD-DFT, toy safety, natural dyes, HOMO-LUMO gap, anthraquinone, intramolecular hydrogen bonding, QSAR toxicity prediction, sustainable manufacturing, UV-Vis spectroscopy, NMR chemical shifts