গাণিতিক রসায়ন

কোয়াণ্টাম বলবিজ্ঞান-এর ইতিহাসে প্রাথমিক আবিষ্কার এবং মৌলিক নীতিগুলির উপর ভিত্তি করে ১৯২৭ সালে রসায়ন-এ প্রথম তাত্ত্বিক গণনা করেছিলেন ওয়াল্টার হীটলার ও ফ্রীত্জ লন্ডন। গাণিতিক কোয়াণ্টাম রসায়নের প্রাথমিক অগ্রগতিতে যেসব বই গুরুত্বপূর্ণ ভূমিকা পালন করেছিল সেগুলি হল- লিনাস পাওলিং ও ই. ব্রাইট উইলসন-এর "Introduction to Quantum Mechanics – with Applications to Chemistry" (১৯৩৫), আইরিং, ওয়াল্টার ও কিমবল-এর "Quantum Chemistry" (১৯৪৪), হীটলার-এর "Elementary Wave Mechanics – with Applications to Quantum Chemistry" (১৯৪৫), এবং পরবর্তীকালে কুলসন-এর "Valence" (১৯৫২)। এগুলোর প্রত্যেকটাই পরবর্তী দশকগুলিতে রসায়নবিদদের জন্য প্রাথমিক উল্যেক্ষ্য হিসেবে কাজ করেছিল।

১৯৪০ সালে দক্ষ কম্পিউটার প্রযুক্তির ক্রমবিকাশের সাথে সাথে জটিল পারমাণবিক সিস্টেমে তরঙ্গ-সমীকরণের বিশদ সমাধান সম্ভবপর হতে থাকে। পঞ্চাশের দশকের গোড়ার দিকে প্রথম সেমি-এম্পিরিকাল পারমাণবিক অর্বিটাল গণনা করা হয়। প্রথম ডিজিটাল কম্পিউটারগুলির প্রধান ব্যবহারকারী ছিলেন তাত্ত্বিক রসায়নবিদরা। যুক্তরাজ্যে এই ঘটনার বিস্তৃত বিবরণ দিয়েছেন স্মিথ ও সাটক্লিফ। ১৯৫৬ সালে স্লেটার-অর্বিটাল ভিত্তি-সেট ব্যবহার করে MIT তে প্রথমবার দ্বিপরমাণুক অণুর এব-ইনিশিও হারট্রি-ফক গণনা করা হয়। ১৯৬০ সালে দ্বিপরমাণুক অণুর লঘিষ্ঠ ভিত্তি-সেটে যথাযথ বিশ্লেষণ ও বৃহত্তর ভিত্তি-সেটে প্রথম গণনা প্রকাশ করেন যথাক্রমে রান্সিল ও নেসবেট। গাউসিয়ান অর্বিটাল ব্যবহার করে বহুপরমাণুক অণুর প্রথম গণনা করা হয় ১৯৫০ সালের শেষের দিকে। গাউসিয়ান অর্বিটাল ব্যবহার করে প্রথম কনফিগারেশন ইন্টারেক্শন গণনা হয় ক্যামব্রিজের EDSAC কম্পিউটারে ১৯৫০ সালে, করেছিলেন বয়েজ ও অন্যান্যরা। ১৯৭১ সালে এব-ইনিশিও গণনার জে গ্রন্থ-পঞ্জি প্রকাশিত হয় তাতে বৃহত্তম অণু ছিল ন্যাপথ্যালিন ও অ্যাজুলিন। অন্যান্য অনেক এব-ইনিশিও গণনার সারাংশ প্রকাশ করেছেন শেফার।

১৯৬৪ সালে বিউটাডাইন ও বেঞ্জিন থেকে শুরু করে ওভালিন পর্যন্ত নানা জটিল ও সহজ অণুর হুকেল-পদ্ধতি গণনা (সংযুক্ত হাইড্রোকার্বন সিস্টেমের π-আণবিক অর্বিটালে ইলেক্ট্রনিক শক্তি নির্ণয়ের জন্য সরলীকৃত পারমাণবিক অর্বিটালের রৈখিক সংযোগ (LCAO) পদ্ধতি) করা হয় বার্কলি ও অক্সফোর্ডের কম্পিউটারে। ১৯৬০ সালে এই পদ্ধতির জায়গায় CNDO-র মত সেমি-এম্পিরিকাল পদ্ধতি ব্যবহার শুরু হয়।

১৯৭০ সামের শুরুর দিকে আণবিক অর্বিটালের এব-ইনিশিও গণনা দ্রুততর করার জন্য ATMOL, গাউসিয়ান, IBMOL, POLYAYTOM ইত্যাদি দক্ষ এব-ইনিশিও কম্পিউটার প্রোগ্রাম ব্যবহার শুরু হয়। এই চারটির মধ্যে কেবল গাউসিয়ানের উন্নত সংস্করণ এখনো ব্যবহৃত হয়, যদিও এখন অন্য অনেক প্রোগ্রাম ব্যবহারের চল হয়েছে। এই একই সময়ে MM2 বলক্ষেত্রের মত আণবিক বলবিজ্ঞান পদ্ধতির উন্নতি ঘটে নরম্যান এলিংগারের হাত ধরে।

'গাণিতিক রসায়ন' শব্দবন্ধের প্রথম প্রয়োগ পাওয়া যায় সিডনি ফার্নবাক ও আব্রাহাম হাস্কেল টাউব এর ১৯৭০ সালে প্রকাশিত বই "Computers and Their Role in the Physical Sciences"-এ। তারা লিখেছেন, "...তাই মনে হচ্ছে 'গাণিতিক রসায়ন' ধীরে ধীরে বাস্তবে পরিণত হচ্ছে।" সত্তরের দশকে গাণিতিক রসায়ন নামের এই নতুন বিদ্যায় বিবিধ পদ্ধতি প্রয়োগ হতে থাকে। "Journal of Computational Chemistry" প্রথম প্রকাশ হয় ১৯৮০ সালে।

গাণিতিক রসায়নশাস্ত্রে বেশ কিছু নোবেল পুরষ্কার এসেছে যার মধ্যে সবচেয়ে উল্লেখযোগ্য হল ১৯৯৮ ও ২০১৩। "ডেন্সিটি-ফাংশনাল তত্ত্বের উদ্ভাবনের জন্য" ওয়াল্টার কোন, এবং "কোয়ান্টাম রসায়নের নতুন গাণিতিক পদ্ধতি উদ্ভাবনের জন্য" জন পপ্লে ১৯৯৮ সালে রসায়নে নোবেল পুরষ্কার পান। ২০১৩ সালে "জটিল রাসায়নিক সিস্টেমের মাল্টিস্কেল মডেলের জন্য" রসায়নে নোবেল পান মার্টিন কারপ্লাস, মাইকেল লেভিট ও এরিয়া ওয়ারশেল।

'তাত্ত্বিক রসায়ন' শব্দবন্ধ ব্যবহার করা যায় রসায়নের গানিতিক বিশ্লেষণধর্মী ব্যাখ্যার জন্য, অন্যদিকে গাণিতিক পদ্ধতি যখন কম্পিউটারে স্বয়ংক্রিয় প্রয়োগের উপযুক্ত হয় তখন তাকে 'গাণিতিক রসায়ন' বলে। তাত্ত্বিক রসায়নে রসায়নবিদরা, ভৌতবিজ্ঞানীরা এবং গণিতবিদেরা পারমাণবিক ও আণবিক ধর্ম ও রাসায়নিক বিক্রিয়ার বিক্রিয়াপথ নির্ণয়ের জন্য অ্যালগরিদম ও কম্পিউটার প্রোগ্রাম রচনা করেন। অন্যদিকে গাণিতিক রসায়নবিদরা সেই কম্পিউটার প্রোগ্রাম ও পদ্ধতিসমূহ প্রয়োগ করেন নির্দিষ্ট রাসায়নিক সিস্টেমে।

গাণিতিক রসায়নের দুরকমভাবে প্রয়োগ হয়ঃ

• পরীক্ষাগারে সংশ্লেষণের জন্য প্রাথমিক যৌগ নির্ণয় করা বা পরীক্ষালব্ধ তথ্য অনুধাবন করা, যেমন বর্ণালীর পীক-এর স্থান ও উৎস অনুসন্ধান করা।
• পরীক্ষার মাধ্যমে নিরীক্ষিত হয়নি এমন অজানা অণু বা বিক্রিয়াপথের অনুসন্ধান করা।

এভাবে গাণিতিক রসায়ন পরীক্ষামূলক রসায়নবিদের সাহায্যও করতে পারে, বা পরীক্ষামূলক রসায়নবিদকে নতুন রাসায়নিক বস্তু সন্ধানের দিকে চালিত-ও করতে পারে।

গাণিতিক রসায়নের বেশ কিছু আলাদা আলাদা ক্ষেত্র আছে, যেমনঃ

• নিউক্লিয়াস এর স্থান পরিবর্তন করে করে শক্তি-তল নির্ণয়ের মাধ্যমে আন্তর-পারমাণবিক বল সিমুলেট করে বা যথার্থ কোয়াণ্টাম রসায়ন পদ্ধতির মাধ্যমে অণুর গঠন নির্ধারন করা।
• কোনো রাসায়নিক বস্তুর ব্যাপারে তথ্য সংগৃহীত রাখা ও সন্ধান করা (দেখুন রাযায়নিক তথ্যভাণ্ডার)।
• আণবিক গঠন ও রাসায়নিক ধর্মের মধ্যে সংগতি স্থাপন করা (quantitative structure–activity relationship (QSAR))।
• রাসায়নিক যৌগের দ্রুত ও দক্ষ সংশ্লেষণের জন্য গণনা করা।
• গাণিতিকভাবে এমন অণু তৈরি করা যা অন্য কোনো অণুর সাথে পূর্বনির্ধারিত উপায়ে বিক্রিয়া করে (যেমন ঔষধ ডিজাইন ও অনুঘটন)।

নিখুঁত ও অভ্রান্ত শব্দগুলি এখানে প্রয়োগ হয় না কারণ রসায়নে খুব কম বিষয়্ই ত্রুটি ছাড়া গণনা করা যায়। কিন্তু রসায়নের সব বিষয়েরই গুণগত ও আসন্ন মানগত গাণিতিক বর্ণনা দেওয়া সম্ভব।

যেহেতু নিউক্লিয়াস ও ইলেক্ট্রন নিয়ে অণু গঠিত, তাই কোয়াণ্টাম বলবিজ্ঞানের পদ্ধতিগুলি এখানে প্রযোজ্য। গাণিতিক রসায়নবিদরা অনেকসময় অ-আপেক্ষিক শ্রোডিঙার সমীকরণের সাথে আপেক্ষিকতা সংক্রান্ত সংশোধন যগ করে সমাধান করেন, যদিও সম্পূর্ণরূপ আপেক্ষিক ডিরাক সমীকরণ সমাধানের দিকে কিছুটা অগ্রগতি হয়েছে। খাতায় কলমে সিস্টেম অনুযায়ী শ্রোডিঙার সমীকরণ সময়-নির্ভর ও সময়-অনির্ভর দুইরকম ক্ষেত্রেই সমাধান করা যায়, কিন্তু বাস্তবে খুব ছোটো সিস্টেম ছাড়া সেটা সম্ভব নয়। এরকম ভাবে অভ্রান্ততা ও গাণিতিক জটিলতার মধ্যে ভারসাম্য আনার জন্য প্রচুর সংখ্যক আসন্ন পদ্ধতি তৈরি হয়েছে।

গাণিতিক জটিলতা বাড়ালে গণনার যথার্থতা সবসময়ই বৃদ্ধি পায়। সম্পূর্ণ-আপেক্ষিক গণনার জটিলতার কারণে বহু-ইলেক্ট্রন এব-ইনিশিও মডেলে বেশ কিছু ত্রুটি থাকতে পারে। এই অক্ষমতার কারণে অবস্থান্তর ধাতু ও তাদের অনুঘটক ক্রিয়ার মত উচ্চ পারমাণবিক ভর-সম্পন্ন পরমাণুর সাথে অণুর ক্রিয়ার বিশ্লেষণ কঠিন হয়ে ওঠে। বর্তমানে গাণিতিক রসায়নে যেসব অ্যালগরিদম আছে সেগুলি ৪০ ইলেক্ট্রন অব্দি ধারণকারী অণুর ধর্ম জটিলতা ছাড়া সূক্ষ্মভাবে গণনা করতে পারে।শক্তির মানে ত্রুটির পরিমাণ হতে পারে সামান্য কিছু কিলোজুল/মোল এর থেকেও কম। আণবিক গঠনের ক্ষেত্রে বন্ধন-দূরত্ব কয়েক পিকোমিটার ও বন্ধন-কোণ ০.৫ ডিগ্রির মধ্যে নির্নয় করা যায়। এর থেকে বেশি ইলেক্ট্রন আছে এরকম বৃহত্তর অণুর গণনার জন্য ডেন্সিটি-ফাংশনাল তত্ত্বের মত আসন্ন পদ্ধতি ব্যবহার করা হয়।

জৈবরাসায়নিক বিক্রিয়ার মত জটিল বিক্রিয়ায় এসব আসন্ন পদ্ধতি প্রয়োগের উপযুক্ত কিনা এবিষয়ে বিতর্ক আছে। বৃহত অণুর বিশ্লেষণে সেমি-এম্পিরিকাল আসন্ন পদ্ধতি ব্যবহার করা যায়। অত্যন্ত বড় অণুর ক্ষেত্রে আণবিক বলবিজ্ঞান (Molecular Mechanics (MM)) নামক এক চিরায়ত বলবিজ্ঞান পদ্ধতিরও ব্যবহার আছে। QM-MM পদ্ধতিতে বৃহত সিস্টেমের ক্ষুদ্র অংশগুলিকে কোয়াণ্টাম বলবিজ্ঞানের মাধ্যমে ও বাকিটাকে আসন্নভাবে গণনা করা হয়।

একটি আণবিক সূত্রের সাহায্যে অনেক আণবিক সমাবয়বকে বোঝানো যায়। প্রত্যেক সমাবয়বই হল নিউক্লীয় কো-অর্ডিনেটেরফাংশন সম্পূর্ণ শক্তি (ইলেক্ট্রনীয় শক্তি ও নিউক্লীয় বিকর্ষণ শক্তির যোগফল) দ্বারা রচিত বিভবশক্তি-তলের স্থানীয় মিনিমা। যদি কোনো নির্দিষ্ট আণবিক গঠনের জন্য সমস্ত নিউক্লীয় সরণের সাপেক্ষে শক্তির অন্তরজ শূন্য হয় তবে তাকে স্থানীয় বিন্দু বলে। যে স্থানীয় বিন্দুতে যেকোনো সরণের ফলে শক্তি বৃদ্ধি পায় সেই স্থানীয় বিন্দু হল স্থানীয় মিনিমা। সর্বনিম্ন স্থানীয় মিনিমাকে সার্বজনীন মিনিমা বলে এবং এটি অণুর সবচেয়ে সুস্থিত সমাবয়ব। যদি কোনো একটি নির্দিষ্ট কো-অর্ডিনেটে সরণের ফলে উভয় দিকে শক্তি হ্রাস পায়, তাহলে এটি দশান্তর গঠন এবং সেই কো-অর্ডিনেটটি হল বিক্রিয়া-কোঅর্ডিনেট। স্থায়ী বিন্দু সন্ধানের এই প্রক্রিয়াকে জ্যামিতি অপ্টিমাইজেশন বলে।

সকল পারমাণবিক কো-অর্ডিনেটের সাপেক্ষে শক্তির প্রথম অন্তরজ হিসেবের জন্য দক্ষ পদ্ধতি আবিষ্কারের পরে জ্যামিতি-অপ্টিমাইজেশনের মাধ্যমে আণবিক গঠন নির্ণয়ের চল শুরু হয়। পর্যাবৃত্ত গতি হিসেবে নিলে দ্বিতীয় অন্তরজ থেকে আণবিক কম্পনের কম্পনাঙ্ক নির্ধারণ করা যায়।দ্বিতীয় অন্তরজ নিয়ে যে হেসিয়ান ম্যাট্রিক্স তৈরি হয় তার আইগেনভ্যালুর সাথে এই কম্পনাঙ্ক সম্পর্কিত। সব আইগেনভ্যালুর মান ধনাত্মক হলে সব কম্পনাঙ্ক বাস্তব সংখ্যা হয় এবং স্থানীয় বিন্দুটি স্থানীয় মিনিমা হয়। যদি একটি আইগেনভ্যালুর মান ঋণাত্মক হয় (কম্পনাঙ্কের কাল্পনিক মান), তবে স্থায়ী বিন্দুটি একটি দশান্তর গঠন। একের বেশি আইগেনভ্যালুর মান ঋণাত্মক হলে স্থায়ী বিন্দুর প্রকৃতি আরো জটিল হয় এবং এ নিয়ে বেশি চর্চা নেই। এরকম কোনো বিন্দু পাওয়া গেলে স্থানীয় মিনিমা বা দশান্তর গঠনের সন্ধান এর থেকে দূরে গিয়ে করতে হবে।

আপেক্ষিকতা উপেক্ষা করে ও শ্রোডিঙার সমীকরণ সরলীকরণের জন্য বর্ন-ওপেনহাইমার আসন্নায়নের মাধ্যমে ইলেক্ট্রনীয় ও নিউক্লীয় চলণ বিচ্ছিন্ন করার পরে যে সময়-অনির্ভর শ্রোডিঙার সমীকরণ থাকে তার আসন্ন সমাধান থেকে মোট শক্তি হিসাব করা হয়। এই প্রক্রিয়ায় মোট শক্তিকে স্থায়ী নিউক্লীয় কো-অর্ডিনেটে ইলেক্ট্রনীয় শক্তি, এবং নিউক্লীয় বিকর্ষণ শক্তি- এই দুইয়ের যোগফল হিসাবে লেখা যায়। এক উল্লেখযোগ্য ব্যতিক্রম হল প্রত্যক্ষ কোয়াণ্টাম রসায়ন, যাতে নিউক্লীয়াস ও ইলেক্ট্রনকে অভিন্নভাবে দেখা হয়। ডেন্সিটি-ফাংশনাল পদ্ধতি ও সেমি-এম্পিরিকাল পদ্ধতি হল এক-ই মূল ভাবনার ভিন্ন ভিন্ন রূপ। খুব বড় সিস্টেমের ক্ষেত্রে আণবিক বলবিজ্ঞান দ্বারা মোট শক্তি নির্ণয় করা যায়।আণবিক গঠন নির্ধারণের জন্য মোট শক্তি হিসাব করার কিছু পদ্ধতি হলঃ

এব-ইনিশিও পদ্ধতি

গাণিতিক রসায়নে যেসব প্রোগ্রাম ব্যবহার হয় সেগুলো সবই আণবিক হ্যামিইল্টোনিয়ান থেকে যে আণবিক শ্রোডিঙার সমীকরণ আসে তার সমাধানের জন্য নানা কোয়াণ্টাম-রাসায়নিক পদ্ধতি। যেসব পদ্ধতিতে কোনো এম্পিরিকাল বা সেমি-এম্পিরিকাল প্যারামিটার নেই এবং যেগুলি কোনো পরীক্ষামূলক তথ্য ছাড়াই শুধুমাত্র তাত্ত্বিক সূত্রগুলি থেকে উতপন্ন হয়েছে সেগুলিকে বলা হয় এব-ইনিশিও পদ্ধতি। যদিও এর থেকে এটা বোঝায় না যে এই সমাধানগুলি অভ্রান্ত, কারণ সবকটিই আসন্ন কোয়াণ্টাম-বলবৈজ্ঞানিক গণনা। এর থেকে শুধু এটাই বোঝায় যে যে আসন্নায়নটি ব্যবহার হয়েছে সেটি কোয়াণ্টাম তত্ত্বের প্রথম সূত্রগুলির পরিপ্রেক্ষিতে যথাযথরূপে সংজ্ঞায়িত এবং সর্বোচ্চ ত্রুটির পরিমাণ গুণগতভাবে জানা। যদি গাণিতিক পুনরাবৃত্তি ব্যবহার হয় তবে যন্ত্রের পূর্ণ ক্ষমতা (কম্পিউটারের সসীম "শব্দ-দৈর্ঘ্য" ও গাণিতিক বা ভৌত আসন্নায়ন মাথায় রেখে) ব্যবহার অব্দি পুনরাবৃত্তি করাই গণনার লক্ষ্য।

বিভিন্ন এব-ইনিশিও ইলেক্ট্রনিক গঠন গণনার আপেক্ষিক শক্তির চিত্র। দূরত্ব স্কেলে নয়।

সরলতম এব-ইনিশিও ইলেক্ট্রনিক গঠন গণনা হল হারট্রী-ফক (HF) পদ্ধতি, যা আণবিক অর্বিটাল তত্ত্বের সংযোজিত সংস্করন এবং যাতে যথার্থ আন্তর-ইলেক্ট্রনীয় বিকর্ষণের বদলে এর গড় মান নিয়ে গণনা করা হয়। ভিত্তি-সেটের আয়তন বৃদ্ধি করলে শক্তি ও তরঙ্গ-ফাংশন হারট্রী-ফক সীমা নামক সীমার দিকে অগ্রসর হয়। পোস্ট হারট্রী-ফক নামের গণনায় হারট্রী-ফক গণনার ফলাফলকে ভিত্তি করে আন্তর-ইলেক্ট্রনীয় বিকর্ষণের জন্য সংশোধন করা একে, যাকে ইলেক্ট্রনিক কোরিলেশন বলে। এই পদ্ধতিকে নিখুঁত করতে থাকলে এই সমাধান আপেক্ষিকতাহীন শ্রোডিঙার সমীকরণের অভ্রান্ত সমাধানের দিকে অগ্রসর হয়। পরীক্ষামূলক ফলাফলের সাথে সম্মতি আনতে হলে আপেক্ষিক ও স্পিন-অর্বিট অংশগুলিকে হিসাবে আনতে হবে, যা ভারি পরামাণুর ক্ষেত্রে গুরুত্বপূর্ণ। এই সমস্ত গণনাতেই পদ্ধতি নির্বাচনের সাথে সাথে ভিত্তি-সেটও বাছাই করতে লাগে। এটি হল কিছু ফাংশনের সেট যা অণুর বিভিন্ন পরমাণুর উপর কেন্দ্রীভূত। পারমাণবিক অর্বিটালের রৈখিক সংযোগ (LCAO) জাতীয় আণবিক অর্বিটাল গণনায় আণবিক অর্বিটালগুলিকে এই ভিত্তি সেটে সম্প্রসারণ করা হয়। তাই যেকোনো এব-ইনিশিও পদ্ধতিতে নির্দিষ্ট পদ্ধতিটি এবং ভিত্তি সেট নির্দেশ করতে লাগে।

হারট্রী-ফক তরঙ্গ-ফাংশন হল শুধু একটি ইলেক্ট্রনীয় বণ্টনের নির্ণায়ক। বন্ধন ভাঙার মত কিছু প্রক্রিয়ায় এটি অসম্পূর্ণ, এবং বিভিন্ন ইলেক্ট্রনীয় বণ্টন ব্যবহার করা দরকার। এখানে নির্ণায়কের গুণাঙ্ক এবং ভিত্তি ফাংশনের গুণাঙ্ক একসাথে অপ্টিমাইজ করা হয়।

আণবিক গঠনের ফাংশন হিসেবে মোট আণবিক শক্তি নির্ধারণ করা যায়, যা আসলে বিভবশক্তি-তল। বিক্রিয়া গতিবিজ্ঞানে এই তল ব্যবহার হয়। বিভিন্ন সমাবয়ব ও তাদের মধ্যে দশান্তর গঠন নির্ণয় করা যায় এই তলের স্থায়ী বিন্দুগুলি থেকে। যদিও এর জন্য সম্পুর্ণ বিভবতল গঠন করতে লাগে না।

একটি গুরত্বপূর্ণ লক্ষ্য হল গাণিতিক তাপগতিবিদ্যা, যা থেকে গঠন এনথালপির মত তাপরাসায়নিক রাশি রাসায়নিক মাপকাঠিতে অভ্রান্তভাবে হিসেব করা যায়। রাসায়নিকভাবে অভ্রান্ত বলতে বোঝায় যাতে ত্রুটির ফলে বাস্তব রাসায়নিক গণনায় গুণগত ভুল না হয়, এবং এর মান সাধারণত ১ কিলোক্যালরি/মোল বা ৪ কিলোজুল/মোল। ত্রুটিকে এই সীমায় বাস্তবে নিয়ে যেতে হলে অবশ্যই পোস্ট হারট্রী-ফক গণনা করে তার ফলাফল হিসাবে আনতে হবে। এই পদ্ধতিগুলিকে বলে কোয়াণ্টাম-রসায়ন যৌগিক পদ্ধতি।

রিচার্ড বাডের-এর অণুর মধ্যে পরমাণু (QTAIM) মডেল তৈরি হয়েছিল অণুর তরঙ্গ-ফাংশনের মত কোয়াণ্টাম ধারণার সাথে, অণুর মধ্যে পরমাণু,কার্যকরী মূলক,রাসায়নিক বন্ধন, লুইস জোড় তত্ত্ব, ভ্যালেন্স বন্ড মডেল ইত্যাদি রাসায়নিকভাবে কার্যকর ধারণার মেলবন্ধন ঘটানোর জন্য। বাডের দেখিয়েছিলেন যে এইসব পরীক্ষামূলকভাবে কার্যকর ধারণাগুলির সাথে যোগাযোগ আছে পর্যবেক্ষণযোগ্য তড়িত্ ঘনত্ব বিন্যাসের টপোগণিতের সাথে, যা পরীক্ষামূলকভাবে নির্ণয় করা যায় বা কোয়াণ্টাম তরঙ্গ-ফাংশন থেকে গণনা করা যায়। AIMAll সফটওয়্যার প্যাকেজে আণবিক তরঙ্গ-ফাংশনের QTAIM বিশ্লেষণ করা হয়েছে।

প্রচুর স্বনির্ভর গাণিতিক রসায়ন সফটওয়্যার প্যাকেজ পাওয়া যায়। তাদের কিছু কিছুতে অনেক রকম পদ্ধতি প্রয়োগ করা আছে, আবার কিছু খুবই স্বল্প, এমনকি একটা পদ্ধতিও প্রয়োগ করেছে। এগুলির বেশিরভাগই হলঃ

• জৈব অণু মডেলিং প্রোগ্রাম, যেমন প্রোটিন ও নিউক্লিক অ্যাসিড।
• আণবিক বলবিজ্ঞান প্রোগ্রাম
• কোয়াণ্টাম রসায়ন ও কঠিন অবস্থা ভৌতবিজ্ঞান সফটওয়্যার
• আণবিক ডিজাইন সফটওয়্যার
• সেমি-এম্পিরিকাল প্রোগ্রাম
• ভ্যালেন্স বন্ড প্রোগ্রাম

• জৈব তথ্যবিজ্ঞান

• C. J. Cramer Essentials of Computational Chemistry, John Wiley & Sons (2002).
• T. Clark A Handbook of Computational Chemistry, Wiley, New York (1985).
• R. Dronskowski Computational Chemistry of Solid State Materials, Wiley-VCH (2005).
• A.K. Hartmann, Practical Guide to Computer Simulations, World Scientific (2009)
• F. Jensen Introduction to Computational Chemistry, John Wiley & Sons (1999).
• K.I. Ramachandran, G Deepa and Krishnan Namboori. P.K. Computational Chemistry and Molecular Modeling Principles and applications Springer-Verlag GmbH আইএসবিএন ৯৭৮-৩-৫৪০-৭৭৩০২-৩.
• D. Rogers Computational Chemistry Using the PC, 3rd Edition, John Wiley & Sons (2003).
• P. v. R. Schleyer (Editor-in-Chief). Encyclopedia of Computational Chemistry. Wiley, 1998. আইএসবিএন ০-৪৭১-৯৬৫৮৮-X.
• D. Sherrill. Notes on Quantum Mechanics and Computational Chemistry.
• J. Simons An introduction to Theoretical Chemistry, Cambridge (2003) আইএসবিএন ৯৭৮-০-৫২১-৫৩০৪৭-৭.
• A. Szabo, N.S. Ostlund, Modern Quantum Chemistry, McGraw-Hill (1982).
• D. Young Computational Chemistry: A Practical Guide for Applying Techniques to Real World Problems, John Wiley & Sons (2001).
• D. Young's Introduction to Computational Chemistry.
• Lewars, Errol G. (২০১১)। Computational Chemistry। Heidelberg: Springer। doi:10.1007/978-90-481-3862-3। আইএসবিএন 978-90-481-3860-9।

• Reviews in Computational Chemistry
• Journal of Computational Chemistry
• Journal of Chemical Information and Modeling
• Journal of Computer-aided Molecular Design
• Journal of Chemical Information and Modeling
• Journal of Chemical Theory and Computation
• Computational and Theoretical Polymer Science
• Computational and Theoretical Chemistry
• Journal of Theoretical and Computational Chemistry
• Journal of Cheminformatics
• Journal of Computer Chemistry Japan
• Annual Reports in Computational Chemistry
• Computers & Chemical Engineering
• Journal of Chemical Software
• Molecular Informatics
• Journal of Computer Aided Chemistry
• Theoretical Chemistry Accounts

• NIST Computational Chemistry Comparison and Benchmark DataBase – প্রচুর সিস্টেমের গাণিতিক ও পরীক্ষামূলক ফলাফলের তথ্যভাণ্ডার।
• American Chemical Society Division of Computers in Chemistry – American Chemical Society Computers in Chemistry Division, resources for grants, awards, contacts and meetings.
• CSTB report Mathematical Research in Materials Science: Opportunities and Perspectives – CSTB Report
• 3.320 Atomistic Computer Modeling of Materials (SMA 5107) Free MIT Course
• Chem 4021/8021 Computational Chemistry Free University of Minnesota Course
• Technology Roadmap for Computational Chemistry
• Applications of molecular and materials modelling.
• Impact of Advances in Computing and Communications Technologies on Chemical Science and Technology CSTB Report
• MD and Computational Chemistry applications on GPUs